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Abstract  Report  1  Scalar  Scattering  from  Turbulent  Plasmas  by  a  Modified 


Method  of  Iteration 

The  scattering  of  scalar  waves  from  s  turbulent  plasma  is  investigated 
by  a  modification  of  the  iteration  procedure  used  in  deriving  the  so-called 
first  order  smoothing  equation.  Electron  density  fluctuations  are  modeled 
by  locally  homogeneous  turbulence. 

A  wave  number  space  representation  is  used  for  the  random  scattering 
integral  equation  which  is  decomposed  into  coupled  integral  equations  for 
coherent  and  fluctuating  scattering  amplitudes.  These  equations  are 
solved  by  invoking  the  "large  scatterer  approximation".  This  density 
independent  approximation  is  used  to  express  the  coherent  field  in  terms 
of  the  fluctuating  field.  Successive  iteration  then  yields  a  solution 
for  the  fluctuating  field.  The  Incoherent  radar  cross  section  is  calcu¬ 
lated  to  second  order  for  an  axially  symmetric  wake  section  Illuminated 
nose-on.  The  turbulence  is  modeled  by  an  isotropic  exponential  correlation 
function  and  a  turbulent  intensity  which  is  axially  uniform  with  a  Gaussian 
fall-off  radially.  The  present  model  predicts  a  cross  section  saturation 
with  increasing  turbulent  fluctuations. 

Abstract  Report  2  Vector  Scattering  from  a  Turbulent  Plasma  by  a  Modified 

Method  of  Iteration 

The  analysis  of  Report  1  is  generalized  to  the  scattering  of  electro¬ 
magnetic  waves.  For  vector  scattering  we  must  solve  two  coupled  vector 
integral  aquations  for  the  coherent  and  fluctuating  fields.  The  "large- 
scatterer  approximation"  (scatter  dimensions  large  compared  to  correlation 


length)  leads  to  a  vector  Integral  equation  for  the  coherent  field  In  terms 
of  the  fluctuating  field  involving  a  non-symmetric  tensor  whose  elements 
depend  on  scatterer  characteristics  and  the  dyadic  free  space  Green's 
function.  This  is  used  to  obtain  an  integral  equation  for  the  fluctuating 
field  which  is  solved  by  iteration  to  second  order  in  the  density.  Using 
the  quasi-normal  hypothesis  a  general  expression  is  derived  for  the 
Incoherent  bistatic  scattering  cross  section  of  an  axi-symmetric  wake 
section  illuminated  at  non-normal  incidence  and  arbitrary  polarization. 

Abstract  Report  3  Application  of  the  Modified  Method  of  Iteration  for 

Scattering  to  SRI  Experiment 

The  general  expression  for  the  incoherent  scattering  cross  section 
derived  in  Report  2  is  specialized  to  compute  the  backseat  ter  cross  sec¬ 
tion  for  direct  and  cross  polarized  directions.  The  SRI  turbulent  flame 
is  illuminated  at  normal  incidence  with  direct  polarization  along  the 
flone  axis.  The  turbulent  intensity  is  modeled  by  a  Gaussian  radial 
and  axially  uniform  distribution  with  a  correlation  function  of  exponential 
form.  These  models  are  close  approximations  to  the  measured  SRI  scatterer 
characteristics.  Direct  and  cross  polarized  backscatter  cross  sections  are 
calculated  for  a  wide  range  of  RMS  electron  density  values  including  under- 
dense  and  near  critical  density  regions.  Agreement  between  experimental 
and  theoretical  absolute  cross  sections  fall  well  within  the  3  db  SRI 
measurement  uncertainty.  For  the  very  underdense  case  the  theoretical  model 
yields  results  in  agreement  with  first  Born  theory  for  direct  polarization 
and  second  Born  theory  (Ruffive-DeWolse)  for  cross  polarization.  As  critical 
density  is  approached  the  model  predicts  cross  section  saturation  observed 
experimentally  (SRI,  RCA  Montreal,  Bell  Labs). 


Abstract  Report  4  Wave  Scattering  from  a  Unidimensional  Layer  with  Strong 


Random  Irregularities  by  the  Method  of  Smoothing 

Under  study  is  the  multiple  scattering  of  a  plane  wave  normally  incident 
oa  a  layer  filled  with  a  strongly  fluctuating  plasma.  The  electron  density 
fluctuations  are  locally  homogeneous,  and  the  mean  dielectric  constant  of 
the  layer  may  be  appreciably  different  from  the  imbedding  medium.  The 
formulation  of  the  problem  is  based  on  the  Dyson-Smoothing  equations  in  the 
bilocal  approximation.  The  approximate  treatment  of  the  mean  dielectric 
constant  as  homogeneous  throughout  is  shown  to  be  generally  inconsistent. 
Exact  treatment  of  the  problem  is  developed  in  the  limit  of  small-scale 
fluctuations;  the  dependence  of  the  scattered  power  on  a  broad  range  of 
plasma  parameters  is  displayed  and  discussed  for  typical  casses  of  physical 
Interest. 

In  the  case  of  a  mean  dielectric  constant  throughout „  a  new  mathematical 
formulation  is  presented.  The  resulting  Integral  equation  is  singular  and 
belongs  to  the  standard  form  studied  in  the  literature.  This  formulation 
is  closer  to  reality  since  it  avoids  the  assumption  of  sharp  boundaries 
for  the  fluctuations. 

Abstract  Report  5  Scattering  Coefficient  of  One -dimensional  Plasmas  of 
Epstein-type  Profiles  with  Random  Irregularities 

The  problem  of  electromagnetic  scattering  by  Inhomogeneous  one- 
dimensional  plasmas  with  electron  density  of  Epstein-type  profiles, 
superimposed  with  random  irregularities  is  studied  by  means  of  a  pertur¬ 
bation  method.  An  exact  Green's  function  is  derived  for  the  problem,  and 
the  resulting  Integral  equation  is  solved  approximately  by  the  Neumann 


Iterative  technique.  Coherent  and  Incoherent  power  reflection  and  transmission 
coefficients  are  obtained  to  a  consistent  second  order  accuracy,  and  an  expres¬ 
sion  is  given  to  define  the  parameter  validity  region  of  the  solution.  In 
addition,  the  incoherent  reflection  coefficient  is  graphically  displayed  for 
an  interesting  parameter  range. 

Where  meaningful,  the  special  case  of  the  half space  of  plasma  is  compared 
to  an  exact  Monte  Carlo  computer  experiment  for  a  plasma  slab.  It  is  shown 
that  for  the  given  region  of  validity,  the  solution  compares  favorably  with 
experiment. 


REPORT  1 


Scalar  Scattering  from  Turbulent  Plasmas  by  a  Modified  Method  of  Iteration 


1.  Introduction 

In  this  rep’ort  we  investigate  the  scattering  of  electromagnetic  waves 
from  a  turbulent  plasma  by  a  modification  of  the  iteration  procedure  in  the 
method  of  smoothing.  The  integral  equation  for  the  scattered  field  is  trans¬ 
formed  into  an  integral  equation  for  the  scattering  amplitude.  Coupled  in¬ 
tegral  equations  for  the  coherent  and  fluctuating  amplitudes  are  formed  and 
solved  by  a  method  which  invokes  the  "large  scatterer  approximation"  during 
the  iteration. 

2.  Analysis 

We  consider  the  scalar  scattering  problem  defined  by  the  integral  equation, 
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where 


i|>(r)  is  the  total  field,  i p  (r)  is  the  scattered  field 
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the  normalized  plasma  frequency  based  on  a  reference 
RMS  electron  density 
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n(ir, )  ■  - electron  density  fluctuation  normalized  on  n_..„. 

1  nRMS  RMS 


V  is  the  scattering  volume. 

In  the  far  field  Eq.  (1)  tor  the  scattered 
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We  identify  the  bracket  term  as  the  scattering  amplitude  in  the  r  direction, 
F(kr)  ■  F(<).  We  derive  an  integral  equation  for  this  quantity  by  multiply¬ 
ing  (1)  thru  by  n(r)  exp(i<*r)  and  integrating  over  V. 
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where  we  have  used  the  spectral  representation 
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equation  (3)  becomes 
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This  is  an  integral  equation  for  the  scattering  function  I(k)  and  according 
to  (2)  the  scattering  amplitude  is 


F(k) 


k2n2 
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(8) 


The  quantity  n(<)  is  a  spectral  representation  for  the  electron  density 
fluctuation.  We  assume  that  scatterer  dimensions  are  large  compared  to  the 
fluctuation  scale  (correlation  length  £).  The  scatterer  is  modeled  by  an 
electron  density  fluctuation  super-imposed  on  a  mean  electron  density  back¬ 
ground.  The  fluctuations  are  described  by  a  locally  homogeneous  turbulence 


-  3  - 


whose  intensity  is  proportional  to  the  square  of  the  mean  density.  At  the 
outset  (eq.  (1))  we  have  assumed  that  the  mean  background  contribution  to 
the  scattering  is  negligible.  If  in  Eq.  (7)  we  measure  wavelength  in  units 
of  correlation  length  then  all  quau-xties  can  be  interpreted  as  non-dimensional 
and  we  indentify 
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as  the  natural  expansion  parameter  for  the  problem. 

Next  we  form  the  equations  for  the  coherent  (Iq(k))  and  fluctuating 
(I^(ie))  scattering  functions  by  applying  the  ensemble  averaging  operator 
P  to  Eq.  (7) . 


PICk)  -  Io(k)  -  Pn(x-k)  +  t  G(p)  Pn(ie-p)  I(p) 


Iq((c)  ■  e  G(p)  Pn(x-p)  I^p)  (10) 

where  we  used  Pq  ■  0,  I  ■  I  +  I,  ,  PI,  *  0, 
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>00  — 

G(p)  =  ^2 — 2  *  and  vector  arrows  were  suppressed, 

k  -p 


If  Eq.  (10)  is  subtracted  from  (7)  we  obtain 

^(•0  -  n(K-k)  +  e  G(p)  n(<-p)  I0(p)  +  e  G(p)  n(ic-p)  ^(p) 

-  e  G(p)  P  n(ic-p)  I1  (p)  (11) 

We  shall  now  solve  pair  (10)  and  (11)  by  iteration  with  c  as  an  expansion 
parameter.  First  use  (11)  to  form  P  n(<-p)  I^(p)  for  use  in  (10). 

P  I^p)  h(k-p)  ■  P  n(p-k)n(<-p)  +  e  G(p)  P  n(p-Pj)nOe-p)  Iq(p1) 

+  e  G(p)  P  n(p-p1)n(>c-p)  i1(p)  (12) 

In  a  previous  study^  it  has  been  shown  that 

prp2 

P  n(P1)  n(p2)  ■  A(p1+p2)  ♦(-■• y- ■■) 

pl+p2 

P  n(p1)  n*(p2)  *  a(pi-p2)  . *  (-■  2~)  (13) 

1  J.Jarem, "Studies  in  Electromagnetic  Scattering  from  Turbulent  Wakes"  Drexel 
Univ.,ProJ.  349  Report, Sept.  1969,  Contract  DAH-CO-4-67C-0072JARPA  Order  1009. 


of  the  turbulent  in- 
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where  A(p)  and  (p)  are  the  Fourier  transforms 
tensity  and  correlation  functions  respectively.  The  scale  over  which  the 
Intensity  changes  is  large  compared  to  the  correlation  length.  In  the  limit 
of  a  large  scatterer  A(p)  approaches  a  delta  function.  We  shall  use  this 
property  in  evaluating  certain  integrals  in  the  iteration.  When  the  average? 
in  (12)  are  expressed  according  to  (13)  we  obtain 

,x+k 


PI1(p)n(K-p)  -  A(K-k)$(~2 - p) 

<+p 

+  e  G(p1)IQ(p1)A(ic-p1)'f (— 2 - p) 


+  e  G(p1)Pn(ic-p)n(p-p1)I1(p1) 

Viewing  A(x-p^)  as  a  delta  function  centered  about  p^ 
term  is  integrated  approximately  to 
e  Xq( <)  $(k-p)  G(p:)  A(x-p1) 

When  the  above  results  are  used  in  (10)  we  obtain 

Iq(k)  *  e  G(p)  A(x-k)  $(— ^  -  p) 
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we  can  solve  for  I  (< )  to  obtain 
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I  (K)  -  - -  {  e  A(K-k)  G(p)  *  (~  -  p) 

l-E^g1(K)g2(K) 

+  e2  G(p)  G(p1)  P  n(<-p)n(p-P1)  ^(p^) 


(19) 
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Next  we  use  (19)  to  eliminate  the  coherent  scattering  function  I  from  (11). 

o 

,  A(p-k)  G(p.)  $  (£~  -  p.) 
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3  G(p  )  G(p2)  P  n(p-p.)  n(p,-P2)  I,(P2) 

&  c  G(p)  n(ic-p)  - - 5 - 

1  -  c  gx (p)  g2(p) 


+  e  G(p)  n(ic-p)  I^(p)  -  e  G(p)  P  n(<-p)  Ij^tp)  (20) 

The  first  and  second  RHS  terms  can  be  combined  after  approximate  integration 
to 

n(tc-h)  +  t2  ,G(p) —  A(p.k)  G(p  )  $(k-p  )  -  n(x-k) - ^ -  (21) 

l-tg1(k)g2(k)  1  -  c  gx(k)g2(k) 

and  (20)  can  be  written 
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To  order  e  we  can  neglect  the  last  term  in  (22)  which  contributes  at  order 

3 

e  .  A  single  iteration  of  the  remaining  terms  yields 
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The  incoherent  differential  scattering  cross  section  is  related  to  the 


scattering  amplitude  by 
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Making  use  of  (23)  and  Pn  -  Pn 
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0  the  cross  section  can  be  written 
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+  c2G(p)G*(p1)Pn(ic-p)n(p-k)n*(<-p1)n*(p1-k) 


-  c2C(p)  Pn(»c-p)n(p-k)G*(p1)Pn*(>c-p1)n*(p1-k)) 

For  quasi-normal  turbulence  the  fourth  order  moment  in  (25)  can  be  expressed 
as  the  following  sum: 

Pn(«-p)n(p-k)  Pn*(x-  p1)n*(p1-k)+Pn(«-p)n*(ic-  Pj)Pn(p-k)n*(Pj-k) 


+  PnOc-p)n*(Pj-k)  P  n(p-k)n*(x-P1)  (26) 

We  note  that  the  first  term  in  (26)  when  used  in  (25)  cancels  the  last  term 
and  using  (13)  to  reduce  the  remaining  terms  leads  to 
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If  A  is  again  viewed  as  a  delta  function  for  the  integration  we  shall 
finally  obtain,  with  A(p)  and  $(p)  even  and  real, 


(27) 
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A(0)  k*  IT. 

o(k)  - - J - J  (♦(•e-k) 

(4ir)  1 1-c  glg2|  +  {,2  G(p)  *(p_k)  G*(pi)  A(Pl-p) 

+  c2  G(p)  $(p-»0  $(p-k)  G*(p2)  A(p1-(ie+k-p))  J  (28) 

3.  Backscattering  Example 

For  our  first  example  we  shall  compute  the  backacatter  cross  section  of 
an  axially  symmetric  wake  section  illuminated  nose-on.  The  turbulence  is 

R 

modeled  by  an  isotropic  exponential  correlation  function  exp(-  j)  with  cor¬ 
respond  ir-.t’,  spectrum 
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and  a  turbulent  intensity  which  is  uniform  axially  over  a  length  2b  and 
has  a  Gaussian  fall-off  radially  with  an  e-folding  radius  a//2*.  For  such 


a  model 


A(p)  ■  Aq  e 
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-  sin  bp 
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To  carry  out  the  integrals  called  for  in  (28)  we  choose  the  following 
corrdinate  system:  ^ 
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The  (x,z)  plane  is  defined  by  vectors  (k,z)  with  k  ■  i  cosa  +  z  sina  where 

*  AAA 

a  ■  angle  of  Incidence  and  z  falls  along  the  wake  axis.  (i,j,z)  form  a 

A  A 

rectangular  base.  For  nose  illumination  a  ■  90*  and  k  ■  z.  For  backscatter 
k  »  kr  ■  k(-k)  ■  -kz.  From  (17) 
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In  spherical  corrdinates  (p,8,$)  with  )c*p  ■  kp  cose  we  have  for  isotropic 


turbulence 
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The  integration  over  p  can  be  expressed  as 
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Using  model  (29)  in  the  first  RHS  tern  in  (31)  we  find  by  contour 


integration 


1  f  P  dp _  0  „  - o  J 

2ku  J  (p-2ku+ie  [1+(pt)2]2  k3  ^  ^ 

When  this  result  is  used  in  (30)  we  obtain 


'1  2n*o  u  +  fi((2u)2  - 


(  (2u) 2  +  -jV)2 

kV 


g ,00  -  ( 2it ) 3  i2  -2-~| 

1  1  +  (2ki) 


Next  we  compute  pj^kz)  according  to  (18). 

.  (n .  r  dp  a(°>  .  r  £  *<?> .  r 

j  p-(p+2lT)+ic  J  p.2T+lt  j  lk»3+l' 

With  model  (30)  for  A(p)  we  write 

1  2 

t-  p.(ap  )  2ir  r»  sin  bp. 

PpdPp  V  ds  dp3  bp3 

*  O 


sin  bpj  j 

bpj  2kpj  +  lc 


_  o  _ 


Plemelj's  formula  yields 

11.  r  .  *ln  bh  i _  i, 

e-*o  J  ap3  bpj  2kp^  +  ie  *  2k 

When  this  result  Is  used  In  (38)  we  obtain  after  Integration 


82(k) 


-iv(2ir)3  b 
— 


(39) 


(40) 


We  complete  the  evaluation  of  (28)  by  first  combining  the  last  two  terms. 

For  backscatter  (k  ■  -k)  we  have,  calling  this  sum  C2, 

C2  -  2c2G(p)  *(p+k)  *(p-k)  G*(px)  A(p1-p)  (41) 


For  model  (32)  the  above  integration  reduces  toJ 


,  ,  /If  *  '2* 

C2  -  2"3‘2  -**■ 


rir 


d* 


d0  $(pk+k)  0(pk"k) 


(42) 


For  the  present  turbulent  model  the  product  spectrum  is  written 


(l+k2i2|p+k|2)2  (l+k2i2|p-k|2)2 


_ ?! _ 

[l+(2ki)2)4[  1- (-(--— }--7) 2  cos20]2 
l+(2ki)Z 


(43) 


If  this  is  used  in  (42)  we  will  obtain 


C.  -  2V,2  /ST  CL  .  t6  2  *  4  <2^>2  *  .W>* 

1  T  [1  +  (2ki)fc]J  (1  +  2(2kOZ  r'2 


If  (44),  (40)  and  (36)  are  used  in  (25)  we  finally  obtain 

^3  o  _  /T  „3/„. A „2,  ,4 

d+8fc)fc  ~  *  u+n(i+2e‘) 


c  l  CTafc(kb)  ^2,2  Uh2  32  CT(ka) „2, 3/2  * 


.2  ir 


RC 


S<kb>  S3  .  .  .2  CT(kb>  S2  ,* 


i_  + 

5  8  iv  8 


1+6 


where 


_  «  2  “RMS 

5  -  0  ,  - 


Pi  n  • . 

c  r  1 1 . 


(44) 


8  ■  2kl,  t  ■  correlation  length 


|  ,  constant  of  proportionality  relating  m.s.  fluctuation 
to  square  of  mean  electron  density. 

■  length  of  wake  section 

■  effective  radius  of  wake,  e  folding  distance  in  Gaussian  radial 
fall-off. 

V 

■  Radar  cross  section 
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Vector  Scattering  from  a  Turbulent  Plasma  by  a  Modified  Method  of  Iteration 


1.  Introduction 

In  Report  1  we  investigated  the  scalar  scattering  of  electromagnetic 
waves  from  a  turbulent  plasma  by  a  modification  of  the  usual  iteration  pro¬ 
cedure  used  to  derive  the  integral  equation  for  the  coherent  field  in  the 
•o  called  method  of  smoothing.  In  this  report  we  generalize  the  modifica¬ 
tion  to  Include  vector  scattering  for  the  purpose  of  investigating  polariza¬ 
tion  effects. 

2.  Analysis 

We  start  with  the  vector  integral  equations  for  the  electric  field*- 


E(r)  -  t0<7>  +  k2^  |  dr,  ^(r-r^) *E(r^)  nU^t)  (1) 

—  —  *  —  i  V  •  r 

where  E  is  the  total  field;  E  *  q  e  is  the  ireident  field;  2  .  is 

o  pi 

^  „  r 

the  plasma  frequency  based  on  a  reference  RMS  electron  density:  G(r,rj) 

dyadic  Green's  function  for  free  space:  n(r,t)  electron  density  fluctuation 

normalized  on  n^g.  Integral  Eq.  (1)  is  transformed  to  an  integral  equation 

_ 1  _  * 

for  the  scattering  function  I(k)  ;k  •  kr, 


I(k)  ■  q  o(K-k)  + 


k2fi2. 

_ El 

(2ir)  3 


dp  (5(p)«i(p)  n (k-p) 


(2) 


where 

T(p)  »  [  dr  eip*r  n(r)  E(r)  (3) 

h 

n(»0  »  [  dr  eip*r  n(r)  (4) 

'V 


1  J.  Jarem,  "Studies  in  Electromagnetic  Scattering  from  Turbulent  Wakes". 
Drexel  Univ.  Project  349  Report.  Sept.  1969,  Contract  DAH-CO-4-67C-0072, 
ARPA  Order  1009. 


y 


G(p)  -  |  dR  eip‘*  G(R)  -  EP-kV--  +  -K 


2  2  '2 

p -It +ie  3kr 


(5) 


$(R) 


-ikR 
e _ 

-4nR 


,  1  i(R) 

3k2 


If  Che  integral  equation  can  be  solved  for  I(k)  then  the  scattering 
amplitude  F(*c)  is  given  by 


2  2 

kV 


F(k)  „  -  -•^1  $  -  rr)  •!(<) 


(6) 


(7) 


and  the  radar  and  differential  scattering  cross  sections  are 


respectively^ 


o  (s e)  ■  Air  o(k)  ■  4tt  <F(<)*F*(<)>  (8) 

Rv 

2  2 

k  n 

Tc  solve  (2)  we  consider  c  ■  - as  an  expansion  parameter  and  we  form 

(2it)J  _ 

the  coupled  integral  equations  for  the  coherent  (I^GO)  and  fluctuating 
(1^ (>c ) )  scattering  functions  in  the  usual  way  to  obtain: 


Ic(<)  - 

c  dp^pVP  I^p)  n(<-p) 

mm 

(9) 

\(7)  - 

q  n(«-k)  +  e  j"dp  £(p)«Io<p)  n(ic-p) 

(10) 

+ 

c  j"dp  <2(p)*(T1(p)-p  T^Cp))  n(x-p) 

Ue  have  used 

P  1(7)  -  <!(<)>  -  I  (7) ,  P  T.  (7)  -  0,  and 

0  1 

P  n(7)  ■  0.  If  we  let 

^(p)  - 

• 

dp  '£(p) 

then  (9)  and  (10)  can  be  written,  suppressing  arrows  over  vector  arguments, 

y<>  - 

c  ^G(p) *P  l^(p)  n(<-p) 

(ID 

17  p.M.  Morse  &  H.  Feshbach,  Methods  of  Theoretical  Physicsy  Part  II, 
McGraw-Hill,  1953. 


s 

-  xo- 


I^x)  •  q  n(*-k)  +  c  ^G(p)‘l0(p)  n(K-p)  ♦  c^(p)(l-P)*l1(p)  h(k-p)  02) 

We  use  expression  (12)  for  1^  in  Eq.  (11)  Co  evaluate  P  I^(p)  n(<-p).  Eq.  (11) 
becomes 

Tq(<)  »  c  fyp^q  P  n(<-p)  n(p-k) 

+  c2v<2(p)*<5(p1)*To(p1)  P  n(<-p)  n(p-p1) 

+  c2  ^(p)«^(p1)*P  n(«-p)  (I-P)  ^(Pj^  n(p-P1)  (13) 

The  correlations  can  be  expressed^as 

P  n(ic-p)  n(p-k)  »  A(x-k)  4(^  -  p) 

•c+p. 

P  n(x-p)  n(p-p1)  •  A(<-p1)  ♦(~2 - p) 

where  A(p)  and  *(p)  are  the  Fourier  transforms  of  the  turbulent  Intensity 
and  correlation  functions  respectively.  For  large  scatterer  A(p)  approaches 
a  delta  function  and  using  this  property  (13)  becomes 

IQ(x)  *  c  A(x-k)  *<$ (p )  •  q  «(*y^  -  p) 

+  e2  2(p)*^(p1)*T)(k)  A(x-p1)  *(x-p) 

+  e2  ^(p)*<S(p1)*P  n(x-p)  n(p-P1)  ^(Pj)  (14) 

The  symmetry  of  the  dyadic  operator  ^(p)  enables  us  to  write  the  second  RHS 
term  in  (14)  as 

e2  I0(ie)*<2(p1)  A(ie-p^)  *8(p)  «(«-p) 
and  the  equation  for  the  coherent  scattering  function  becomes 
T(<)-(d-c2^(Pl)  A(ie-p1)*?(p)  «(p-«))  •  ro(ic)*$(ie) 

«  c  A(x-k)  <*(p)*q  ♦(“Y^~'P) 

+  e2  ?(p)*4S(p1)*P  ^(Pj)  nGc-p)  n(p-pA) 


(15) 


where  the  T(ic)  tensor  defined  In  (15)  is  not  symmetric  in  general. 
Introduce  the  tensor  *5(0  such  that 
<fOc).‘Sto  -1 
In  terms  of  T(k)  we  have 


T21T32~T22T31  T31T12"T11T32  T11T22'T21T12  i 

L  -i 

where  T  ■  det(T(tc)). 

Multiply  each  side  of  (15)  from  the  right  by ''S(ie)  to  obtain 
I0(tc)  ■  e  A(<-k)  G(p)*q  $(^“  -  p)*S(<) 

+  c2  G(p)*G(p1)*P  ijiPj)  n(x-p)  n(p-p1),S(K)  (17) 

Use  this  expression  to  eliminate  Ifl(p)  in  (12) 

Ij^k)  -  q  n(<-k)  +  c2G(p)*n(<-p)  A(p-k)  H(p^)  »q  $(^~  -p^^p) 

+  c3G(p)*n(K-p)  B(p1)*^(p2)*p  I1(p2)  n(p-P1)  n(p1-p2)*^(p) 

+  e  £(p)  (I-P)*I1(p)  n(K-p)  (18) 


Ue  again  make  use  of  the  delta  function  like  property  of  A(p-k)  to 
simplify  the  second  RHS  term  in  (18).  This  term  is  approximately 


e2  n(ic-k)  G(p)  •  A(p-k)  Gfp^q  $ (k-p^  '^(k) 


Combining  this  term  with  the  first  RHS  term  and  letting 


Q(k)  =*  q  +  e2  1T(p)  A(p-k) *q.^T(pi  )$(pi -k)  'V(k) 


(19) 


we  can  write  (18)  as 


iiOO  -  Qn(ic-k)  +  e  G  (p)«  (<s-r)i1(p)n(<-p) 


(20) 


+  e3  7f(p)-"cf(pi  )‘^(p2)*n(K-p)Pn(p-Pi)n(pi-P2ni  (P2 )  •T’(p) 


We  now  use  an  iteration  technique  to  solve  (2q).  Let 


IlOO  *  Jq(<)  +  £  Ji(k)  +  c2J2(^)  +  +  0(e4) 


(21) 


Using  this  perturbation  expansion  in  (20)  and  comparing  terms  of  same  order 
in  e  we  obtain  the  following  iterative  series: 

JoOO  *  Qn(x-k) 


A  .  . 

Jl(K) 


J2(k) 


*^3  OO 


G  (p) • (o-P)n(K-p)Jp(p) 

-A 


^ (p)*  (6-P)n(ic-p)Ji(p) 

t 

(p).  (6-P)n(ic-p)J2(p)  +  g"*(p)  (pi)*"cT (P2)n(>c-p)Pn  (p-Pi)n(Pi~P2) 


(22) 


•Jq(P2) *  S  (p) 


The  scattering  amplitude  Fj(<)  to  order  e3  is 
(kfi  )2 

Fj (k )  ■  -  — -4"-  (5  -rr)»(Jp(<)  +  eJjOO  +  e2J2(<)  +  e 3 J3 (< ) ) 


(23) 


In  the  direct  (q)  polarization  direction  the  scattering  amplitude  is  given 
by 


F (< )  *  q* ?i (k) 


)2  »  *  ^ 

- Vi/ —  (q-q-rr)*!!  (tc)  =  d*  Ii  (<) 


In  the  cross  polarized  direction  (q^  ■  q  x  k)  the  scattering  amplitude 


is  written 


Fj_(<)  -  qj_*  Fi  (k  ) 


(kJ-p  )2  »  A  **  _ 

— ^ —  (qj_-q|Tr)*Ii(K)  5  c-Ii(k) 


To  form  the  cross  section  in  the  direct  polarization  direction  we  square 
and  averse  (24)  and  use  the  following  expressions  for  the  second  and 
fourth  order  moments:  (odd  moments  are  assumed  to  vanish) 


Pn(pi)n(n2)  =  A(p1+p2)<'(  ~—~2) 

Pn*(pi)n(p2)  =  A(p!-p2)v  (P-y--2) 

Pn(pi)n(p2)n(P3)n(P4)  *  Pn(Pi)p(P2)Pn(p3)n(p4> 


(26) 


+  Pn(pi)n(P3)Pn(p2)n(P4) 


+  Pn  (pi)n(P4)Pn(P2>n(P3) 
o(k)  *  FF(<)F*(ic)  *  P  j  d  •  I  x  (< )  |  2 

■  P  |  d*0n(<-k)  +  ed*^T(p)  •Q(n(K-p)n(p-k)-A(K“k)(&(^”- -p)) 

+  e2d-  (p)*’<(r(pi),Q(n(<-p)n(p-Pi)n(Pi-k)-n(ic-p)A(p-k)$(~-  -pj)) 


+  c3d«  (p)* g  (Pi)*  G  ^P2)*Q(n(<-p)n(p-Pi)n(pi-P2)n(P2-k) 

-  n(»e-p)n(p-pi)A(pi-k)«(^-—  -p2)  (27) 

-  Pn(*c-p)n(p-pi)n(pi-p2)n(p2-k) 

\ 

+  Pn(<-p)n(p-Pi)A(p1-k)<t(^^  -p2))  I2 
After  considerable  reduction  we  finally  obtain 
o(ic)  “  Ao$(ic-k)|  d*Q  |2 

+  ez2*(  -k)Real(d**Q*)  (d*^T (p)  .*0  (pj)  *Q) 

X$(tc-p)  [A2 (pj-k)  +  A2'T>i-p-k+*c)J  (28) 

+  e*(d.T(p)*0)(d.'G>(J'*q)*Q*) 

X$(tc-p)$(p-k)  [A2  (q-p)  +  A2(q-K-k+p)J 

Equation  (28)  is  the  direct  polarization  differential  cross  section  in  the 

A 

scattering  direction  k  ■  kr.  The  result  neglects  terns  of  order  c4  or  higher. 
The  corss  polarized  cross  section  is  readily  obtained  by  replacing  d  by  c 
in  (28). 
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APPLICATION  OF  THE  MODIFIED  METHOD  OF  ITERATION 
FOR  SCATTERING  TO  SRI  EXPERIMENT 

I.  Introduction 

- j- 

In  Report  2  the  modified  method  of  smoothing  was  applied  to  calculate 
the  direct  and  cross  polarization  cross  sections  for  scattering  from  a 
turbulent  plasma.  In  this  report  we  specialize  those  results  to  compute 
the  backscatter  cross  sections  for  both  polarizations  for  the  SRI  flame 
scattering  experiment  at  9.4  CHz.  The  scattering  geometry  is  shown  below: 


The  transmitter-receiver  is  located  about  50  cm  above  the  nozzle  and  the 

A 

flame  is  illuminated  broadside  with  direct  polarization  (q)  along  the  z 
axis.  The  SRI  jet  has  a  non-homogeneoua  electron  density  distribution  in 


which  the  turbulent  electron  density  fluctuations  are  characterized  by  a 


locally  homogeneous  turbulence  of  the  form 


<6n(r1,t)  fin(r2,t)>  ■  (n A(— =y^)  B(rj-r2) 


(1) 


where  n^  is  th$  rms  value  of  electron  density  at  the  reference  point  in  the 
flame  (z*0,  pa0),  A  is  the  normalized  turbulence  intensity  function,  and 
B(r^-r2)  is  the  normalized  turbulence  correlation  function.  The  turbulence 
intensity  is  modeled  by  a  Gaussian  radial  and  axially  uniform  distribution 
-¥■ 


A(R>  - 


Z(t) 


(2) 


where  R>pp+Zz> 


rl+r2 


a  *  3  to  5  cm 


Z(z)  ■  1,  -b<Z<b,  b  ^  5  to  10  cm 


•  0  otherwise 


For  the  correlation  function  we  assume  the  analytically  simple  isotropic 
exponential  function 


B(rl  *  r2*  *  • 


rrr2i 

t 


(3) 


where  the  correlation  length  l  *  1.75  to  2.5  cm.  Measurements  indicate 
that  the  radial  turbulence  scale  is  about  1.75  cm  and  the  longitudinal 
scale  is  2.5  cm.  The  spectrum  representing  B(r)  is 


♦00 


U  +~0ci)T? 


where  $  ■  8ni' 


(4) 


3 


We  will  also  need  the  spectrum  of  the  turbulence  intensity  distribution 


A(p) 


|  dR  A(R)e1|>'R 


4  *t"',3b  .  2. 

A  e  - r - ,  A  ■  ira  b 

o  p-b  *  o 


(5) 


2.  Cross  Sectio..  Functions  for  the  SRI  Model 

In  this  section  we  derive  the  explicit  form  of  the  functions  necessary 
to  compute  the  scattering  cross  sections  for  the  SRI  model.  In  (R2.17)  we 
need  the  vector 


Q(k)  -  q  +  e2  G(p)  A(p-k)*q*G(71)  $(Pj-k)  *S(k)  (6) 

/ 

where  J?(k)  is  defined  in  terms  of  the  elements  of  tensor  T(k).  The 
latter  is  given  by  (R2.1S) 

T(k)  ■  6  -  e2  G^)  A(k-p1)  *G(p)  O(p-k)  (7) 

Introducing  the  tensors 


and 


#00  —  ~ 

G(p)  #(p-k)  -  dp  $(p-k)  (P-£-  y--~-  +  -—•) 

j  p2-k2+ie  3k^ 


(8) 


f(k)  -  G(p)  A(p-k) 


(9) 


We  can  write  (7)  as 


T(k)  -  t  -  e2  f <lc) -g(lc) 


(10) 


We  proceed  to  the  evaluation  of  g  and  f  for  models  (4)  and  (5).  Let 


2  ("dp  *(p-k)  pip. 

’*i  "  ^  2  1 

p  -k*-lc 


(11) 


-  4  - 


In  spherical  coordinates  with 


palp  sin8  cos$  +  j  P  sin0  sin^  +  *  cos6 


2n  cos2< 


»u  -  k2  |"dp  p*  J  d*  j  ain6  d0  4>(  |  p-k| ) 

V 

Taking  the  polar  axis  along  k  we  have  for  model  (4) 


r  4  fn  d0  sin©  C0i 

‘u  •  7?  I.  "  "  lo  (pW-lcXp2 


+k2+t2-2kp  ccs0)2 


r  4  /' 

-  -f? :  dp  p 

2kV  i. 


d0  sln8  cos  8 


o(p2-k2-ic)(pZ-2kp  cos0  +  k“+Tfc) 


2 .  _2.2 


2  *2 
where  t  •  t 


Letting  u  ■  cose  we  can  write  (15)  as 


*11  ■  VkV 


fu2  d„  f-y-ja 

J ,  J  (P  -k  - 


PJL _ 

* — i  2  2 

ic)(p“-2kup+k'+t  ) 


'Calling  the  integral  over  p,  Ij^  we  write 

A,  B,  C, 


h  -  U"  fdp  [  p-(^U)  +  p-(-WO  Vca+IB)  )2  +  P‘(“+18) 

r-*n  *  ■  * 


€"►0 


U1  ,„  + _ X- 

+  ,2  p-(o-iS) 


(p-(a-iB)  ]‘ 


■] 


where  a  -  ku,  3  -  /t2  +  k2(l  -  “2>  »  °  and 

.*  e 

t  »  u2 


(12) 


(13) 


*  2k(k-ni>  (k-p2)‘ 


(pj-u2)  (p^r  pj-k' 


.  .  k2  +  P2 _ s 

2  (_2  ,.2  P1-P2)’ 


(14) 


(15) 


(16) 


(17) 


(18) 


P1  •  a  +  i$t  P2 


a  -  18. 


-  5  - 


Integral  (17)  can  be  evaluated  with  the  contour  C  ■  T  +  C  show  below: 


P,  -  a  -  IS 


i  dp  — E -  -  f  dp  (  )  +  f  dp  (  )  -  2tt i ( A,  +  C,) 

l  <  )<  >:  l  l  12 


We  have  since  the  contribution  over  C_  vanishes 

CP 


Ij  •  2nl(A^  +  C2)  and 


•ll  " 


du  2*1 


i  2(k-Pl)^(k-p2)fc 


*4o  fl 

slu 

.jd _ (-ii+A-,1 

(PJ-k2)(Pl-P2)2  pj-k2  pl"r2  J 


Calling  the  contribution  from  Aj,,  I we  have 


(19) 


(20) 


ltT  «  k  rl  2 
o 


kll 


2t 


/  u_  du _ 

[j  (k-a-iS)2(k- 


a+18) ' 


(21) 


-  6  - 


Carrying  out  this  integration  and  letting  1  ■  kl  we  obtain 

o 


-i"  ♦  k  r  2  2  2  d+212)2  *1 

I„  - - f—  (l+21‘)ln(l+41‘)  -  2r<l  + - f--) 

11  8l6  1-  °  °  0  1+4 12  J 


(22) 


With  1^2  (he  contribution  over  integral  C2  we  write 


*12 


-n‘*  2i 
0 

kV 


-jiliE-pi _ +  Jl—) 

(p2-k2)(p1-p2)2  p2“k2  Pl“P2 


(23) 


It  we  let  x  •  p.  «  a  +  16  »  ku  +  i  /k2+T2-k2u2 
I  i 

then  Integral  (23)  is  transformed  to 


-it2* 


X*2  4hV 


k+it 


2  2  2  2  2  2  2 
xdx(xW+xV  k  k>T  . 

.  2  .2W  2  ,2  2*  ^  2.2  2  ,  2  2; 
(x  -k  ) (x  -k  -T  )  x  -k  x  -k  -T 


-k+it 

and  integration  by  partial  fractions  yields 


(24) 


-in2k* 


X12- 


81 


1-121 


.4 


<1+21‘)  In  +  141  (l+l2  +  -^—)] 

o  l+i2lo  o  o  1+4£2  J 


(25) 


Combining  (22)  and  (25)  we  obtain 


*11  ’ 


■in2*  k 
o 

81? 


f(l+2i2)  ln(l+4l2)  - 
L  o  o 


2  <l+2*o)2 

2t2u  + - f— ) 

0  1+412 


2  1-i2lo  2  1 

+  (1+20  +  141  (1  +  r  +  -5-r) 

o  W*21o  o  o  l+4l2J 


(26) 


-  7  - 


Next  we  calculate  a^  according  to 

2ff  fit 


*22 


k‘2  f dp  p2  f"p  fde  sinQ  «(|p-k]) 

*  o  'o  Jo  p  -k-ic 


Jk2.* 


f« >p‘{ 


ffd0  sin8  sin20 


2  2  2  2  2  2 
(p  -k  -ie)  (p  -2kp  cos0  +  t  ; 


(27) 


with  u  *  cos 9  equation  (27)  becomes 


■•o  -i  2.  f-  dp  i 

22  *  “TT  du'1‘u  >  ?  2 

22  2k  V  J_L  (p  -k  - 


le)(p2-2kup+k2+T2)2 


(26) 


The  above  Integral  can  be  broken  up  into  the  following  Integrals: 


iff  #  k  fl 
+ - - du - . 

2l4  Jx  (k-Pl)2(k-p2)‘ 


a22  ’  "  ®11 


JsL  f1  _du.  Pi _ (_JsL  +  JjL-) 

t*  i  (pi-k2HPl-p2)2  (p2rk-’ 


2i*V  rl 

“72 


(29) 


Calling  the  middle  term  I2  we  have 
rl 


iff2k# 


X21" 


21 


du 


iff2k* 


(2k2+T2-2k2u)2 


21 


T2(T2+4k2) 


-1 


(30) 


with  1  ■  ki  this  becomes 

o 


iff2  k  ♦ 


X21- 


1+41 


(31) 


Calling  the  last  term  in  (29)  I22  and  letting  x  ■  px 


-2iff2* 
_ 0 

1 

fk+lt  d,  K3 

1  *»"  4 
kre 

2k 

(x2-k2)(x2-k2-T2) 

-k+iT 

2  2  2 
V  k  +“ 

( - - —  +  *  1  —  -  ) 

'  1  9  T  1  1  1' 


(32) 


-  8  - 


Carrying  out  this  integration  we  finally  obtain 


-i2ir2*  .  2. -.2 

o  i  t_+5k_ 

.2.4  4t  "2.  ..2 

hi  T  +4k 


I  ■ 

x22 


ir2k*  1+512 

T  .  O  _ _ 0_ 

X22  3  2 

21  l+41x 

0  0 

v 


Using  (31)  and  (33)  in  (29)  there  results 


iir2k$  ir2k$  1+512 

*22  ■  "*11  +  - f  + - r - 2~ 

1  11  1+41 Z  21J  1+4  fc 

00  o 

Combining  terms  in  (26)  a^  simplified  to 


"i*  k*0  r 


1  (l+4i2+2l4)-i(l+5i2+5i4)’ 

o  0  o  _  0  o 


- - <1+2C)  in(l-i2l  )  -  21  - 

11  4i6  L  0  00  1  +  UJ2 


When  this  is  used  in  (29)  there  results 
l»lMn  .  , 

•22  - -JT-  [  <1+2to>  -  2W 


Next  we  evaluate  a.,  according  to 


-  k  2  [  dp  p2  [  d$  [  sin9  dO  «*» ( |  p’-k |  )  P\—Cy~— 
0  *  0  ' 0  p  -k-ic 


2 _ 2, 


Comparing  this  Integral  with  (27)  we  observe  after  integration  over  $  that 


*33  ■ 2  *u 


Since  a^(l^j)  involve  sin$  or  cos$  Integrations  over  (o,2ir)  these  com¬ 
ponents  vanish,  or, 

»ij  ■  °»  i  ^  J  (3 

The  next  type  of  integral  to  evaluate  £  in  (6)  is 


am  • 


f  dp _ 

*  (p*-k*-ic) (t  +|p-k  *) 


-9- 


This  Integral  is  evaluated  in  Morse  and  Feshbach  (pp.  1083). 
-w2k$  (1  +  i2i.  ) 

-u-  Tt77?T 

o  o 

The  last  integral  in  (8),  call  it  c^,  is 


'ij 


♦>  rcm 

3k 


<S>(  I  p-k| )  - 


«->> 
A  it  6 


3k 


dp  p2  4*(  |  p  |  ) 


c 


ij 


An  $  6 


f"  X2  dx 
(x2+t 


2)2 


n2k  4> 


3 1' 


ij 


(Al) 


(A2) 


(A3) 


Combining  (35),  (Al)  and  (A3)  we  obtain  ?u. 


«11  "  a*ll  +  bll  +  C11 


n2k$ 


II _ i-f(l«l2)  tn(l+12l  )  - 

uib  L  °  0  i+Ar  J 


2  3 

l+r+i6r 

o  o 

3i3(l+A£2) 
o  o 


(AA) 


similarity  (36),  (Al)  and  (A3)  yield 


822 

2 

n*k$ 


*  * 

a22  +  b22 


+  c 


22 


-i(l+2i")  1  (l+i2i  ) 
o  n  o 

«i6 


(l-ll  )  l«2+i6t’ 

0  .  O _ o 


21: 


3i2(l+Ai2) 
o  o 


(45) 


-10- 


Final  ly  (38),  (41)  and  (43)  yield 


*33 


ir2k* 


(l+’r)  in  (1+2$,  ) 
o  o 


1  +  l2  +  i6i3 
_ o _ o 

3i3(l  +  4i2) 
o  o 


-  (i  (l+4«2+2iJ]  +  i[l+5£2+5£*))] 

1+4  a 2  0  0  0  O  O  J 


(46) 


Next  we  evaluate  f ^ (k)  according  to  (9) . 


"  {  <iP  A(p-H  [  -  -PA-PJ.  .{.1J  + 

13  J  L  p  -k  -ie  3k  J 


(47) 


For  model  (5)  the  last  integral  in  (47)  is 


°ii  r  -  -  -  6ii  r 

— 2  dp  A(p_k>  “  — 2  pft  dpr 

3k2  3k2  P  C 

ft 


1 2ir  r* 

d<j> 


2  2 
po  a 

- -  Sin  bp 

dp3  *„«  3 


bP' 


(2ir) 3  6 


£1 


3k 


(48) 


The  first  integral  in  (47)  is  evaluated  approximately  treating  A(p-k)  as 
a  quasi  delta  function. 


f  dp  A(p-k)  [p^  k~2-6±^  ] 


2  I2  4 

p  -k  -ie 


(kikj-  *<j) 


r  -  -- 

dp  A(p-k) 
2Tc*  (p-k)-ie 


(49) 


-11- 


For  model  (5)  the  last  integral  in  (49)  is 

_  _  .co  .  \ 

dp  A(p)  u  dpidp2dp3A(je  a  *¥  sin  bp3 

2k*p-ie  2kpicosa  +2kp3sina-ie  bp3 
•  00 

Integrating  over  P3  and  recognizing  that  the  principal  value 
vanishes  we  obtain 


dp  A(q)  _  iirAo 

* 

2lc*^-ie  2k  sina 

dp2 

dPle-f-(^2> 


sin(bpicota) 

bpicota 


with 


>O0 

dP2e 

-OO 


2/2 7 


0 

equation  (51)  reduces  to,  x  1  bp  cotcx  ,  q  * - tana 

1  /8b 


Itt i/2tt  Aa  2tana 
ka  sina  b 


.  siux  -q2x2 
dx - e  M 


But  it  can  be  shown  that 


sinx  -q2x2 
dx - e 


I  erf<2?> 


Using  this  result  in  (53)  and  recalling  that  Aq  ■  ira2b 


-  i^~(2Tr)3kb 
4k2sina 


,  A  b  x 
ecf  a~  Tana 

(  &  -  ) 
a  tana 


(50) 

contribution 


(51) 


(52) 


(53) 


(54) 


(55) 


-12- 


We  note  that 


„  —  and  erf(Z)-— ->  1. 

2  Z  +  0  /;  Z-- 


Using  (48),  (49),  and  (55)  in  (47)  we  obtain 


7( k)  -  4-f-  +  (7-kk)f(a)} 

t,2  •> 


(56) 


where  we  wrote 


f(a)  - 


i/H  kb  er  a  tana 


4  sina 


a  tana 


(57) 


Expanding  (56)  we  write 


f  00 


(2tt  ) ' 


v  +  sin2af(a) 


-  sinacosaf(a) 


3+  f(a) 


0 


sinacosaf (a) 


+  cos2af(a) 


(58) 


In  terms  of  g  and  f  the  tensor  (7)  becomes 


T  00- 


1  -  e2f 


ll8ll 


1  -  e2f 


22g22 


f13gll 


f13g33 


1  -  c2f 


33833 


(59) 


13- 


T  -  det(T)  -  TuT22T33  -  Ti3T22T31 


(60) 


(1"^  ^22^22^  (1— ^  fnRn)(l-e 


11*11' 33h33  *11*33  13 


The  elements  of  S  (Report  2,  eq.  (16))  are 


1 

*S  "  T 


(l_e 2 ^22^*22^  ^1_e2f33R33^ 


0  *^13^33^*c^22^22^ 


0  (1“e2fll*ll)(1“c2f33833)‘fi38llR33  0 


'f13Rll(1'e"f33S22) 


0  (l“c^f 22822^  ^”c2^11811^ 


(61) 


For  the  STl  experiment  the  flame  is  illuninatcd  broadsides  so  that  a  “  90* 
and  f1j(«  ■  0)  *  0.  Also,  from  (57) 


t(o) 


i/flag 

4/2 


so  that 


1 

3 


7  .  (2»)3 

k2 


0 


0 


0 


1^  ,  i^ao 
3  4/2 


0 


(62) 


0 


—  + 

3  17F 


-14- 


.IM 


3  .. 


s .n,  ...  .  ;:wi  .  ikbx  .  **  is 

2  ^  JJ)(j  ^  2  ^  +  **  3^ 


(62) 


The  elemenCs  of  simplify  to 


!  1  “ c  fii8n 


su- 


1  “C  f22g22 


i  -  C  £33833 


ii 


l“*  fll8ll 


J± 


zz 


l-c‘f„g 


22*22 


1— C  £33833 


(63) 


When  (62)  and  (63)  are  uped  in  (6)  we  obtain  for  Q  the  following  result 


for  q  •  z: 

'i  •>  „*  •  '  •  - 

Q(k)  -  q  +  c2(U  fn  +  jj  f22  +  zz  f„)»q(ii  gn  +  Jj  g„  +  zz  g„)*S 


33 


'll  T  *22 


*33 


1  “  e  £33833 


-  Z  Q 


(64) 


3.  Direct  Polarization  Backscatter  Cross  Section 

In  this  section  we  use  the  general  result  of  Report  2  (Eq.  27)  to 
calculate  the  backscatter  differential  cross  section  for  the  direct 

A  A 

polarization  case.  He  recall  that  q  »  z  (direct  polarization  direction), 

*  a 

a*0  (normal  Illumination),  k  ■  i  (direction  of  incident  illumination), 

A 

ic  ■  kr  ■  -k  (backseat tering) .  For  backscattering  and  direct  polarization 


we  have 


-13- 


k2P.2  . 

5 - 4^  Z 


Q  •  z  Q  ■  — 

1  f33R33 

k2h2 

d  •  o'  •  —  --P-1  Q 

4* 


k2i>  o 


d*  C  (?)•  G  (pi)*Q  *  ”  5k^4n‘ 


dp [3  ("• p)p+Z (p2-4k2  )  ) 
p2-k2+ie 


f  tlPl  ^(Z-p'x  ^x^Cpj-Ak2)] 

pf-k2+ie 


k2fi2  Q 

d*^T(p)*o  * - 

3k24ir 


dp{3(Z*jp)2  +  P2"4k^]_ 

p2-k2+ie 


When  these  results  arc  used  In  (23)  the  direct  polarization  backscatter  cross 
section  becomes 


o(-k) 


k4r>piA0i(2k)|o|2 

(4n)2 


•Je2 

{1  +  —  Real 

9k4A0 


dp_  i'iJjP )'  (p+k) 
p2-k2+ie 


do  A2 (q+lQ 
q2-k2+ie 


+  -iiL  Real 
9k4A0 


dpjijjn )$(p+k) 
p2-k2+ie 


dq  A2(q-p-2k) 
q2-k2+ie 


(65) 


+ 


_2ef_ 

9k% 


^(Zk) 


dp*  h3(p)fl(p+k)$(P’k) 
p2-k2+ic 


‘dq  A2  (q  ~P  )  j 

q2-k2-ie 


-16- 


Call  Che  second,  third  and  fourth  terms  in  the  brackets  of  (65)  01,  02,  and 
03,  respectively,  and  where  the  polarization  dependent  functions  are  defined 
by 

hi(p)  -  [3(Z-V)p+Z(p2-4k?)]-{3(Z.k)k+Z(k2-4k2)] 

(66) 

-  3k2(4k2-p2-3(Z.i)2) 

h2(p)  -  (3(Z*p)p+Z(p2-4k2) J  •  [3  Z«  (p+2k)  (p+2lT)+Z(p+2k)2-Z4k2] 

(67) 

■'  (p2*41  ’)  (p2+4k*p)+3(Z*  p)2  (5p2+lGk .p-4k2) 
h3(p)  -  (3(Z*p)2+p2-4k2)2  (63) 


Next  we  proceed  to  the  evaluation  of  trj. 


2  c2  „  , 
0 1  - - -  Real 

9k4A0 


dp  3k2(-p2-3(Z-p)+4k2)  f  _ dg  A2 (a) 

p2-k2+ie  J  (q-k)2-k2+ie 


(69) 


- -  ‘  - 


We  evaluate  the  integral  over  q  (call  it  o12)  approximately. 


°12 


dqA2(q)  . 

d<L  A2Xq).  .  h2 

-  —  a2 
dq ie  4  q  1 

_  iisi 

dq2e  4 

1 

dq3 

sin  bq  3 

q2-2k.q+ic 

* 

-2k*q+ic 

. 

-2kqi+ie 

b93 

in  IT  9,  X9  in 3 /tT  . 

"oir  — —  r  (» a  b)^  -  -  — g—  Ao 


2k  a  b 


(70) 


17- 


The  integral  over  p(on)  is  evaluated  next.  Take  i  as  the  polar  axis  with 
k*p  ■  pcosO  ,  Z*p  *  psin0sir>$,  u  ■  cos9. 


2ti 

d<P  3k2  (-3p2sin20sin2i{>-p2+4k2) 

0 

(1  +  i^(p^+2kpcos3+U*)P 


(71) 


01  "  "  RJX  8  a°*'0  [1*+  2?§(1+u)]2"  *  ^Pl  8*aolo™l(lo)  (73) 

-1 

The  last  integral  is  evaluated  to  yield 


[(£q)  - 


i  +  <»-»* 

(1  +  2 1?*)7  dx 


4lo  +  41*  +  1 

^<TV4i2) 


-  ^  1"  (H4i§) 


-18- 


Next  ve  evaluate  02  approximately. 


o2 


2e2 

9k4A0 


Real 

1 


dp  h2(p)$(p+k) 
p2-k2+ie 


_ dq  A2(q) 

(q+p+2k)2-k2+ie 


(75) 


Evaluate  the  integral  over  q(<?22)  first.  Neglecting  small  q  terms  in  the 

denominator  we  have  approximately 


022 


' 

4 


dq  A2 (q) _ 

4k  •  <i+p2+4k  •  p+3k2+ie 


(76) 


With  4k *q 

■  4kqi  and  po 

P'2+4iT.  p+3k 

"  "4  k 

_  cj  2 

.2  2 

°22  " 

d«Ue"jr_lL 

-  a  92 
dq2e  4 

4k  , 

qi+po+ie 

we  write 


dq3 


sinhqj 

bq3 


ir2a/^Ao 

2k 


p  .  -°M 

}  -lire  4  +  P 


r  * 


a2  ..2 


d<lie _ 4_ 


qi+PO 


i 


(76) 


But  it  can  be  shown  that  the  principal  value  integral  is  expressible  in 
terms  of  the  Dawson  Integral 


P 


dqie 


91+PO 


2/?D(^) 


(77) 


022 


*^An 


2k 


[  _llTe 


(78) 
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When  this  result  is  used  into  (75)  and  we  introduce  spherical  coordinates 
we  have  after  i  integrations 


o2 


2  e2 
9k4A0 


"cal 

V 


1 

2 


p-k+ie 


-00 


fl 

du[27t(p2+4kpu)  (p2-4k2)+np2  (l-u2)(15p2+30'*pu-12k2)] 

-1 


Jifl. 


(l+i2(p2+2kpu+k2)]2 


2k 


2  f 

pi  -ine 


Carrying  out  the  p  integration  and  letting,  x  ■  p£,  Iq  •  kH,  ao  ■  ka, 
we  obtain. 


02  .  ii-—  f  — —  m2i(ao,fo)  +  — • tn??(ao ,lo)? 
9ku  C  4  2ki 5  J 


(K» 


where 


ft21 (a0  *ao) 


10u3+u2-2u+l 


J.j  [1+2  e2  (l+«)) 2  4 


du  - - exp 


(81) 


®22 (a0 


du  P 


1-1 


dx  x  ^ao  {x2+4Hqux+32.§  ]  ^ 
x-^o  8^0 


(82) 


x2  (1-u2)  (15x2+30: qxu-1 2/,2)+2(x2-4g.S)  (x2+41qxu) 
(l+x2+2f-oxu+lo)2 


For  the  SRI  experiment  ao  *  7.8  and  Iq  “  4.4.  Numerical  integration  yields 
m2i(7.8,4.4)  -  -0.12,  n22(7. 8,4.4)  -  203,  and  m!(4.4)  -  0.0475. 


The  fourth  term  in  (65),  03,  is  now  evaluated.  Treating  A(p-q)  as  a 
quasi  delta  function  and  writing 

_ JL _ _  1  f _ 1 _ _J. _ 

(p2-k2+ie) (q2-k2-ie)  (q-p) • (q+p)-ic  (  p2-k2+ie  q2-k2-ic 


we  obtain  approximately 


03 


2c2v~1 (2k) 
9k“  Ao 


h3(p) :(p+k) V(n-k)( - - - 

p2-k24*ie 


l 


2-ic 


)  ir  ;‘  2,(q  -p) 

'  .k  .A 

; 2n. (n-n)-ie 


(33) 


Carrying  out  the  03  integration  and  recognizing  that  the  principal  value 
integral  vanishes,  A(q)  ■  A(-q),  v;e  find  the  q  integration  reduces  to 


dq  A2(p-q) 

Jk  A  -A 

2p* (q-p)-ic 


dqidq2A2(qi ,q2» 


_  Pi3l."P2S2.) 

P3 


(84) 


But 


1 _ 1 

p2-k2+ic  p2-k2-ic 


[5(p-k)+5(p+k)] 


(85) 


and  model  (5)  for  A(p)  in  (84)  yields 


dq  A2(q )  m  iit2</ffA$ 
2p*q-ic  ab 


(86) 
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V/hen  this  result  is  used  into  (75)  and  we  introduce  spherical  coordinates 
we  have  after  $  integrations 


o2 


n  ,  1 
- Heal  v 

9k4A0  v 


dp  — - 

p-k+ie 


f  du(2Tr(p*’+4kpu)  (p2-4k2)+:ip2  (l-u? )  (15p2+30kpu-12k2)  ] 


_  _ 

*  U+T2(p~2f2kpu+k2)]2 


Trfa^ai 

2k 


2 

+  2^D(-a-^) 


l 


Carrying  out  the  p  integration  and  letting  x  ■  pi,  £q  *  ki,  ao  *  ka, 
we  obtain. 


o2 


ic  r 

9kk  C  i 


m21 (ao  .?o) 


a 


2kS 


,  ?q) 


i 


(80) 


where 


®2 1 (a0  »lp) 


.1 

.-1 


10u3+u2-2u+l  ,  a§(l+u)2v 

- -  exp  (_  - ) 

[1+2^q  (1-Hi)  ]2  4 


n,*  .  *  •  * 


(81) 


®22 (a0 .f-p) 


-f 


du  P 


k 


dx  x 
x-Iq 


»( 


ap  [x2+4S,pux+3ft§  ] . 


m 


(82) 


xz  ( 1  -  u ' )  ( 13x2+30i pxu-1 2.' j’)-f-2(x2-4?J5 )  (x2+4S.qxu) 
(l+x2+2i0xu+t§)2 


\ 


For  the  SRI  experiment  ag  ■  7.8  and  &q  ■  4.4.  Numerical  integration  yields 
m21  (7. 8,4. 4)  -  -0.12,  m22(7.8,4.4)  -  203,  and  1^(4. 4)  -  0.0475. 


The  fourth  term  in  (65),  03,  is  now  evaluated.  Treating  A(p-q)  as  a 
quasi  delta  function  and  writing 

(p2-k2+ie)  (q2-k2-ie )  (q~p)**(q+p)*-ie  ^  p2-k2+iE  q2-k2-ie 


we  obtain  approximately 


03  .  2e24>~1  (2k) 
9k4  A0 


dp  h3(p)  !  (p+k)  v(n-k)  ( - -1 - - 1 - ) 

p2-k2+ie  p2-k2-ie 


dr  .\2(q-o) 


(33) 


2p. (n-n)-ix 


Carrying  out  the  03  integration  and  recognizing  that  the  principal  value 
integral  vanishes,  A(q)  =  A(-q) ,  v;e  find  the  q  integration  reduces  to 


dq  A2 (p-q ) 

jk  j.  Jk 

2p* (q-p)-ie 


ijr 

2P3 


t 


dqidq2A2  (qi  ,q2  »  “  -E-^-^-2--2-) 

P3 


(8-0 


But 


1 _ 1 _ 

p2-k2+ie  p2-k2-ie 


~  [s<p-k)+<5(p+k)] 


(85) 


and  model  (5)  for  A(p)  in  (84)  yields 


dq  A2(q)  m  1tt2»|/ttA§ 
2p.q-ie  ab  pp 


(86) 
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When  (35)  and  (8b)  are  used  (S3)  we  obtain  v;ith  model  (5)  for  <i>(p)  and 


Aq  -  ira2b,  >j>q  ■  8tt fi. 3 ,  e 


k2i:2. 


- ao  *  ka,  *o  a  6  *  (~ ~ — )‘ 

(2tr)3  1+2  £q 


(87) 


03 


fi4Pl  H  ao*o 
1  8 


16J(B)(l+4fcfo 

*(H2l5)4 


where 


(88) 


J(8) 


r  tt/2 

dO  sin40 


v/2 

dt  (l-3sin2Scos2<>)”2 


'  0  1  0 

Making  use  of  the  substitution  x  =  l-3sin20cos2 if>  we  can  carr>  out  the 
double  integration  in  terms  of  complete  elliptic  integrals: 


(89) 


J(3)  - - -n- -  ((2S2-3 3+4)i:(3)  +  (B-4)(1-3)K(0)]  (90) 

123 2 (1-3) 

Recognizing  that 


e(B)  ■  §  (i  -  |  -  5  (yf£>2e2  - 


6+0  2 


K(3)  -  \  (1  +  (y)2P  +  (|~)202  + 


S-K)  2 


lim  r/a*  .  lim  „ ,  /  16x^2 

g+l  E(B)  ■  1*  ^(3)  ■  In  (j__g ) 
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we  obtain  the  following  limiting  values: 


J(0) 


3  IT2  .  J(3)  TT 

32  m‘a  3-*-l  =  4(1-3) 


A  curve  of  J(6)  is  shown  in  Fig.  1.  We  note  that 


lim 

3->l 


16j(8)  (1+41.  §)  2 
ir(l+2*§)4 


4(1+4 
(1+2J ll)2 


Next  we  use  results  (73),  (80)  and  (83)  in  equation  (65)  for  the  direct 
polarization  cross  section.  With  the  definitions  in  (37)  and  writing  the 
first  Born  radar  cross  section 

4wkV  A00(2k)  4m!4 p  a§b0i3 

o.  - - PI - U -  (93.) 

(4tt)  2  2k2  (1+4  J,2) 2 

We  obtain 


0(-k) 


ob!qI2[i  +  ^4Pi 


~  ag  ^  o  l"*-  + 


m2?  +  16J(8)  ( 1+4 il q ) 2  j 

9ir3/2^  it  (1+2S.2)4  J 


(92) 


According  to  (64),  (62)  and  (46) 


1  -  i  - 


~  2~C2T)  j 
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■)  ir2k$og(£o) 


(93) 
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where 


gUo)  ■  ~  J(l+2?.J)ln(l+i2£0)  -  --2— (£0  [l+4£§+2£j]+i[l+5i J+54 $])( 
2^o  1 J 


l+i02+i6  5,Q3 

+  - - - - - 

3£3(l+4£2) 


or  separating  into  real  and  imaginary  components 
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When  definitions  (87)  are  used  in  (93)  we  obtain 
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1  J  4/2 


For  ao  *  7.85,  bo  ■  15.7,  Iq  ■  4.4,  2.Qg(£  )  *  1.05  +  i  0.282, 
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38^8 

1+2Z2q 
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39.8 

0.95,  J(0.95)  =  15.5 


(98) 


Applying  the  numerical  values  to  (92)  we  obtain  for  the  direct 
polarization  radar  cross  section 


ct(rc) 
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_ Pj _ 
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(99) 


4 .  Cross  polarization  backscatter  cross  section 

Next  we  specialize  the  general  cross  section  formula  (Report  2,  eq.  28) 
to  determine  the  cross  polarization  backscatter  cross  section.  Setting 

_  _  _  k2Sl2p  .  ^  * 

ic  ■  -kt  letting  d  c  * - L  j  and  recalling  that  Q  -  ZQ  we  observe 

4tt 


-25- 


that  the  first  two  terms  in  (P.2.28)  vanish  and  we  have  for  the  cross 
polarization  cross  section  (o^) 


Oj.(-k) 


2c2(c*  G*(p)  -Q)  (c-  C  *(q),Q*K(p-H'  )*(p-k)A2(q-p) 
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dp  hj;  (p)  v(p-k)  >  (p+k) 


do  A2 (q-p ) 
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Jk  *  *  -S 

where  h1(p)  »  (j*p)(z*p) 


For  model  (5)  the  integral  terns  in  (100)  reduce  to 
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(101) 


With  p  expressed  in  spherical  coordinates 

h£(pk)  *  k1*  sin20cos20sin2p  (102 

When  turbulence  model  (4)  and  (102)  are  used  In  (101)  and  the  result  is 
used  in  (100)  there  results 
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where 


Ja(6) 
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(104) 


Jx(P)  Is  shown  in  Fig.  1  for  0.7  <  fi  <  1.  For  the  numerical  values  used 
in  the  previous  section  the  cross  polarized  radar  cross  section  is,  with 


J  (0.95)  *  0.235, 


Vc) 

0a(rc) 


(1+0. 34n40  )2+7.m8_  Pl  /? 
1 1  Pi 
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Plots  of  equations  (99)  and  (105)  are  shown  in  Fig.  2.  These  results 
compare  favorably  with  SRI  experimental  cross  sections  when  the  measure¬ 
ment's  uncertainty  is  approximately  3  db. 
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PREFACE 


In  recent  years,  an  increasing  amount  of  interest 
has  been  focused  upon  the  problem  of  wave  propagation  and 
scattering  in  random  media,  both  from  the  viewpoint  of 
theory  and!  measurement.  The  underlying  reasons  for  this 
interest  are  the  numerous  instances  of  wa^ve  propagation 
through  turbulent  regions  such:  as:  light  or  short  wave 
communication  in  the  atmosphere,  under  water  scattering,  and 
scattering  from  reentry  objects. 

To  investigate  these  problems  analytically,  one  has  to 
utilize  a  mathematical  model.  Unfortunately,  the  model 
which  presents  the  best  physical  description  is  mathematically 
intractable.  This  is  due  to  the  complicated  dependence 
of  the  medium  upon  its  physical  parameters ,  :.i.e.  temperature, 
pressure,  humidity,  electron  density,  velocity,  etc.. 

In  Chapter  I,  the  choice  of  a  simpler  model  which  may 
describe  the  prominent  characteristics  of  the  wave  medium 
interaction  is  adopted.  'In  Section  B.,  the  medium  is 
mathematically  defined;  the  mean,  rms,  and  correlation 
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functions  are  presented  as  adequate  measures  of  the 
randomness  in  the  medium.  In  Section  C. ,  the  problem 
explored  in  this  thesis  is  mathematically  formulated  for  a 
layer  of  infinite  extent  in  the  yz  direction  and  of  thickness 
x#  in  the  x  direction.  This  layer  is  illuminated  by  a  plane 
wave  propagating  normally  to  the  fluctuations  in  the  index 
of  refraction.  This  layer  is  also  characterized  by  a 
collision  frequency,  a  homogeneous  mean  and  a  strong 
small-scale  fluctuation  in  the  electron  density.  The 
resulting  random  differential  equation  of  the  field  is 
formalized  and  its  difficulties  are  explained.  Since  a 
straight  Neumann  series  solution  fails  under  the  present 
circumstances,  selective  summation  techniques  are  used. 

Chapter  II.  develops  these  techniques  (diagram  and 
smoothing)  with  the  appropriate  approximation  to  account 
for  large  fluctuations  and  multiple  scattering.  Though 
both  techniques  lead  to  the  same  integro-dif ferential 
equation,  the  derivations  are  completely  different; 
consequently,  these  derivations  shed  light  on  different 
aspects  of  the  problem  and  on  the  interpretation  of  the 
results. 

The  preceding  techniques  have  been  used  in  the  study 
of  wave  propagation  in  infinite  media.  In  1962,  Bourret 
*  (5),  (6),  derived  the  first  order  Dyson  equation  ( 1 5 ) 

through  the  hypothesis  of  local  independence;  he  applied 
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the  resulting  equation  to  wave  propagation  in  an  infinite 
homogeneous  random  medium.  In  a  series  of  papers,  Keller  re¬ 
derived  (2>6)  ,  (28) ,  the  same  results  and  avoided  the  above 
assumption.  In  1963  Tatarski  and  Gercenstein  (41)  derived 
Bourret's  results  using  a  variant  of  the  smoothing  method, 
tfhe  resulting  integro-differential  equation  was  applied 
to  scalar  wave  propagation  in  an  infinite  homogeneous  medium 
with  strong  small-scale  fluctuations.  In  1964  Tatarski  (4$) 
studied  the  electromagnetic  wave  propagation  in  an  infinite 
medium  with  strong  dielectric-  constant  fluctuation. 

Macrakis  (32.)  .  (1965)  gave  a  more  straightforward  derivation 
of  some  of  Tatarski's  results.  Brown  (8)  (1967)  applied 
Tatarski' s  formulation  to  the  propagation  problem;  he 
interpreted  the  results  in  terms  of  coherent  and  incoherent 
scattering,  and  the  theory  of  dielectric.  Frisch  (16)  and 
later  Bassinini  (2)  (1967)  obtained  the  known  solution  of 
the  classical  random  oscillator  through  the  smoothing 
method. 

The  application  of  the  Dyson  or  smoothing  formulations 
to  a  bounded  or  inhomogeneous  random  medium  has  received 
little  attention.  Bassinini  et  al.  (  I  )  (1967)  used  the 
infinite  medium  effective  wave  number  parameter  given  by 
Tatarski  to  study  scattering  from  bounded  media.  Recent 
•papers  (Rosenbaum,  1969;  Kupiec  et  al.,  1969;  Collin,  1970) 
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studied  some  aspects  of  interface  effect  on  coherent 
wave  motion.  Rosenbaum  (4/)  obtained  approximately  the 
coherent  wave  scattered  by  a  random  half  space  whose 
mean  dielectric  constant  is  different  from  the  uniform 
nonrandom  medium  in  the  conjugate  half  space.  For  a 
homogeneous  background  throughout,  Kupiec  et  al.  (30) 
solved  the  scalar  Dyson  equation  for  a  normally  incident  wave 
on  a  random  medium  occupying  a  half  space  and  slab  geometries. 
Collin  considered  an  obliquely  incident  wave  on  a  random 
half  space  superimposed  on  a  uniform  medium  throughout. 

It  is  difficult  to  assess  the  results  obtained  since  the 
randomness  was  assumed  to  be  homogeneous  near  the  boundary 
A  complete  investigation  of  the  bounded  random  or 
inhomogeneous  random  problems  is  difficult.  In  this  thesis 
some  aspects  of  these  problems  are  considered  and  applied 
to  a  randomly  fluctuating  plasma  layer  whose  uniform 
background  is  generfsilly  different  from  the  outside  medium. 

In  Chapter  III,  the  conjecture  made  by  Bassinini  et  al. 
about  the  effective  wave  number  is  tested  through  ap¬ 
plication  to  the  layer  problem.  The  reflected  coherent 
and  incoherent  powers  are  derived.  The  comparison  with 
the  successive  approximation  method  is  also  made.  In 
Chapter  IV,  the  exact  dependence  of  the  scatter  power  on 
large  electron  density  fluctuations  is  investigated  in  the 


fp 
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limit  of  small-scale  fluctuations.  A  locally  homogeneous 
randomness,  which  appropriately  describes  the  medium 
near  the  boundary,  is  used,  Typical  cases  are  studied  for 
a  wide  range  of  plasma  parameters.  In  Chapter  V,  a  new 
formulation  of  the  scattering  problem  is  presented  for  a 
Utiiform  i  background  throughout  with  a  locally  homo¬ 
geneous  randomness,  whose  mean  characteristics  change 
smoothly.  This  property  is  advantageous  in  the  Fourier 
transform  space,  and  a  singular  integral  equation  of 
standard  form  results.  This  description  appears  to  be 
closer  to  reality  since  it  avoids  the  presence  of  a 
sharp  boundary  between  the  stochastic  and  nonstochastic 


medium. 


I.  FORMULATION  OF  THE  PROBLEM 


A.  Introduction 

In  this  chapter,  the  medium  and  its  interaction  with 
an  incident  wave  is  mathematically  described.  This 
description  leads  to  linear  differential  equations  whose 
coefficients  are  dependent  upon  the  macroscopic  inhomogenei¬ 
ties  of  the  medium. 


B.  Description  of  the  Medium 

Two  types  of  inhomogeneities  exist,  deterministic  and 


random.  It  is  deterministic  if  an  identical  measurement 


is  performed  many  times  and  the  results  obtained  are  always 
alike.  If,  however,  all  conditions  under  the  control  of 
the  experimenter  remained  the  same,  the  results  continually 
differ  from  each  other,  the  medium  is  said  to  be  random. 

For  a  "temperate"  and  isotropic  plasma  with  electron 
density  fluctuations,  the  complex  dielectric  constant  is(iq) 


where  JT-  ,  are  respectively  the  normalized  plasma 

and  collision  frequencies;  the  quantity  is  the  electron 
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density  distribution  normalized  on  the  peak  or  mean  value. 

i 

Here,  the  mean  value  of  is  unity  and  the  fluctuation^ 

is  random. 

In  the  applications  of  interest  to  the  present  in- 
% 

vestigation,  the  fluctuations  in  electron  density  are 
aSsum&cfc]  to  be  related  to  turbulent  mixing  (48)  ,  (50)  . 

The  mean  and  second  moments  of  the  fluctuations  usually 
suffice  for  most  scattering  calculations.  The  mean  in¬ 
dicates  a  steady  value  of  the  electron  density  over  the 
period  of  observation.  The  second  moment  measures  the 
correlation  between  electron  density  fluctuations  at 
neighboring  points.  Ensemble  averaging  over  a  set  of 
realizations  is  assumed  throughout; 

For  homogeneous  turbulence  the  second  moments  are  given 
by 

=  a£U,)  B ^2) 


For  locally  homogeneous  turbulence,  we  have 


(3) 
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The  exponential  and  Gaussian  correlative  functions  for 
isotropic  turbulence, 


B(  v  O 


.  I*i« 

e  — 


(4) 


B(x,- 


(5) 


ate  useful  in  many  applications. 

C.  Derivation  of  the  V7ave  Equation 

Consider  a  plane  wave  of  the  form  s>eJ 
normally  incident  on ‘a  plasma  layer  (fig.  1).  The  slab 
is  characterized  by  a  complex  dielectric  constant  which 
is  a  function  of  space.  The  medium  is  assumed  to  be 
linear  and  isotropic.  The  electron  density  possesses  a 
homogeneous  mean  and  a  strong  small-scale’  fluctuation. 

Maxwell's  equations  for  all  space  can  be  written  as 


(6) 


:  -  jwf**  H 

V  x  H  &  Jwi#  K(x)  E 


(7) 


Fig.  1  Random,  fluctuations  of  electron  density. 
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where  K(x)  is  the  "effective  dielectric  constant",  equal 
to  one  outside  the  scatterer  and  varying  with  x  inside  the 
scatterer. 

The  wave  equation  for  E  is  found  by  vector  manipulation 
of  (6)  and  (7).  Let  =  “'V* fc#  ,  then, 

KU)  1  ' 

In  cases  where  the  electric  field  vector  is  normal  to  the 
variation  in  the  refractive  index  the  cross-coupling  term 
in  (8)  vanishes  ;  let  be  a  component  of  the  electric 
field  vector,  then 


♦  k'kco)  <fr  *o 


(9) 


Complete  mathematical  solution  of  this  differential 
equation  with  random  coefficients  is  still  lacking.  Here, 
the  study  is  confined  to  the  investigation  of  the  statistical 
properties  of  the  solution,  $(*)  ,  of  the  random  equation. 

To  facilitate  the  analysis,  equation  (9)  maybe  put  in 
operational  form: 


11 


L  <f>  -  o 


(10) 


where 

\ 


L  sL+L 


£ 

<Ufc 


l,  . 


k'  ~  VK.  *  C»  -  N1 


fit  -  ( *  +  J  ^ 

o  ♦-*;) 


(id 


together  with  the  continuity  conditions  at  the  boundaries, 
and  radiation  conditions  at 

Equation  10  is  the  reduced  wave  equation  which 
occurs  in  classical  and  quantum  fields  ($6) ,  (17). 

The  solution  of  the1 reduced  wave  equation  must  be 
found  approximately  because  calculation  of  L  involves  a 
.  nonuniform  component  in  L#  and  the  random  component  L,  . 


.12 


Fortunately,  physical  interest  is  focused  not  on  the 
functional  dependence  of  ,  but  on  the  statistical 

moments  of  the  wave  solution.  The  moments  are  in  fact  the 
observable  features  of  the  physical  phenomena;  it  will  be 
our  purpose  to  investigate  the  coherent  field  <$(*)>  and 
the  fluctuating  field,  .  Specifically,  this  thesis  is 
concerned  with  the  dependence  of  the  scattered  field  on  the 
large  fluctuations  of  the  refractive  index. 

The  procedure  employed  consists  in  the  calculation  of 
the  inverse  operator,  L"' .  This  formulation  will  be  the 
subject  of  the  next  chapter,  and  reduces  to  finding  the 
effective  characteristics  of  the  density  fluctuations,  and 
subsequent  determination  of  the  coherent  and  incoherent 
scattered  power. 
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II.  THE  METHODS  OF  PARTIAL  SUMMATION 

A.  Introduction 

The  validity  of  the  widely  used  successive  approximation 

method  is  restricted  to  bounded  media  and  small  strength  " 

random’  perturbations  (Appendix  I)  .  The  problem  under 

investigation  does  not  satisfy  these  conditions;  and  in 
• 

addition  the  summation  of  all  terms  up  to  second  or  third 
order  does  not  wprk  because  the  series  c  onverges  too 
slowly  to  be  practical.  In  these  cases,  two  approximation 
methods  (smoothing  and  diagram  methods)  are  available; 
they  will  be  described  in  this  chapter.  The  diagram 
method  (14  )  involves  a  summation  over  a  selected  class  of 
repeated  diagrams.  The  method  of  smoothing  (49)  involves 
the  decoupling  of  the  mean  and  fluctuating  field  through 
the  solution  of  the  latter  field  by  formal  iteration;  an 
equation  for  the  mean  field  is  then  obtained.  Though 
both  methods  yield  the  same  results  (l5)  when  applied 
to  linear  random  equations,  their  methods  of  derivation  shed 
light  on  different  aspects  of  the  problem.  In  Section  Bi 
the  development  leading  to  the  Dyson  equation  (l5)  is 
carried  out.  In  Section  C,  the  smoothing  equation  is 
derived 

X.  The  Diagram  Method 

Let  the  field  $(«0  in  a  random  medium  with  index  of 
refraction,  ,  be  related  to  its  source  J  through 
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<*>(*) 


]<*:)  d*' 


(1) 


The  Green's  function,  ,  satisfies  the  random 

equation  in  integral  form  s 


G  (*,%')  =  6  (*,*')  -  /  G  L,  l**,/-")  GCx,*')  eU, 


(2) 


when  GU,  *'}  is  the  free  space  Green's  function  (57), 


-  i-  ejkU*Vl 

#  -.Tin 


(3) 


The  solution  of  (5  )  by  iteration  gives 


G(x,k') s G t*,>0 •  U, Wcti  +  *»*  t  (4j 

♦  /g^U.hJL, (*)... 


An  obvious  physical  interpretation  of  the  n**1  term  is: 
a  wave  propagates  unhindered  from  n*  to  1 ,  where  it  is 
scattered  by  the  random  inhomogeneity,  the  resultant 
propagates  unhindered  from  1  to  2,  is  scattered  at  2,  etc.. 

If  the  index  of  refraction  is  normally  distributed, 
the  average  of  the  products  of  the  random  functions  in 
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the  terra  is  (8  ) 


<1,(0  LtC>%)  ^  s  2.  <L,<»)L#U>  )•..  h  te/’mj  <  L,  W  } 


or 


<lv(i)Lk(»)-lkN)> 


i ••• 


e- 


*»ftv«n  (6) 

*.dd  (7) 


and 


r:o.x)  =<  m»>  m&)  > 


(8) 


^(».i),the  correlation  of  two  points,  tends  to  zero 
for  separations  large  compared  to  the  correlation  length. 
The  stun  in  (  5  )  is  constructed  from  possible 

permutation  of  the  indices  (8)*  This  summation  is 
overfall  the  combinations .of  the  points  1,  2,  ...n  into 
groups  of  two 


IJtl,!),  ...  f  tV***!,*,) 

To  facilitate  the  writing  of  £be<.. functions 
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<G(»,  *0  >  and 


\ 

equivalent  diagrams  of  each  term  in  the  mathematical  series 
is  drawn.  A  dictionary  of  the  fundamental  symbols 

^  ) 

may  be  constructed  as  follows:  We  represent  each  G#f 
by  a  horizontal  line  segment;  if  we  have  two  G*** 
we  join  the  two  lines  and  put  a  dot  at  the  junction.  The 
dot  indicates  the  coefficient  of  the  term  times  the 
integration  over  the  random  inhomogeneity;  A  two  point 
correlation  is  identified  by  a  dashed  line  connecting  the 
two  random  inhomogeneities.  The  multiplication  of  two 
diagrams  is  indicated  by  placing  one  on  top  of  the  other. 

The  Neumann  series  of  the  Green's  function  Gc*,*') 
may  be  diagrammed  as  follows: 


Its  product  with  itself  will  be 


* 

f 


Since  only  L,  indicates  randomness,  the  means  of (9)and  (10) 
are: 


♦ 


+ 


•***'*  •**'v 

f  T  ■  6  ■  ■  fc . i  ■  »ip 


(ID 


where 


<&(«,«')  >  -  . 


<G  *  G*  >  u 


It  is  proven  (Appendix  A.)  that  the  Green's  function  series 
is  convergent  for  a  sufficiently  small  perturbation: 
the  terms  in  the  expansion  are  then  successively  smaller 


with  increasing  order.  In  this  case  we  can  estimate  the 
accuracy  of  the  approximation  from  the  highest  order 
term.  For  large  perturbations*  the  above  process  fails. 
In  this  case,  the  terms  that  appear  to  give  the  largest 
contribution  are  chosen  and  summed  to  an  equivalent 
deterministic  integral  equation.  Physical  arguments  at 
times  enable  us  to  identify  the  terms  to  be  chosen. 
Though  the  diagram  technique  is  lacking,  rigorous 
mathematical  justification,  it  is  useful  since  it  yields 
results  which  compare  favorably  in  some  cases,  with 
exact  methods  (  2-  )  . 

To  perform  the  summation,  each  term  in  the  diagram 
equation  of(llJis  dismantled  into  noncorrelated  parts. 

For  instance 


Factorization  of  noncorrelated  parts  in  the  diagrams  lead 

to  the  Dyson  (12.)  and  Bethe-Salpeter  (4*)  equations  of 
•  • 
field  theoretic  methods: 


In  operational  form,  we  have 


where 


where 


+ 


+ 


% 


The 'effective  wave  number  operator", |  |,  and  "intensity 

,  are  the  sum  of  the  diagrams  that  cannot 


operator". 


be  further  dismantled  without  breaking  dotted  lines. 

To  construct  the  first  approximation  to  the  Dyson 
equation,  the  property 


t|l  -  xi )  •  s  o 


U-*t  >  l 


(16) 


reve&ls  that  the  correlation  between  neighboring  points 


the 


(18) 

It  is  to  be  noted  that  the  solution  of  this  equation 
should  not  be  performed  by  iteration  as  this  process 
leads  bacfc  to  a  divergent  series.  This  is  expected  as 
the  Dyson  equation  is  a  selective  restructuring  of  the 
Neumann  series. 

C.  The  Method  of  Smoothing 

The  smoothing  method  (l5)  to  be  discussed  below  is 
an  alternative  development  to  the  diagram  method.  For 
random  equations,  the  two  methods  lead  to  the  same  results 
Some  ad  hoc  procedures  are  also  discussed  in  Appendix  B*. 

To  facilitate  the  an;  ysis,  the  operational  form  of 
the  wave  equation  is  adopted: 


where 


4>  =  P$  + 

m 

(20) 

F  %  <  > 

£*  5  Ci-  p)  * 

By  applying  successively  the  projection,  P  ,  and 
fluctuation,  $$  ,  operators  to  equation  (19)  we  obtain 

<  ♦  >  -  J  -  G^  P>^  5  ♦ 

(21) 

i  $  s  .  G^  Cl-  P)t,  (<*>  +  J*) 


Solution  of (20)  by  formal  iteration  yields 
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If  we  replace  Sty  in (20)  by  its  value  in (22)  ,the  coherent 
field  becomes: 


<*>  S  Gy  J  ♦  M  <  $> 


with  the  mass  operator 


(23) 


M  .  £  PSC- C.  <t -OL,)  P 

«UI  ' 

,.£p«..  <-Cl.  >"p 


(24) 


Equation  (23)  represents  the  Dyson  equation  for  the 
coherent  field.  Above  order  one,the  Dyson  equation  is 
usually  intractable  because  of  the  increasing  complexity 
of  the  mass  operator.  In  what  follows,  the  analysis  and 
discussion  are  concerned  only  with  that  first  approximation; 
for  n*l,  equations  (21)  and  (23)  reduce  to 


V,  <  ♦  }  *  <  ><  ♦  >  *  ° 


(25) 


(26) 


The  solution  of  Equation  (25)  has  been  shown  (41)  to  be 
asymptotic  to  the  exact  solution  of  the  classical  random 
oscillator  for  the  condition  of  strong  small-scale  fluctu¬ 
ations  . 


* 

D.  Some  Observations  on  the  Partial  Summation  Method 
,  The  methods  of  partial  summations  transformed  the 
stochastic  integral  equation  (12)  into  a  nonstochastic 
integro-differential  equation,  (25),  which  possesses  the 
mathematical  form 


-  Ik<*,**>  <  >  d*'  so 


(27) 


Where  J  is  the  kernel  and  <  *  >  ,  the  unknown 

coherent  field. 

A  general  solution  of  Equation  (27)  is  not  known.  If 
a  straightforward  iteration  procedure  is  used,  this  leads 
to  Equation  ( 18) ,  which  is  known  to  diverge  in  the  case 
of  large  fluctuations.  When  the  kernel  is  of  the  convolu¬ 
tion  type (4  4) 


the  application  of  Fourier  transform  leads  to  .a  solution 
of  the  wave  equation.  Also  if  the  integral  is  of  the 
Wiener-Hopf  type {44 ) 
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d* c*  >  +  <  P  )  •  f  *')  <  &(*')  }  ApS  i  o 

J*'  ‘  #  (27) 

.«  <•<« 

% 

a  technique  is  known  for  its  solution.  Both  solutions  have 
• 

been  discussed  in  the  literature.  Physically,  both 
cases  required,  among  other  things,  the  assumption  of 
random  homogeneity;  this  seems  hardly  justifiable  near 
the  boundaries  of  bounded  scatterers.  In  addition,  the', 
appropriate  Green's  function  that  satisfies  the  radiation 
condition  at  «»  and  the  boundary  conditions  of  the 

scatterer,  is  not  of  the  convolution  type.  These  physical 
modifications  add  to  the  complexity  of  the  scattering 
problem  considered  in  this  thesis. 

E.  Conclusion 

Further  work  on  the  mathematical  foundations  of  the 
Dyson  and  smoothing  equations  are  needed  (>3) .  At 
present,  the  capacity  and  limitations  of  the  first  order 
Dyson  or  smoothing  equations  are  best  gauged  by  the 
examination  of  the  results  of  their  application  to  / 

prototype  problems.  The  resulting  integro-differential 
equation  is  difficult  to  investigate  in  the  case  of  random# 
inhomogeneous  or  bounded,  media . 
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III,  WAVE  SCATTERING  FROM  A  RANDOM  UNIDIMENSIONAL  LAYER  ' 
WITH  AN  APPROXIMATE  EFFECTIVE  INDEX  OF  REFRACTION 

A.  Introduction 

Recent  papers  (S  ),  (32),  (21),  considered  the 

\r 

propagation  in  an  infinite  random  medium  with  strong 

small-scale  fluctuations.  For  this  case,  the  resulting 
• 

integro-dif ferential  equation  of  the  first  smoothing 
equation  is  of  the  convolution  type;  its  solution  is 
readily  obtained  through  Fourier  transformation. 

However,  examination  of  a  bounded  or  an  inhomogeneous 
medium  requires  extreme  simplification  of  the  stochastic 
problem.  In  a  recent  paper,  Bassanini  et  al.  (  t  ) 
asserted  that  the  interface  effects  on  the  coherent 
effective  wave  number  parameter  are  negligible  in  the 
limit  of  small-scale  fluctuations;  they  concluded  that 
the  development  made  by  Tatarski  (41)  in  the  evaluation 
of  an  effective  parameter  for  an  infinite  random  medium 
with  uniform  background  is  applicable  to  bounded  media. 

,  Since  the  mathematical  foundation  of  this  procedure  is  not 
wall  set,  its  results  are  subject  to  doubt.  This  chapter 
may  be  viewed  as  a  numerical  testing  experiment  to  see 
if  the  above  assertion  has  not  produced  a  violation  of  the 
principle  of  conservation  of  electromagnetic  energy.  The 

■  reflected  coherent  and  incoherent  powers  are  analytically 

* 

obtained  as  a  function  of  the  strength  of  the  fluctuation. 
Comparison  to  the  successive  iteration  method  is  also  made. 
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B.  Infinite  Medium  Effective  Index  Parameter 

In  this  section  we  reduce  the  original  stochastic 
problem  in  its  coherent  part  to  an  equivalent  deterministic 
boundary  value  problem  through  an  estimate  of  the  effective 
characteristics  of  an  infinite  medium. 

The  Green's  function  is  taken  as 


t-Jfc  |  »  -  »  > 


(1) 


The  randomness  is  assumed  to  be  homogeneous  and  the 
correlation  is  exponential.  The  field  within  the  plasma 
is  governed  by  an  integro-differential  equation  of  the 
convolution  type;  its  solution  has  been  carried  out  through 
the  application  of  Fourier  transformation  (41)  or  through 
the  assumption  of  a  wave  solution  (1$)  of  the  form 


<  *«> 


•i  V.  .«• 

>  -v  ■ 


Following  the  second  alternative,  one  obtains,  after  some 
manipulations 


i 


*  k  (  *C4  .  j  fa  ) 


(2) 


where 


*  [  la  ♦  0*) 


y')/ifv 


+(A‘*  U%)VtJ/t 


*,  ♦  (r>,  ♦ 


r* 


It  -Ap  c*ft.  Cl  -  -A*.  >/(•  ♦  -n£  ^ 


«  .  /  ,.t  wi 

g  s  akA^f  ft.-*1*  a  i  ♦  a. ) 

«,*  r<  *,♦<** ♦  <>Vi>/i  jVt 

.  jt,  s  U-*.  ♦  K4+0  \£  ^ 

*•»*'•  AtAt*  ^ 

•*5  Ar Jl*  ^ 

For  Vs®  ,  Equation  ( 2  )  reduces  to  the  deterministic 
plasma  wave  number;  for  c  t  o  »  the  turbulence  attenu¬ 
ates  the  wave;  it  transforms  the  coherent  energy  of  the 
incident  wave  into  incoherent  energy;  this  effect  is  more 
prominent  a#*s  the  normalized  plasma  frequency  -ftp 
the  over-dense  regime. 


nears 
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C.  The  Ooherent  and  Incoherent  Reflected  Powers  for  a_  slab 
1»  Ooherent  Raflected  Field 
The  analysis  is  confined  to  the  simple  bounded  slab 
filled^ with  a  randomly  fluctuating  plasma  and  described 
in  Chapter  I  (fig.l  ).  Assuming  that  Bassanini's  work  is 
correct,  the  coherent  field  equation,  within  the  slab, 
may  be  written  as 


i,4 ^ >  4  !»*<♦>  sp 


The  general  solution  of  the  coherent  field  equation  is 


>r  At 


n* 


♦  &  * 


(5) 


»  .V  T-  Wfc  A  .11  f  V.rv 


where 


r  r**»  • 


k  (  ♦  j  ) 


(6) 


The  boundary  conditions  are 


!♦<*;>  *  <  > 


i  -  - 


.  JL  &.<+(•)> 

jw  a* 


(7) 
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<4j>  c",l,<#r  <  4>C«.>  > 


(8) 


-jk«. 

- 


i.  4.  <  > 

)k  dy. 


\ 


Solution  of  the  preceding  algebraic  equations  will  deter 
nine  the  four  unknowns.  A,  B,  <4>m  >  *  <  %> 


A*  »  ** 


©  S  -  j  ir*  *• 

W,  .  .  w-<  *  •>)*.' 

<♦»>  s  i'  4 


(9) 


Where 

t  :  I-  WjW 

i,  v-  *■;*''*• 


(10) 
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2.  Incoherent  Reflected  Field 
The  incoherent  field  equation  may  be  written  as 


A*  -  1,4  «►  > 

v  a> 


(ii) 


The  solution  of  Equation  (u)  by  the  "variation  of  paramet¬ 
ers"  ctl )  Is 


S<f>  r  A  $*<»>♦ 


(12) 


♦  i  fhi.fr 

W(U„S^) 


I Lc+>  dx 
n(b,Mh) 


y  k  „r*. 

where  bfrt  r  e.  ,  $  s  e  ere  the 

homogeneous  solutions  of  (H)  ,  W  is  the  Wronskian. 

Apply  xr*^  Cw'.i  w xr/ui  w jf*  Cw/tditions  to  ?12)  #*  the  in**~~ 
coherent  field,  bp  #  within  the  plasma  becomes 


IfrC*)  =  Ci  c 


r* 


KK * \’r)  *  «*%t**JeT* 


.  1,  I?  eUr^ 


e^Vt^,,  y.-  y> 
A.C^  e*V,t»j  (u^y.y^  )  +  (*#v*fc*  *’*  *■*  *  *  } 


(13) 


where 


3.  Reflection  Coefficients 

The  reflected  power  is  now  found  as  a  coherent  and 
an  incoherent  contribution.  The  coherent  reflected  power 
ie  easily  calculated 


|<*>l\  |Lkl^£m2|' 


(14) 


The  incoherent  reflected  power  may  be  written  as 


(IS) 


For  the  exponential  correlation  function,  statistical 
averaging  reduces  (15  )  to 

=1 ttY'  -Mo-lN  | L  |* k  (- r -r,  -r’.  y ) 

1  YJt‘  ti*  -a‘.)  1  l 


- 1;  4’  6.;  t1*'**  K  <r,'.  y;-r.  -r>  -tVrt’ «■' 
kU'.r/.vr))  ♦ 


•».  v  yvll'“-k  (H  -r.-r.’-r*)  +£i>-.t.lV  ‘'V*'-*; 


Mt-V,  r,’.  r')J 4lt>/f‘^*r'J“-k|vr,r'.r.*> 


-i,  v  v>y1Y'*  k  i-v;  .y, ,.4>  t  r  y  *,•«**>;• 


“  ttf  -  Y*,r-  5,lt  14,  tJ  V  *c,**,'>  fr*-r.’,  r-t  > 
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V-<re) 

A‘  4'  i;  e"x<'f,*tr  4  ’i’**  kty*.r*,  >.  y.) 
A‘  A.  A.  t-  *  »’♦  *  *  4 »  *•  u  >•.  r*  ,  r  4  V.  ; 


|ij»  e.«T  ,  y,,J 


(16) 


where 


+(«., O »  --! - [4—*  —]*  [4- 

c,*ct  t-c,  £.et  C|*  *♦*. 


J — J  ■+  « 


C<i«  : 


+  «L 


(17) 


4+‘i  <*-«,>  C^ej  («♦«,)(*-  O 


For  a  nonstochastic  medium,  CT  s  o  ,  the.  incoherent 
reflected  power  zeroes  and  the  coherent  reflected  power 
reduces  to  the  well  known  deterministic  result. 

If  Y«re  the  reflected  coherent  and  incoherent 
powers  in  Equation  (14)  and  (16)  reduce  to  those  obtained 
by  Jarem  (**•)  using  the  successive  approximation  method; 
this  condition  is  approached  in  practice  for  very  small 
random  variations.  When  the  difference  between  the  mean 


values  of  the  slab  and  the  surrounding  medium  is  not 
too  small,  the  conjecture  made  by  Bassinini  et  al. 
becomes  unreliable  since  it  is  inconsistent  with  the 
continuity  conditions  at  the  boundaries.  In  fact,  the 
principle  of  conservation  of  electromagnetic  energy  is 
violated  (fig.  2). 

t 

D.  Conclusion 

The  basic  conclusion  drawn  from  the  analysis  of  the 
obtained  results  points  out  that  the  neglect  of  the 
interface  effect  on  the  coherent  field  is  not  generally 
close  to  reality:  1.  Some  observed  phenomena,  such  as 
the  saturation  of  the  scattered  power  before  the  critical 
density  is  reached  ,  are  not  displayed.  2.  The  conserva 
tion  of  electromagnetic  energy,  near  the  transition 
regime,  is  violated. 
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IV.  WAVE  SCATTERING  FROM  A  RANDOM  UNIDIMENSIONAL  LAYER 
WITH  AN  APPROPRIATE  EFFECTIVE  INDEX  OF  REFRACTION 

A.  Introduction 

In  the  previous  chapter,  use  of  the  infinite  medium 
Green's  function  was  shown  to  be  improper  for  a  slab, 
inhere  the  difference  between  the  mean  value  of  its  di¬ 
electric  constant  and  imbedding  medium  is  not  small.  In 
this  chapter,  exact  solution  of  the  problem  is  obtained 
in  the  limit  of  small-scale  fluctuations.  We  consider  a 
uniform  plasma  slab  with  superimposed  turbulence  whose 
intensity  and  short  scale  behavior  follow  a  Gaussian  form. 

Numerical  evaluations  of  typical  cases  are  displayed 
as  a  function  of  the  mean  square  fluctuations. 

B.  Exact  Effective  Index  of  Refraction 

The  coherent  field,  within  the  slab,  is  governed  by 
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The  Green's  function,  is  that  of  the  background 

medium  wave  equation 
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The  general  solution  of  Equation  (2  )  can  be  written  in 
terms  of  known  functions;  we  let  and  ^  be 

the  solutions  of  Equation  (2  )  which  satisfy  the  rad¬ 
iation  conditions  at  x-»  «  e©  and  the  boundary  conditions 
at  x«»i  *0  •  The  Green's  function  may  then  be 

constructed  (l41;(&5),  (44) : 
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W  is  the  Wronskian;  is  the  greater  of  the  two 

numbers  x  and  x',  and  5^.  is  the  lesser  of  x  and  x1.  Let 
us  determine  the  functions  and  for  the 

turbulence  free  uniform  slab  (fig. 2  ).  The  determination 
follows  from  the  standard  solution  of  the  slab  problem  , 
together  with  the  boundary  conditions  and  the  properties 
of  the  Green's  function;  this  leads  to  • 
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To  pursue  the  analysis  of  Equation  (  1 )  we  consider  the 
mean  value  of  the  field  that  has  a  scale  of  inhomogeneity 


much  larger  than  the  correlation  scale  of  the  perturbation. 
The  correlation  function  manifests  itself  then  as  a 
local  effect  (l>)#  (41) 


£(«,*)=  $U-  *f)  * 
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This  means  that  we  are  interested  in  the  variation  of  the 
coherent  field  on  a  scale  that  is  large  compared  to  the 
correlation  length.  In  this  case,  the  first  order 
smoothing  equation  becomes  for  al.locally  homogeneous 
fluctuation  whose  intensity  or  strength  follow  a 
Gaussian  form: 


c£<4>>  +  k*  Mt*>  e  o 


where 


H(t)  ?  Ci-iv)  -  GCic.x) 


This  equation  reduces  to  that  obtained  elsewhere  (Jo)  for 
the  case  of  an  infinite  uniform  background  and  homogeneous 
fluctuations  for  which  an  analytical  solution  can  be 
obtained  in  terms  of  known  functions.  However,  the 
solution  of  Equation  (8  )’  for  a  finite  slab  thickness  is 
not  known  analytically  and  must  be  integrated  numerically 
for  typical  cases.  No  difficulty  is  encountered  as  any 
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singularity  of  the  "effective  wave  number"  is  shifted  off 

the  real  axis  by  the  presence  of  the  randomness.  This 

shift,  which  appears  as  a  collisional  effect,  describes  the 

transformation  of  the  energy  in  the  propagating  wave  into 

fluctuation  energy.  The  transmitted  coherent  field  at  the 
• 

boundary  is  chosen  to  be  of  unit  magnitude  and  zero 
phase.  The  Runge-Kutta  method  (21)  is  used  to  carry  the 
integration  to  the  input  boundary,  and  the  reflection 
and  transmission  coefficient  are  then  automatically 
evaluated. 


C.  Scattered  Power 

1.  Coherent  Scattered  Power 

Consider  the  following  scattering  geometry  (fig.  2)  . 
The  change  in  the  positive  reference  direction  is  made  for 
eomontaHonel  convenience.  -  , 


The  boundary  conditions  are: 

x-x* 
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Let  the  complex  coherent  field  be  represented  in  terms  of 
real  and  imaginary  parts  as 


t  * 
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Let 
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Equation  (8  )  may  be  decomposed  into  a  coupled  set  of 
second  order  differential  equations  with  variable  coefficients 
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Let 

<4>%U)  >  =  ^.  an*  d<4>.C*)  >  =  ?c  (13) 

dx  At 


Equation  (12)  takes  the  form  (  ll  ) : 


Initial  values  may  be  obtained  if  we  take  the  wave  in 
the  slab  at  x=x#  as  the  reference;  for  convenience  let 
<♦<*>>=  I.  there  results  from  the  boundary  conditions. 
Equation  (9  ) 
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.  1J  (  x#)  so 


-  *$(».)  =  k 


We  solve  for  the  incident  reflected  and  transmitted 
coherent  field  as  a  function  of  the  field  within  the 
plasma  at  x=0;  we  obtain 
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The  coherent  reflection  and  transmission  coefficients  are 
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The  respective  phase  shifts  are 
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The  incoherent  reflected  ,  Rt  z  <  154^1*  >  ,  and 

transmitted,  s  <  !  S  4^  I  *  >  #  powers  are  given  formally 
by  (11,26)  .  in  the  one-dimensional  case  this  leads  to 
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Manipulation  of  (7  ),  (19)  and  (20)  give 
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D.  Discussion  of  the  Results 

The  first  Born  approximation  is  widely  used  in  the 
literature  (*fr$)  ;  this  approximation  assumes  that  the 
incident  wave  propagates  to  and  from  its  scattering  center 
unhindered  by  the  rest  of  the  slab.  Subject  to  this 
assumption  the  analytical  results  show  a  iinear  relationship 
between  the  scattered  incoherent  powers  and  the  mean 
square  fluctuations,  9-*  .  On  the  other  hand,  experimen¬ 

tal  observations  display  a  saturating  effect  as  9** 
increases;  this  relays  a  breakdown  of  single  scattering 
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theory  due  to  the  stronger  interaction  between  the 
incident  wave  and  the  electron  fluctuations.  The 
present  analysis  brings  in  the  effects  of  multiple 
scattering  of  the  incident  wave  and  weighs  their  effects. 

The  results  show  a  linear  relationship  between  the  scattered 

v  * 

incoherent  powers  for  small  mean  square  fluctuations,  o*  j 

ap  the  value  of  O'1,  increases,  the  scattered  power 
begins  to  saturate  (fig. 4  ,5  ).  This  saturation  is 
accentuated  as  the  normalized  plasma  frequency,  , 

nears  the  transition  regime  (fig.  6,7  ).  This  is  expected 
since  the  mean  index  of  refraction  appears  in  the  de¬ 
nominator  of  the  coherent  equation's  effective  index  of 
refraction.  The  saturation  eventually  levels  out, 
tending  generally  to  some  limiting  value;  this  qual¬ 
itative  behavior  is  supported  by  experimental  observations 
(40). 

If  a  slight  damping  is  introduced  weak  and  subdued 

*  '  ’  ’  ^  r  .  i  •  ••  i  .  . 

incoherent  scattering  results  (  fig. 8  ,  9,10),  since  its 
source,  the  coherent  field,  suffers  attenuation  as  it 
illuminates  the  layer.  In  addition,  no  appreciable 
departure  from  the  power  scattering  law  <ISR|>~  a*1 
with  increasing  vx  is  noticed.  However,  it  should  not 
be  surmised  that  the  single  scattering  theory  gives  a 
physical  description  of  the  interaction. 

A  comparison  of  the  homogeneous  fluctuational 
effect  with  the  inhomogeneous  one  show  a  relatively  higher 
incoherent  scattering  in  the  former  case  throughout  the 
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variation  of  the  mean  square  fluctuations  .(  fig.  l£)  ;  a 
less  definite  dependence  appears  in  the  coherent  scattering 
case  (fig.  12,13) .  For  a  locally  homogeneous  randomness, 
the  Gaussian  correlation  shows  a  larger  scattering  power 
than  t{ie  exponential  one  for  the  cases  considered  (  fig.  14  ) 
Typical  examples  for  the  dependence  of  the  coherent 
reflected  and  transmitted  powers  on  the  mean  square 
fluctuations  are  shown  in  Figurel5,16.  The  results  show  a 
definite  dependence  on  large  c*1,  .  This  correlates  with 

the  theoretical  formulations  since  g**  appears  as  an 
attenuation  effect  on  the  coherent  field;  this  attenuation 
describes  analytically  the  energy  transfer  from  the 
coherent  to  the  incoherent  field.  For  slight  damping  the 
variations  in  the  coherent  scattered  power  as  a  function 
of  C*U  are  subdued  in  comparison  to  zero  collision;  for 
the  coherent  power  is  comparatively  insensitive  to 


*»V  ^  ^ 


.  r  i.v 
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The  results  seem  to  indicate  that  a  nonrandom  approach  should 
be  sufficient  to  describe  the  coherent  return  with  an 
accuracy  of  at  least  3  db  . 

It  is  to  be  noted  that  the  present  formulation,  contrary 
to  the  preceding  one  in  Chapter  III,  obeys  the  principle 
of  conservation  of  energy  as  is  shown  in  the  varied  cases 
considered. 
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Incoherent  transmitted  power  as  a  function  of  the 
mean  square  fluctuations  with  normalized  plasma 
frequency  as  a  parameter. 
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Rj  X  10 


Fig.  10  The  incoherent  reflected  power  in  underdense-overdensc 

regime. 


Rj  X  10 
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Comparative  effect  of  correlation  function  on  the 
reflected  power. 


V.  TRANSFORM  METHODS  AND  THE  SMOOTHING  EQUATION 


A.  Introduction 

In  the  preceding  chapter  we  examined  the  scattering 
problem  in  the  limit  of  small  scale  fluctuations.  In  this 
chapter  the  problem  will  be  formulated  as  a  transform¬ 
smoothing  equation  in  order  to  achieve  a  more  general 
description  of  the  fluctuational ' s  scale.  A  singular 
integral  equation  of  the  standard  form  is  obtained. 

Formal  solutions  of  this  equation  are  produced  and 
discussed.  Direct  numerical  simulation  is  also  investi¬ 
gated. 

B.  Analysis 

For  a  plane  incident  wave,  the  integral  formulation 
of  Equation  (11,1)  takes  the  form: 

-  4>U)  =  t  *  +  £  fo  y\(x')  fU')dxf  ^ 


This  equation  may  be  transformed  into  wave  number-space 
upon  multiplying  it  by  *  and  integrating 

over  the  whole  space: 
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Equation  (2  )  may  be  put  in  operational  form 


^K>  '  V*'^  +  €  4>Cp)  (4) 

The  Fourier  transform  of  the  correlation  function  (1,2  ) 
may  be  written  as  (  13) 


<’ICP',’jW  >  =  A(l>,rp.)  tCf.-r,) 


(5) 
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not  reproducible 

where  and  p,  )  are  respectively  the 

Fourier  transform  of  the  intensity  and  the  small  scale 
behavior  of  the  plasma.  The  first  order  smoothing  equation 
(  4  )  for  the  coherent  field  gives  then 


>  =  f(K,  k)  +  i 

where 


$  ^  P  ^  ~  ^  A  (  Jc-  p)  ^  ( i<  +  p  ) 

yu)  -  G0fc  p; 
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'  cuguouxOii  v"*  J  “'oy  u6  Vvixc ceii  explicitly  as 
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Replacing  by  +k  in  equation  (8  ) ,  we  obtain  the 

forward  scattering  amplitude: 

I  fCk.v)  ♦  Cirri'*,)'] 


The  back  scattering  amplitude  is 


1  K&C'L)  7  =  +  L  f-^F 

U I  p-  *•  +jfl 


(-  u,  p)  <  <$>  (p )  >  f  <  #(-?>?) 
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The  sum  and  difference  of  (  9)  and  (lo)  are  respectively 


where 
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C.  Singular  Integral  Equations 

The  application  of  Plemmelj  formulae 


/a  ECpL. 
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to  equation  (11)  leads  to: 
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Equation  (14)  may  be  set  in  normal  form;  for  example, 
the  Kernel  may  be  expanded  in  a  Taylor's  series; 

then  we  have 


i 

p  ,  =  Uh  k2  +  k  c  it,  p) 

k  f-  1c.  (16) 

where 


KCfc,p)  = 

(17) 


manipulation  of  Equation  (14)  and  (17)  leads  to 
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or  symbolically 
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where 
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The  expression 
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is  said  to  be  the  regular  part  of  the  equation. 


(21) 


D.  Solution 


singular  integral 


A  general  method  of  solution  of  the 
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•  / 

equation  has  been  advanced  by  Carleman  and  developed  by 
Vekua  (36).  The  method  ((6)  consists  in  the  separation 
of  the  characteristic  part  from  the  equation  and  in  its 
solution.  This  reduces  the  singular  integral  equation 

to  a  weakly  singular  one.  We  shall  now  consider  the 

• 

characteristic  equation.  Its  solution  is  well  known,  (SB)  • 
We  have 


+(*0  =  -  ilk)  HijUf.  {(?) 

■wj  Vd  i (pi 
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If  the  coefficients  of  Equation  (ij.8)  do  not  satisfy  the 
condition 


'  aHw)  -  4‘C  k)  =  I 


(24) 


we  must  divide  both  sides  of  (18)  by 


[  a?C  U)  .  ir'(u) 
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Applying  Equation  (22)  to  (19),  we  obtain,  after  some 
manipulation. 


^(L)  + 


(p,k) 


&(  p?  dp 


(25) 


where  the  Fredholm  kernel  is  defined  by 


K’^Cbjk)  sacu)  b(k)  lit)  L  <r  Vp,, p ) 


(26) 


and  where  the  free  term  has  the  form 


Tfj  WrV 


(27) 
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The  Carleman-Vekua  method  forms  a  basis  of  the  theory  of 
singular  integral  equations.  In  the  case  the  kernel  of 

to 

Equation  (10)  is  degenerate  /<*(!»,/>)  s  %.  of.  (It)  &• ([) 

hi  '  ' 

the  solution  of  Equation  (18)  is  given  in  closed  form. 


(fCk):  U,(f)  f(p1dp  (28) 

/  * 

This  is  a  Fredholm  equation  with  degenerate  kernel.  It 
may  be  reduced  to  an  algebraic  system  of  equations  (2>4)  . 

In  fact,  a  singular  integral  equation  of  general  type 
may  be  reduced  to  an  equation  :with;  degenerate  kernel. 

By  virtue  of  Weiertrass’  theorm,  kw(lt,p)  can  be 
approximated  by  a  degenerate  kernel 


K'U.p)  =  +  'fCk.p') 


(29) 


Where  the  norm  ||  y(  k,  p)  }f  in  the  space  i-t  may  be  made 
as  small  as  we  please. 

In  the  case  where  the  random  intensity  varies, slowly 
in  real  space  with  respect  to  the  coherent  '.field,  a 
simplification  of  equation  (6  )  results?  the  coherent 
field  under  the  integral  sign  is  considered  to  be 
sampled  at  the  wave  number  k  ?  equation  (6  )  reduces  to 
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<  <K»0  >  5  -6  Attc.k)  A  ) 
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\ 

The  incoherent  field  is  then 

$4>  z  f  £  G9f(pf)  <<t>(f,)  > 

1  (  (31) 


We  have  concentrated  so  far  on  the  analytical 
methods  of  solutions.  Direct  numerical  simulation  of 
Equation  (n)  should  not  be  discarded;  the  apparent 
difficulty  is  due  to  the  singularity  in  the  kernel;  it 
may  be  avoided  through  the  introduction  of  an  appropriate 
average  value  for  the  inteqral  in  the  vicinity  of  o=k 
or  the  rewriting  of  the  integral  equation  so  that  the 
singular  contribution  cancels  out. 

E.  Conclusion 

The  present  development  appears  to  have  some  definite 
potentials.  A  survey  of  the  methods  for  the  approximate 
solution  of  the  resulting  equation  is  given.  For  slowly 
changing  characteristics  of  the  random  intensity,  the  scat¬ 
tered  field  is  easily  deduced.  In  addition,  the 
formulation  is  close  to  reality  since  it  includes  the 
.  effect  of  inhomogeneous  turbulent  intensity. 
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N0T  REPRODUCIBLE 

VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  Conclusions 

This  thesis  treats  the  case  of  a  random  slab  whose 
mean  refractive  index  is  different  from  the  deterministic 
refractive  index  in  the  exterior.  Previous  studies 
neglected  the  effect  of  the  boundaries  in  the  computation 
of  the  mean  wave;  this  procedure  is  shown  to  be  ener¬ 
getically  inconsistent.  The  present  treatment  is 
consistent  and  evaluates  the  coherent  and  incoherent 
scattered  powers.  The  principal  conclusions  of  the 
analysis — that  the  incoherent  scattered  powers  experience 
saturation  with  a  further  increase  in  the  mean  square 
fluctuations,  and  that  the  coherent  scattered  power  may 
be  adequately  predicted  by  a  nonrandom  approach  for 

graphically  for  some  typical  cases. 

B.  Recommendations 

1.  The  present  work  has  concentrated  on  the  limit  of 
small  scale  fluctuations;  a  measure  of  the  sensitivity  of 
the  power  scattered  as  a  function  of  the  correlation  scale 
is  desirable.  A  proper  formulation  has  been  set  in  Chapter 
V.  of  this  thesis  and  can  be  used  to  investigate  this 
point. 

2.  A  comparable  analysis  to  the  present  work  is 


needed  for  three  dimensional  fluctuations.  One  would 
still  assume  that  the  correlation  function  varies  much 
more  rapidly  than  the  Green* s  function,  and  it  may  be 
viewedias  an  impulse  function  except  in  the  case  where 
divergent  integrals  are  obtained  if  the  unsmeared  function 
Is  used. 

3.  The  present  method  of  solution  is  not  restricted 
to  a  homogeneous  background  layer  with  sharp  boundaries. 

It  may  be  used  in  the  case  of  layers —  Epstein,  linear, 
parabolic  profiles —  which  have  inhomogeneous  backgrounds 
and  diffuse  boundaries.  It  would  be  interesting  to 
study  the  transitional  effect  on  the  scattered  power. 
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APPENDICES 


APPENDIX  A 


ON  THE  VALIDITY  OF  THE  BORN  APPROXIMATION  IN  RANDOM 

WAVE  SCATTERING 

% 

A.  Introduction 

.  The  scattering  of  a  plane  wave  incident  upon  an 
appropriate  random  medium  has  been  widely  treated  in  the 
first  order  Born  approximation^ Evidently,  there  are  two 
points  of  important  consideration:  1.  For  a  given 
scattering  problem,  does  the  Born-Neumann  series  converge? 
2.  If  this  series  converges,  how  large  and  error 
is  incurred  by  truncating  it  after  the  first  term?  In 
this  chapter,  a  simple  quantitative  criterion  is  derived 
that  insures  convergence  and  gives  an  upper  bound  for  the 
truncation  error.  This  bound  is  expressed  in  terms  of 
relevant  parameters  which  characterize  the  convergence 
condition,  the  scattering  volume,  and  the  structure  of  the 
medium.  This  result  is  established  using  the  method  cf 
successive  approximation,  and  the  inequalities  appropriate 
to  the  condition  of  quadratic  summability  of  the  kernel 
in  the  basic  mean  square. 

B.  Statement  of  the  Problem  and  Solution 

A  scattering  problem  may  be  represented  by  the  follow¬ 
ing  integral  equation  (  5  7  )  • 
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=  <t=Jn.)  y  £  JSfx,  xV  (1) 


1 

A  straight  forward  approach  to  solving  the  integral 
equation  is  successive  approximations.  We  begin  with  the 
zero  order  approximation 


(2) 


Let  us  substitute  the  zero  order  approximation  into  the 
original  equation  under  the  integral  sign  to  obtain  a 
first  order  approximation,  and  the  process  is  then  repeated. 
The  resulting  series  for  the  scattered  field  is 


Where 


<j>^(  x)  =  t’"  j  K,  U,  V)  <t>  ( l')  di! 


*  j \ \ )  KS\ < V d\ 


(4) 
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K'ln-X)  -  c  i &(*!**')  ^(h!) 

and  KJl,  n-V  is  called  the  iterated  kernel  in 
relation  to  the  given  kernel.  The  first  term  in  the  series 
solution  may  be  regarded  as  contribution  due  to  the  incident 
wave;  the  second  term,  the  contribution  due  to  a  single 
scattering  of  the  wave  from  each  inhomogeneity  in  the 
scatterer;  the  third  term,  the  contribution  caused  by 
the  double  scattering  of  the  wave,  etc.  The  formal 
integration  of  each  term  in  the  series,  though  random, 
should  amount  to  the  integration  of  functions  when  a 
particular  sample  function  of  the  ensemble  is  consideredCt.4). 
To  each  sample  function  corresponds  a  solution 
ani*  p  «■»  i  i<ai  A  A  tiv , "’r.lv.C  of  the 

solution  is  generally  different  for  each  sample  function, and 
the  solution  thus  takes  a  whole  set  of  values.  The  ensemble 
averages  of  these  values  are  of  interest  in  the  applications. 
The  method  consists  of  the  following  steps:  The  explicit 
solutionC3)of  (1)  is  first  determined  formally  for  each 
member  of  the  ensemble.  Then,  the  mean  and  the  correlation 
function  of  fyfy) «*©  obtained.  As  the  integration 
becomes  exceedingly  complicated,  calculations  have  been 
limited  to  the  first  term.  A  bound  on  the  truncation  error 
of  the  series  and  the  condition  of  convergence  or  that  series 


A 
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are  investigated. 

C.  Convergence  and  Error  Bounds 

The  general  term  in  the  series  solution  for  the 
mean  square  value  of  the  scattered  field  may  be  written  as 


= <1  Ik  iuj 


The  Bunyahovshi  inequality  states  that 


(7) 


l(4,4>  I  « H  41M1411 


(8) 


Using  the  above  relation,  we  get  for  the  second  order 

UO*  pCu,  g-JcL*.  J  £  <“-X?  <,  (9) 

Ws  assume  that  the  random  variable  >|  follows  a  Gaussian 
process  with  mean  zero  and  rms.  Applying  the  relations  in 
(8  )  and  (9  Jtythe  mean  square  value  of  (7  )  ,  we  obtain: 

(10) 


<M>j\  <  ['*)  o'-^/lV/i^wiV 


0 

**  / 
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-  Sur  M 


\)  Wt, ,  i')«k 


Then 


c„  =  s 


«  I  l  ( i  «■■*./*, S )  l*«u.  fn«'„'x';lk  J^' 


<  B 


where 


b'  = 


J I IK(i,y)  tXdvt  ifi! 


Using  the  recurrence  relation  in  ( 13) ,  there  results 
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5 ' "" 5 


(15) 


* 

Since  the  scatterer  has  a  finite  volume,  V,  a 
relation  exists  between  B  and  C,  where 


b1  «  c,v 

(16) 


and  (15)  becomes 


•  I 


(17) 


Poplar (17)  iri  '(10) 


th^Te 


Where 

T>  =  C,  vV 


(18) 


(19) 


The  series  is  uniquely  convergent  for  D  less  than  one 
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Now  we  establish  an  error  bound  on  the  truncation  of 
the  series  after  the  first  term,  thus  estimating  the 
importance  of  multiple  scattering,  in  view  of  the  elementary 
inequality 


(  4;  f  <^1*  «  3  c  W’l'  *■  I*.  1*7 


(2p) 


the  scattered  pov;er  is: 


<i4>K)i%  .<  i  ^Cr)* 


(21) 


A  comparison  of  the  first  Born  approximation  with  {  20) 
gives  a  measure  of  the  relevance  of  multiple  scattering. 
Admitted!'/,  thjj;  mr*af’v.Tn  5**  cn»nc»T.,Vi5>4-  the  ether 

hand,  a  very  stringent  bound  would  not  satisfy  many 
practical  cases. 
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•  APPENDIX  B 

AD  IIOC  DERIVATIONS  OF  FIRST  ORDER  SMOOTHING  EQUATION 
\ 

A.  Formulation 

•  For  wave  propagation  in  random  media,  the  first  order 
smoothing  equation  has  been  derived  in  different  ways  (  5  ) 
(15),  (a6) .  The  scattering  equation  in  operational  form  is 

C-.  +  t.)  <|>  5  O  - 


The  mean  is 


L.  <*>  +<L,H>  *  o  (a) 

The  difference  between  1  and  2 

L.  S  4>  r  l,  <  >  +  f  $  (3) 


where 


s  L, 


•  # 

•  / 


■  Equations  2  and  3  form  a  set  of  coupled  equations  that 
describe  the  problem. 

B.  Local  Independence 

With  the  hypothesis  of  local  independence,  the  wave, 
?nd  the  perturbation, L, ,  are  considered  statistically 
independent,  i.e. 

<L,L*L.  4>  >  =  -ft-,  L,  >  <  <#>>  (4) 


Let  us  multiply  (1)'  byL1  and  average.  The  resulting 
equation  may  bo  further  simplified  by  applying  the  local 
independence  hypothesis;  there  results 


L.4  ♦  » 


-  <  L,  r.  L,  >  <  J.  > 


The  above  technique  leads  to  the  same  formulation  as 
the  first  smoothing  equation.  A  physical-  measure  to  the 
assumption  of  local  independence  may  be  given.  It  is 
known  that  sharp  changes  in  the  refractive  index  makes 
a  wave  dependent  upon  the  medium;  then,  a  necessary 
condition  for  the  hypothesis  to  hold  requires  that  the 
dependence  or  correlation  between  the  medium's  index  of 
refraction  and  the  wave  be  zero.  In  this  case,  there  are 
'  'reasons  to  believe  that  the  functions  are  statistically 


independent,  or  at  least  the  stochastic  dependence 
between  them  is  so  weak  as  to  be  negligible. 


C.  "Distorted"  Wave  Born  Approximation 

In  the  so  called  "distorted"  wave  Born  approximation, 
the  first  smoothing  equations  are  derived  by  neglecting 
the  term  in (2)  ,  with  no  apparent  justification. 
Manipulation  of  ^  and  result  in 

-i 

L.  <  4-  >  -<L,  L,  L,  ><<>>  =  o  (6) 

$  $  ~  ~  (7) 

which  are  exactly  the  same  as  the  first  smoothing  equation. 
Fof'h  1  'Yher  orders/  however,  th.*  tn  rnQ+-vr»dc 

generally  differ.  This  is  expected  since  the  neglected 
term  generates  terms  comparable  in  size  to  the  remain¬ 
ing  ones  in (2) 

D.  Modified  Neumann  Method 

The  first  smoothing  equations  are  here  derived 
directly  from  the  Neumann  series.  let  the  field  £  be 
replaced  in(l)  by  the  incident  field  ^  and  scattered  field 


Solving  for  ^  ,  we  get 


i  =  m  l:’  l,  <f>, 


(9) 


inhere 


M  =  ( 1  ♦  O, )"  - 1 +C  i.,  +  a;'Lj^ ... 


do)  * 


To  first  order  approximation  we  get  for  the  mean  field 


*-•  <  £  >  -  <  L,  L,  L  )  <t>m  -  o 


(ID 


^  .  jn  front  of  the  oner? tor,  is  the_first; 

iteration  of 


<t>  (b 

To  ^ 


t.  <  ^  >  -  <  t  C  t  >  <  ^  >  s  O 


(12) 


The  fluctuation  field  is  easily  solved  for  ,  as  : 
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1.  INTRODUCTION 


In  contemporary  communication  processes  (e.g., 
ionospheric  and  tropospheric  scatter  links;  stellar 
observations;  re-entry  wake  analysis),  electromagnetic 
scattering  plays  an  important  role  and  has  received  interest 
since  the  very  earliest  formulation  of  electromagnetic 
radiation  theory.  The  influence  of  matter  upon  propagating 
electromagnetic  signals  can  be  classified  into  separate 
topics,  such  as,  scattering,  absorption,  depolarization, 
dispersion  [e.g.,  Jor.ee,  1964].  In  this  work,  we  investigat 
the  scattering  of  electromagnetic  waves  by  a  non-uniform 
plasma  with  random  irregularities.  Consideration  will  be 
restricted  to  monochromatic  waves,  which  will  be  allowed  to 
interact  with  a  plasma  medium  characterized  by  the  Lorentz 
model  [Hcald  and  Wharton,  1965],  Such  a  model  has  been 
shown  to  be  valid  for  radio  and  radar  propagation  in  plasmas 
associated  with  rocket  exhausts,  re-entry  wakes,  and  the 
ionosphere,  in  the  absence  of  a  magnetic  field  [Jarem,  1969] 
The  Lorentz  plasma  model  leads  to  an  effective  dielectric 
constant  for  the  plasma  which  will  be  taken  to  be  a  random 
continum,  i.e.,  the  particle  character  of  the  media 
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(electrons,  ions,  neutral  molecules)  will  be  ignored  and  the 
random  deviation  of  the  dielectric  constant  from  the  mean 
(due  to  turbulence)  will  be  continuous,  varying  smoothly  in 
time  and  linear  extent.  Complete  information  about  the  media 
will  never  be  available  and  thus  a  complete  description  of 
the  interaction  with  the  electromagnetic  waves  must  be  made 
in  terms  of  a  statistical  interpretation  [e.g.,  dnWolf,  1968]. 
Such  media  are  mathematical  models  of  real-life  situations  » 
namely,  turbulent  atmospheric  layers,  rocket  exhausts,  and 
re-entry  wakes.  A  good  deal  of  interest  in  these  media  comes 
from  the  irregular  disturbances  of  radio  and  radar  communi¬ 
cation  signals  when  propagating  through  such  physical  media. 
Such  disturbed  signals  may  also  be  used  as  remote  sensors 
for  basic  information  on  turbulent  flows,  and  as  a 
diagnostic  tool  for  re-entry  wake  analysis. 

Solutions  to  the  problem  defined  above  are 
already  known  in  certain  limiting  cases  (e.g.,  low  mean 
electron  density,  and  small  turbulent  fluctuation  levels; 
particularly  simple  geometries).  It  is  the  aim  of  this 
paper  to  study  the  formulated  problem  for  the  one-dimensional 
Epstein  profiles  [Epstein,  1930]  for  an  arbitrary  peak  level 
of  electron  density.  Electron  density  fluctuation  levels 
will  not  be  arbitrary  but  will  be  restricted  by  the 
perturbation  method  here  employed.  Where  applicable,  the 


results  will  be  compared  to  the  statistical  computer 
experiments  of  Hochstin  and  y.avter.a  [1969;  1970],  vrttc  hav 
examined  the  case  of  a  randomly  fluctuating  plasma  slab. 
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2.  BACKGROUND 


2.1  BORN  APPROXIMATION 

Scattering  from  turbulent  media  has  drawn  in¬ 
creasing  interest  since  the  publication  of  work  on 
tropospheric  scattering  by  Booker  and  Gordon  [1950] .  Booker 
and  Gordon  made  use  of  the  perturbational  Born  approximation 
[e.g.,  Dioke  and  Wittke ,  I960],  and  since  that  time,  a 
number  of  papers  have  been  written  to  extend  the  accuracy 
of  that  approximation,  many  of  them  attempting  to  treat 
the  full  three-dimensional  electromagnetic  problem  [e.g., 
Tatar$kHt  1961;  Salpeter  and  Treiman ,  1964;  Menkes ,  1964; 
Ruffine  and  deWolf ,  1965].  Because  of  the  perturbational 
nature  of  the  Born  method,  solution  above  the  lowest  order 
yields  only  slight  additional  information  with  a  great 
increase  in  complexity.* 


*An  exception  to  this  statement  might  be  made  for  the  second 
Born  approximation,  which  yields  cross  polarized  information 
predicted  to  be  zero  by  the  first  Born  approximation 
[Ruffine  and  deWolf,  1965;  Jarcn ,  1969b].  This  information 
still  comes  with  great  complexity. 
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To  apply  the  Born  approximation  we  start  with  a 
wave  equation  for  the  propagation  in  a  turbulent  medium 
(This  equation  will  be  derived  for  a  special  geometry  in 
Section  5.2.). 


♦  kl  [<o  ♦  *ifs)  +  <5 *(•-.)]  >6  =  0 


(1) 


where  the  wave  number  has  been  separated  into  a  constant 
part  k0,  an  inhomogeneous  sure  part  k1  and  a  random  in* 
homogeneous  part  5k.  We  then  write  (1)  with  the  homogeneous 
term  on  the  left 


id* /dr.* 


♦  «<]<{> 


(2) 


Both  inhomogeneous  terms  are  thus  considered  as  source 
terms.  A  free  apace  Green's  function  [Friedman,  1956]  can 
thus  be  utilized  in  the  integral  equation  corresponding 
to  (2) 

00 

*(•-}  *  ■>-■(-)  *  /  /•(s' )*(•’.« ')<?«'  . 


(3) 


7 


where  /(a')  represents  the  right  side  of  (2)  and  <{>.(2),  the 
solution  of  (2)  with  /  ■  0.  An  approximate  solution  to  the 
integral  equation  (3)  can  then  be  obtained  by  iteration  of 
(3)  one  or  more  times.*  The  accuracy  of  such  a  solution 
depends  on  the  size  of  the  entire  perturbation  <^6*. 

Thus,  even  for  small  random  deviations,  one  is 
restricted  by  the  variation  of  from  zero.  For  a  plasma 
medium,  this  restriction  limits  solution  to  the  underdense 
regime  (See  Section  4.).  Even  when  kj  is  constant,  one  is 
restricted  to  a  finite  size  plasma-!*,  because  the  formulation 
of  the  Born  approximation  docs  not  allow  for  attenuation  of 
the  incident  waves. 

2.2  THE  DISTORTED  WAVE  BORN  APPROXIMATION 

Various  formulations  categorized  as  distorted 
wave  Born  approximations  have  been  devised  to  circumvent 
some  limitations  of  the  Born  approximation  [Krsoa,  1968]. 

The  easiest  to  apply  involves  a  local  solution  of  the  wave 


•The  method  becomes  very  unwieldy  for  more  than  two 
iterations. 

tA  propagation  problem  is  thus  excluded. 

SThus  there  are  an  infinite  number  of  scattering  events  of 
equal  importance,  with  a  corresponding  blow-up  in  the 
scattering  coefficients. 


$ 


equation  written  as  follows 


{d2/dr.2  ♦  +  KjC^OIH  *  (4) 

The  Green's  function  is  taken  to  be  of  the  free  space  form, 
but  the  nonrandom  portion  of  the  wave  number  is  allowed  to 
vary  in  accordance  with  its  actual  distribution.  This 
formulation  fails  when  the  plasma  wave  number  has  rapid 
variation  within  the  distance  of  a  wavelength,  for  a  local 
solution  docs  not  include  the  effects  of  the  rate  of 
variation  of  the  inhomogeneity . 

An  exact  formulation  of  the  problem  is  obtained 
through  the  use  of  the  correct  Green's  function  for  (4). 

Such  a  Green's  function  is  available  for  only  simple  plasma 
profiles  such  as  a  slab  or  halfspace  (or  as  we  later  show, 
the  Epstein  profiles).  The  resulting  integral  equation  is 
the  same  as  (3)  except  that  the  Green's  function,  G  is 
correct  for  the  particular  profile  at  hand,  and  f  is  now 
the  right  side  of  (4).  Thus  the  random  inhomogeneity 

• 

represents  a  perturbation  about  the  mean  profile.  Since 
this  perturbation  is  often  small,  a  meaningful  solution  to 
(4)  can  be  obtained  even  when  the  plasma  varies  greatly 
from  free  space  conditions  (The  Green's  function  represents 
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an  exact  solution  when  there  is  no  random  term.)*  The 
difficulty  introduced  by  the  formulation,  is  that  the 
Green's  function  is  no  longer  of  the  convolution  form  as 
was  the  case  of  the  free  space  Green's  function.  Thus, 
Fourier  transformation  (which  is  sometimes  useful  in  formu¬ 
lations  using  a  free  space  Green's  function")  does  not  offer 
any  apparent  benefits  (for  a  non-convolution  form  Green's 
function  problem). 

2.3  SELECTIVE  SUMMATION  TECHNIQUES 

A  number  of  techniques  have  been  devised  to 
overcome  the  limitations  of  the  successive  type  approxi¬ 
mations,  including  smoothing  [Tatarskii  and  Gertsenehtein , 
1963],  renormalization  [Karal  and  Keller,  1964]  and  diagram 
methods  [Bourret ,  1962].  To  the  lowest  solvable  order,  all 
these  techniques  are  equivalent  [Jarem,  1970].  We  shall 
outline  the  smoothing  method,  and  we  refer  the  reader  to 
the  discussion  given  in  a  summary  paper  [e.g.,  Frisoh ,  1968; 
Kreea ,  1968]  for  further  details  of  these  solutions. 

In  the  method  of  smoothing,  we  iterate  in  a 
manner  to  obtain  the  fluctuating  field  in  terms  of  the 


"A  number  of  applications  of  Fourier  transformation  have 
been  given  by  Jarem  [1968,  1969b]. 
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mean  field,  and  an  equation  for  the  mean  field.  We  may 
write  the  scattering  integral  equation  (3)  in  operator 
format  as 


$  *  (5) 

m 

where  a  represents  By  sepaiately  applying  the 

averaging  operator  P  and  the  fluctuation  operator  I-P  (I 
being  the  identity  operator)  to  (S) ,  we  obtain  a  pair  of 
coupled  integral  equations  for  the  coherent  (+Q)  and  the 
fluctuating  (4^)  fields 


F*  :  <*0  ■  ^  ♦  GPO*!  (A) 

(7-P)*  =  ^  -  <7(J-?)o(<>0  ♦  ^).  (ro 

Formal  iteration  of  (6B)  yields 


l 

r.=  1 


0‘ 


(7) 
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Using  this  expression  in  (6A) ,  we  obtain  a  master  equation 
for  the  coherent  field 

m 

*0  ‘  \  +  l  Sro[C7(I-P)c]\c.  (8) 

nm  l 


By  retaining  only  the  first  term  of  the  summation  in  (8),  wc 

* 

obtain  the  first  order  smoothing  equation 


+  o  *  +  JPc5oV 


which  by  iteration  is  equivalent  to 


*0  *  +  l  lGPoGa]n<f<i. 

n*  l 


CP) 


(10) 


Therefore  the  first  order  smoothing  solution  represents  a 

summation  of  terms  which  are  indicative  of  multiple 

scattering  and  are  thus  included  in  $0  obtained  by  solving 

the  integral  equation  (9).  The  fluctuating  field  $  can 

now  be  obtained  to  first  smoothing  order  by  using  (7),  with 

»  ■  1.  The  first  order  smoothing  equation  has  been 

determined  by  Friaoh  [1968]  to  be  valid  for  (Jc.an  )2  <<  1, 

v  P 
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where  a  is  the  correlation  length  (defined  in  Section  5.2), 
and  ftp  is  the  normalized  plasma  frequency  (defined  in 
Section  4).  The  integro-differential  equation  obtained  by 
applying  a  second  order  differential  operator  to  (9)  has 
been  solved  only  approximately,  even  for  facile  geometries 
[e.g.,  Kupicc  et  al.,  1969].  The  most  simplifying 
approximation  to  apply  is  the  assumption  of  a  delta  function 
correlation.  This  assumption  immediately  casts  the  equation 
for  the  coherent  field  into  an  ordinary  differential 
equation. 

2.4  OTHER  SOLUTION  METHODS 

In  addition  to  the  two  categories  of  solution 
techniques  described  above  some  work  has  been  carried  out 
by  energy  transport  theory  and  also  by  information  theory 
concepts  [Kreea,  1968],  The  energy  transport  method 
[ Watson ,  196$]  considers  scattering  by  individual  scatterers, 
and  results  in  a  set  of  coupled  multiple  scattering  equations. 

In  the  information  theory  formulation  [ Kreca , 

1968]  ,  an  attempt  is  rade  to  obtain  closure  of  the 
stochastically  nonlinear  wave  propagation  equation  without 
use  of  a  perturbation  expansion  or  invoking  any  closure 
assumptions.  This  method  has  been  useful  in  statistical 
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3.  STATEMENT  OF  PROBLEM 

In  recent  years,  research  work  in  random  scattering 
has  been  carried  out  for  one-dimensional  plasmas,  not 
necessarily  because  there  is  a  correspondence  to  a  physical 
problem,  but  because  insight  is  sought  for  the  more 
realistic  three-dimensional  geometry.  It  is  in  this  spirit 
that  we  shall  confine  our  investigation  to  a  one-dimensional 
problem,  even  though  there  is  a  direct  correspondence  to  the 
physical  situation  of  reflection  of  vertically  directed 
radio  waves  from  a  turbulent  ionosphere.  The  analysis  may 
in  fact,  be  easily  extended  to  oblique  incidence,  but  will 
not  be  considered  here,  as  it  would  only  introduce  one 
additional  variable  into  a  result  whose  analysis  is  already 
burdened  by  a  large  number  of  parameters.  Propagation  will 
thus  take  place  in  a  layered- inhomogeneous  medium  whose 
properties  only  vary  along  one  axis  of  a  rectangular 
coordinate  system 

Kc  will  consider  both  wave  reflection  and 
transmission  by  the  random  inhomogeneous  layers,  and  to 
obtain  a  rigorous  solution  of  this  problem,  we  solve  the 
one-dimensional  scalar  Helmholtz  equation,  subject  to  the 
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appropriate  boundary  conditions.  For  a  wave  of  unit 
magnitude  incident  from  the  left  these  are  [Stratton ,  1941, 
Chap..  9] 

2  ->  -^0"  (11) 

2  -  -,  *(s)  -  (12) 

where  R  and  T  are  the  reflection  and  transmission 
coefficients,  respectively,  which  are  to  be  determined. 

Because  of  the  difficulties  involved  in  using  the 
smoothing  equation,  we  have  chosen  to  initially  investigate 
the  problem  of  a  very  inhomogeneous  and  random  plasma  by 
the  perturbation  method  with  an  exact  Green’s  function.  In 
addition  to  providing  results  for  a  restricted  parameter 
range,  such  a  technique  should  yield  insight  into  solving 
one  of  the  wider  range  formulations.  We  will  determine  the 
scattering  coefficients  for  a  particular  distribution  of  mean 
layer  properties  (discussed  in  Section  5),  and  for  particular 
values  of  the  relevant  stochastic  parameters. 

As  a  means  of  providing  confidence  in  our  result, 
we  will  compare  our  work  to  the  computer  experiments  of 
Bochstim  and  Martens  [1969,  1970]  for  an  appropriate 
selection  of  parameters. 
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Although  we  have  chosen  to  examine  the  case  of  a 
plas  ma  medium,  the  analysis  may  be  easily  applied  to  the 
propagation  of  acoustic  waves  in  layered  turbulent  media, 
or  the  propagation  of  electromagnetic  waves  in  a  randomly 
fluctuating  real  dielectric,  such  as  might  represent  the 
troposphere . 
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4.  ELECTROMAGNETIC  WAVE  PROPAGATION  IN  PLASMAS 

The  governing  equations  for  electromagnetic  wave 
propagation  arc  the  Maxwell  equations,  [c.g.,  Stratton,  1941, 
Chap.  1;  Ramo  et  al.  ,  1965] 

v  x  r  «  -U/n  (A)  v  •  D  •  p  (C) 

(13) 

v  *  Ti  *  7  *  o n/:t  (ii)  v  •  b  *  o.  (D) 

By  simple  manipulation,  the  Maxwell  equations  can  be  reduced 
to  a  vector  Helmholtz  equation  for  the  electric  field. 

First,  by  assuming  a  time  dependence  of  the  form  e"*a,t ,  and 
taking  the  curl  of  (13A) ,  we  obtain 

7  x  v  x  z  *  £wn0(o  - 

(14) 

■  w2m0(c0  +  ia  I  w)  L' . 

An  effective  dielectric  constant  [Papas,  1965]  can  now  be 
defined  by 

“  co 


c 


+  i c/u . 


(15) 
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liquation  (1*1)  now  becomes 


7*7  x  r 


(10) 


From  the  standard  vector  identity  [c.g.,  Brand,  1957] 


7  x  V  x  /,  «=  7(7  •  ,T) 


r* 

V 


(17) 


liquation  (16)  bccoacs 


7  (7 


u,2n0c 


(IS) 


The  field  vectors  D  and  K  are  related  by  tlie  constitutive 
relation  [5;vra  i  ton ,  1941,  Chap.  1] 


D  * 


cr. 


(19) 


Equation  (15C)  now  becomes,  for  a  charge  free  region* 


v 


0, 


(20) 


*A11  effects  of  plasma  free  charges  arc  contained  in  the 
conductivity  and,  subsequently ,  in  the  effective  dielectric 
constant. 
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or 


Ve  •  E  *  cv  •  £  «  0. 


(21) 


Thus  wc  obtain 


r/c. 


(22) 


Equation  (18)  is  now  written 


V2.V  +  u)2u0c/,'  *  -7[vc  •  i./ c ]  . 


(23) 


The  electromagnetic  wave  number  k  is  defined  by, 


*2 


0C  , 


(24) 


so  that  (23)  now  becomes 


(V2  +  k*)E 


-V[\'k  •  L/k  ]  . 


(25) 


liquation  (25)  is  the  inhomogeneous  vector  Helmholtz  equation. 
For  plane  wave  propagation  perpendicular  to  the  refractive 
index  gradient 

•  T  =  0,  (26) 

and  any  component  of  1  may  be  represented  by  tire  scalar 
Helmholtz  equation 


(V2  ♦  ;:2 ) 6  *  0.  (27) 

The  Lorentz  model  of  the  plasma  conductivity  may  be  derived 
by  considering  the  motion  of  a  charged  particle  in  an 
electric  field*  [licald  and  Vhavtor.,  1965).  A  one-dimensional 
equation  for  the  forces  on  such  a  particle  can  be  written 
as 


ns:  *  -cl 


v  ns 


(28) 


*A  more  exact  derivation  of  the  plasma  conductivity  lias  been 
given  by  Saver,  [1969],  and  includes  a  discussion  of  the 
validity  of  the  Lorentz  model  for  radar  applications. 
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The  second  term  on  the  right  hand  side  corresponds  to 
viscous  damping  by  the  electron-neutral  particle  collisions. 
The  steady  state  solution  to  (28)  for  oscillatory  fields 
is  given  by  [e.g.,  Rainville ,  1964  ,  Chap.  8] 

x  *  eT’/r,(i\  (w  +  i  v)  .  (2P) 


The  current  density,  J,  is  defined  by  [Stratton ,  1941, 

Chap.  1] 

J  *  or:  *=  -nex  (30) 

The  conductivity,  o,  is  then 


o  ■  ne2/n(v  -  iu)  ■  (ne2/r;)[(v  +  i u)/(v2  ♦  w2)],  (31) 


The  effective  dielectric  constant,  defined  by  (IS),  becomes 


(w?/u)2  (Wp/«)2 (v/u) 

1  +  (v/w)2  1  ♦  (v/w)2 


e 


1 


(32) 
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where  the  angular  plasma  frequency,  is  defined  by 
[Kclco ,  1064] 


nc2/;.?c0 . 


(33) 


A  normalized  plasma  frequency  a  ,  and  collision  frequency 

a  arc  defineJ  as  follows 
c 

nr  *  ur/u  ’  no  -  */« .  (34) 

where  u  is  the  angular  wave  frequency.  Plasmas  for  which 
flp  <1  are  usually  referred  to  as  being  underdcnce  ,  while 
for  which  Q  >1  arc  called  ovcrdcnce . 


those 


25 


5.  THE  EPSTEIN  PROFILE 


5.1  THE  DETERMINISTIC  LAYER 


In  this  section,  we  shall  basically  follow  the 
analysis  outlined  by  Brekhovakikh  [I960].  The  aim  of  this 
derivation  is  to  transform  the  wave  equation  (27),  into  an 
equation  whose  solution  is  already  known.  As  will  be  seen, 
when  the  inhomogeneous  distribution  is  that  which  was 
studied  by  Epstein  [1930],  this  transformed  equation  is 
the  hypergcomctric  equation  [Whittaker  and  Watcon,  1965], 
The  logic  of  our  derivation,  however,  will  be  deductive, 
since  we  will  start  with  the  hyporgeometric  equation,  and 
show  how  this  transforms  into  the  wave  equation. 


The  hypergeomctric  equation  is  given  by 


djF  (q+g+1) g  ■  y  dF  o3 

d(*  '  ”1(1-0  dt  '  5(1*0  " 


(35) 


Let  transformations  of  the  independent,  and  dependent 
variables  be  given  by 

F  *  r  (;;)"(;:) 

5  *  ?(s). 


(56j 
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The  derivatives  of  F  occurring  in  (35)  arc  then 


dF 

dr’ 


nr  „ 

j~  “ 

dt. 


d2r  _  rffr  „  +  ,  df"  dZ 
Sp'  JCT  *  “  dZ  dr. 


d2z 
d?  * 


(37) 


Substituting  these  into  (35),  we  obtain 


+ 


dz 

Jr. 


}.cr  ( a+J  ♦  l *  +  yr_  d£ 

c.*?  (1-0  d?  rfl-O  dZm 


+ 


( a  +  P_+ 1 )  +  _  y _ dr  aP.r 

" *(T*  O  '  Iff.  ; 0  O )  dz  Oi-O. 


P. 


(38) 


In  order  that  (35)  be  equivalent  to  the  wave  equation,  the 
coefficient  of  dF/dZ  must  vanish,  or 


dv  <?.*:  +  f.'2;:  (ct+  ?■+ 1  )r  dr.  +  yr _  dr.  m  „ 

‘  dz  dZ  rdF  '  "(1-0  '  <*5  5(1-0  *5 


(59) 


I-quation  (39)  is  a  separable  differential  equation  of  the 
variables  r,  and  z »  nd  is  readily  integrated  to  obtain  the 
trails  formation 


i*  - 


(c;/‘”)  U 


.  ri  (y-ct-  fi- 1  )/2  „  -  y  / 2 


(40) 


The  function  P(a)  has  not  yet  been  specified,  and  is,  in 
fact,  arbitrary.  The  wave  equation  is  then  written 

d2Z/dz 2  ♦  k2{z)Z  -  0,  (41) 


and  now  the  function  k2(z)  is  given  by 


where 

4*1  *  y(v") 

4 ::2  *  i*(a*a)2  ♦  y(y-2)  (43) 

4 K3  *  (a+?-Y)2  *  1. 

In  order  to  make  the  analysis  tractable,  Epetein 


[1930]  chose  for  the  transformation  of  the  independent 
variable,  the  simple  function 


2^ 


The  dielectric  distribution  in  (42)  then  becomes 


(“) 


(Uc*'  •) 


(45) 


where  the  three  new  constants  kQ,  .*/,  i!  have  been  defined  in 
terns  of  .’.j  ,  ■'  ,  ,"3  as 


*  1/4).  i:7  =  *(;:0/r): 


3  =  -4(7:0/k)2.V.  (46) 


In  order  to  further  examine  the  dielectric  distribution,  a 
layer  thickness  parameter  5  will  be  defined  as 


S  *  2k0/n.  (47) 

To  determine  the  parameters  a,  £,  y ,  the  equations 
(45)  arc  first  solved  for  these  in  terms  of  the  constants 
*ii  *2*  A3-  The  parameters  a ,  B ,  y  arc  then  related  to  the 
inhomogeneous  distribution  through  (46)  and  (47)  to  obtain 

a  *  ( 1  ♦  1 1  -  *  U’[l  -  (1  -")•']}/: 

s  *  (1  ♦  [i  -  i r2  ,j-  ♦  i5[i  ♦  (l -:.•)■*  3 }/:  (4S) 

y  =  !  ♦  iC. 
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The  distribution  of  (45)  is  most  useful  when  cither  U  or  U 
is  taken  to  be  zero.  In  the  former  case,  the  distribution 
is  called  a  symmetric  distribution  since  it  represents  a 
space  filled  with  dielectric  symmetrically  distributed  about 
z  *  0,  with  free  space  condition  as 

The  ease  li  =  0,  represents  a  transition  layer  going 
from  free  space  conditions  at  minus  infinity  to  an  arbitrary 
value  of  dielectric  constant  at  plus  infinity.  The  parameter 
5  (defined  by  (47)]  is  representative  of  the  rapidity  of 
this  variation,  or  in  the  ease  of  the  symmetric  layer,  the 
width  of  the  layer.  Both  symmetric  and  transition  layers 
are  displayed  in  Fig.  (1). 

The  transformation  (36)  maps  the  entire  real  axis 
a (-«,«)  into  the  negative  real  axis  £(0,-«),  with  3=0 
corresponding  to  5  *  -1. 

The  hypcrgeometric  equation  has  24  solutions. 

In  the  region  |c|<l  (s<0),  a  linearly  independent  pair  of 
solutions  is  given  by*  [Evdilyi  et  at.  ,  1953] 


*That  F2  is  a  solution  of  (35)  may  be  verified  by  direct 
substitution.  By  virtue  of  the  fact  that  the  solutions 
(49)  arc  two  power  series  which  begin  at  different  powers 
of  £,  they  cannot  be  proportional  to  one  another.  F^  and 
F2  arc  thus  linearly  independent. 


r 


3n 


1  Yf(a-Y  +  l,  C  -  Y  ♦  l  ;  2  -  Y  ;  O  » 


(4?) 


where  f(u,  3;  vi  C)  is  the  hypergeometr ic  function  defined 
by 


•”  («  >  3  ,  Y  i  O 


a?  a  (aH  l)3_(3  +  i)  „ 
v  y  Cy+ i  )  T*:“  ^ 


Sl(cr>  1  )_( a  +  2J  zj  l*JJL i+J  )  3 

y(y  +  1  j  (y  +  2)  + 


(50) 


In  order  to  obtain  the  corresponding  solution  in 
the  region  U|>1  (:;>()),  the  hypergeometr ic  function  must 
be  considered  as  tlic  representation  in  the  region  U|<1  of 
an  entire  function.  In  order  to  find  these  solutions  of 
(55)  corresponding  to  those  in  |f|<l,  it  is  necessary  to 
determine  t lie  representation  of  this  entire  function  in  the 
region  |s|>l.  This  procedure  is  illustrated  in  Appendix  A 
and  it  results  in  tlic  following  analytic  continuation 
formula  for  the  hypergeometr ic  function  in  the  region 

U !>i  (=>o) 
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F(atl\o;r.)  =  C  l  (-s)  ~aF  (a  ,  1  -<?+c  ;  1  -  l+a ;  \/z) 


C2  ( ‘**0  ~  JP  >  1  -c+b  ;  I  -s+l ;  1  /  a)  > 


(51) 


where 


Ci  = 


r(7»)r(e-a) 


>  f 


„  ..  r  (e)r(a-- ) 

"  r(«)f(e-io 


(32) 


The  solution  to  the  reduced  wave  equation  with  an  Epstein 
profile  can  thus  be  written  down  as 


!<P 


.....  _  -V,(y*i)/?m  .  (l  +  ct  +  ^-y)/?  . 

, .  U)  ■  k  C  (!■■')  i  let , ,, ;  Y  ;  f  J 


202.(C)  -  rr4>‘(Y'1)/2(l*C)  (1  +  tt+:'Y)/? 


x  F(ci-y+1>0*y  +  1;2*y;C) 


z  >0 


0  1  + 


(O  -  B-*E(Y*I)/J(l-0(l+a*C'1r)/2  I 


•a 


£j('C)  F(a>l‘Y+a»l"3+a»l/c) 


Z0?t(O 


-  »-V‘rl)/!(l-i)“V3’,)/!! 

£3('0V  °  ^"Cc-y+l ;c.'3  +  !  ;1/e) 

+  5.,(‘0Y  1  .-'Cr.-y*l  ,i;  ,  .’.-a*!  ; !/;)  ;■  , 


(33) 


'*>2 


where 


i  .  LCx) r JiL'  aJ 

r(f,  j  r  (y-a) 


_  rc-yjr_(--aj 

3  f  („-Y  +  l)I-(l-a) 


.  r(y)r  (a-s) 
c2  "  r(a)f(Y-B) 


,  -  r  ( 2  -  v )  r  ( a  -  rQ 
r  (a-Y  +  1)  r  (1  -£i) 


(S4) 


The  asymptotic  forms  of  the  solutions  as 
can  be  easily  evaluated  to  be 


D  1 


.  ( :: ) (-1)  e  u 


(55) 


1 


(  ' \  ”  r  “  t  »■  /  2  i  j  ”  - 1  i  «  . 

(a )-►?'!  ( - 1  j  { n j e  1  *  Bo  1  } 


0  2+  )-»«/.  (-a)  («3e  +  C4e 


ik  I ; 


}  . 


Note  that  as  s-»-»  each  solution  is  in  the  form  of  a 
traveling  wave,  one  traveling  to  the  right,  the  other 
traveling  to  the  left.  As  s-»»,  t lie  solutions  arc  made 
up  of  linear  combinations  of  traveling  waves.* 


*Thc  pair  of  solutions  (49)  picked  as  independent  were  done 
so  for  convenience,  as  they  arc  not  t lie  only  independent 
pair.  The  imposition  of  boundary  conditions  would  force 
the  equivalence  of  all  the  arbitrary  pairs  of  independent 
so  1 u t ions . 
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A  comparison  of  (24),  (32)  and  (45)  will  show 
that  for  the  distribution  profiles  of  electron  density 
the  parameters  M  and  IJ  arc  given  by 

M  -  N  «  n|c[(1~ine)/(1+Q2)] ,  (56) 

where  ft.  is  the  maximum  value  of  plasma  frequency  for  the 
V  o 

distribution.  For  both  types  of  distribution  profiles,  the 
collision  frequency  will  be  considered  to  be  constant,  and 
for  convenience  in  considering  both  distributions,  our 
equations  will  only  contain  .7,  knowing  that  U  *  //. 

Also,  the  a-variation  in  mean  properties  will  be 

written 


ena/(l +eWa)*  =  *(*)  (57) 

where  5*1  for  the  transition  layer,-  and  5*2  for  the 
symmetric  case. 

5.2  THE  TURBULF.UT  LAYER 

In  considering  the  propagation  through  a  turbulent 
plasma  layer,  we  assume  that  the  distribution  of  electron 
density  can  be  regarded  as  a  continuous  fluid  and  that 


i fit  fire  smooth  a iu!  slow  (though  random)  functions  of 

space  and  time  [c.g.,  dcl.'olf,  196S).  By  this,  we  mean  that 
the  tine  scale  of  turbulence  [Hinv.c ,  1959)  is  significantly 
longer  tnan  t lie  tine  necessary  for  propagation  across  a 
dimension  characteristic  of  the  layer  thickness.  IVc  may  then 
wrif'-  '  '  ti  p  r  illation  for  a  particular  realization  of  the 

tii  -‘runiK,  1967]  so  formed.  V.’c  also  assume 

that  f  >h'  is  stationary,  i.c.,  the  statistical 

prop'' ~  r  +'  ■'  medium  are  not  functions  of  time  [Rccknann , 
ior.7i  a  second,  spatial,  ensemble  is  formed  by  the  electron 
density  as  a  function  oi  location  at  a  particular  time. 

‘ '  ‘  .  changes  of  electron  density  arc  governed 

bi  bid  dynai .  i .  .  .  ui  ni»  i  once ,  v;c  assume  that  the  time  and 

snnrr  nseiibU  *o?h  stationary  and  equivalent  [deWolf, 

I- 


The  distribution  of  electron  density  n,  will  be 

•/\ts  as  follows 


if  =  if  ^  ( ~ )  +  n  j  ( a  ,  t ) 


(58) 


vbrr  r» 
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The  brackets  <  >,  represent  an  ensemble  averaging  process 
[Bcakr.ann,  1967].  In  describing  t he  turbulent  medium,  we 
shall  have  use  for  the  following  parameters:  the  turbulent 
intensity  c ,  where 


C  =  <n*>/n0,  (00) 

m 

and  the  normalized  electron  density  fluctuation  n,  where 

t . 

n  e  u ;  .  (61) 

The  fluctuating  component  of  electron  density  is  then 
written 


ni  =  «0(s)cn(.*!)  •  (62) 

We  shall  restrict  this  analysis  to  homogeneous 
turbulence  [//tune,  1958]  by  taking  c  to  be  constant 
throughout  the  plasma  region. 

For  notational  convenience  in  this  section,  we 
will  write  the  squared  wave  number  k2 ,  as  being  composed  of 
a  mean  non-fluctuating  component  and  a  random  component,  so 


w 
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that  we  write 


;-2  = 


] 


(03) 


The  exact  terra  correspondence  of  this  expression  may  be 
obtained  by  a  comparison  of  (33),  (34),  (45),  (56)  and  (58). 
The  wave  equation  (41)  may  thus  be  written  in  the  form 


'/dr. 


o  k  ]  1  C  =  0  . 


(04) 


We  take  the  term  to  be  a  source  term,  so  that  the 

technique  of  Green’s  function  may  be  applied,  i.c.,  after 
writ  ing 


id2 /ua2  ♦ 


'■  n K  }  4 


(05) 


tire  left  side  of  tire  equation  is  just  that  of  an  Epstein 
profile. 

Let  us  now  construct  a  Green's  function  d(a,a') 
for  (65)  which  is  continuous  and  corresponds  to  an  impulsive 
force  at  a  =  a'.  We  subject  the  Green's  function  to  the 
condition  of  outgoing  waves  as  |s  )->•«.  Thus  G  satisfies 
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the  equation 


[d2/dz2  +  k  (a)  }  G  (a,  a'  )  =  5  ( -  a  1  )  , 


(C(.) 


with 


S (-»,-') 

G  ( 00  >  a  *  )  «  c 


-ik, 


(07) 


By  integrating  (06)  with  respect  to  a,  wc  may  determine  that 
6(3,3' )|  is  discontinuous  by  the  value  1.  The  properties 

of  the  Green's  function  can  then  be  summarised  as 


1  . 

2. 

We  are  seeking  the  solution  to  the  problem 


C,  is  continuous  at  •;  -  a':  0v\  -  Gn  •  =  0 

3+  c-  (OS ) 

G'  is  discontinuous  at  a  =  a '  :  6 '  ,  -  <J\  =  1 . 

**  ^  3 


{d2/dz 2  +  k§ic(o)}A  =  /(a),  (09) 

with  an  incident  and  reflected  wave  as  and  a  trans¬ 

mitted  wave  as  s-*-®.  These  boundary  conditions  are 
expressed  as 


where  k  j  nay  cither  be  equal  to  (symmetric  layer)  or 
different,  as  in  the  transition  layer.  Let  ancl  Z2  be  a 
pair  of  independent  solutions  of  the  homogeneous  equation 


( dT/d:: ?  *  ,';s<  =  0, 


(71) 


with  the  following  asymptotic  forms* 


z2m  r.  n: 


iC 


2?(-»0  =  Pe'lKoi 


*  i,c 


■IK 


}. 


(72) 


The  form  of  the  Green's  function  can  then  be  immediately 
written  down  as 


ICS 


0(r.,r.') 


z>: 


J  (*» »"  ) 


<2  1  ^  j  (  ~  )  +  ^  2  “  2  ^  n  ^ 

^  ~  j  (-*  )  + 


(73) 


* The  sc  are  the  asymptotic 
reduced  wave  equation  as 


forms 

given 


of 


the  solutions  to  the 

(55). 


c\  n 
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By  applying  the  outgoing  wave  property  of  G, 


Gi (’-.s') 
G2(»,  s') 


GjZ,  (s)  +  c2Z2  (3)«e 


(")  +  "2  (")a«^M“ 

e1',:i‘'[i1/.£1  ♦  b2DC2] 

*  e  vkl"[by'.r2  +  h2 :rj  , 


and  then 


a 


1 


/«V 


(74) 


(73) 


The  continuity  properties  of  C,  (6S),  lead  to  the  following 
set  of  equations 


MV-"')  ■  f£f  V"’))  •  a2-2 (“ ■  3  *  0 


mzi 1 1  -  fit  si (*■>]  •  «2=;c*’)  *  1. 


(76) 


I 

These  yield 


-2a(a')/r/; 


[-ZjCa')  +  US2/^4)Z2(3-)]/tV 


(77) 


where  .7  ib  the  l.'ronsk i«.n  determinant  [3-.iir.viHc,  1964  , 


Chap.  7j 


-  *  l-  .(•‘O-IOO  -  J (;0  1 _ 

i  *.  1  *-  <•  "  M 


I  # 


(78) 


Since  the  l.’ronsk  i  an  of  the  wave  equation  is  constant*,  it 
may  be  evaluated  at  any  convenient  point.  The  Green's 
function  is  then  given  by 


- 1 


(^V)  [:•  1  C-:  *  )”2 


(a^V.-OVM 


(79) 


?  (•:  '  )-j  (")/,t 


°  2  ( "•  )  ^  2  ^  ‘‘  ^  2  ^  L 'i*  4  ^  * 


To  make  use  of  the  Green's  function,  we  multiply  (69)  by  G, 
and  (66)  by  Tk  to  obtain 


+  k ^ k  (:: )  J v 


(80) 


Subtracting  the  tv.o  equations  and  integrating  over  z  from 
- «  to  the  left  hand  side  may  be  rewritten  as 


1 9  S  6  ] . 


*This  is  a  special  case  of  Abel's  identity 
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no 


while  the  integral  containing  the  delta  function  may  be 
immediately  evaluated.  Therefore 


•  o> 


By  evaluating  the  bracketed  term,  using  the  asymptotic  forms 
of  G  and  $,  and  interchanging  the  variables  a  and  s',  then 


4  (") 


>222 


(S3) 


where 

z<z  ’ 

C(z'  ,z)  «  G^\z)  «  (.2ik0B)-1 

*  {ZjOOZ^a’)/*  -  F2Z2(;0Z2C2')/FB4} 

(84) 

Z>Z  '  _  J 

G(z' tz)  «  ffj(«',»)  *  (2£fe0S) 

X  {Z1(2)Z2(2')/4  -  Z2(2)Z2(s')F2/BS4}. 
Comparison  of  (41)  and  (64),  allows  us  to  write 
down  the  constants  A  and  3,  and  the  exact  formal  solution 
for  the  case  of  the  turbulent,  media.*  Note  that  the 


*For  our  situation,  (83)  does  not  represent  a  solution,  but 
represents  an  integral  equation  in  the  desired  quantity. 
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proper  representation  of  the  homogeneous  solutions  Zj  and 
Z2  depends  upon  whether  their  argument  is  greater  or  less 
than  0.  The  formal  "solution"  is  thus  an  integral  equation 
for  the  field  $ .  To  evaluate  the  reflection  or  transmission 
coefficient  of  the  layer,  the  field  is  evaluated  at  -®  or 
®  respectively. 

To  solve  (83)  we  shall  use  the  Neumann  iteration 
technique  [Green,  1969]  which  generates  a  series  for  the 
solution  of  the  reduced  wave  equation  at  -•»  of  the  form 


<*(-»)  *  eikoz  +  [j?0  +  +  r2  +  ... ]e~iko z.  (85) 


The  terms  within  the  parenthesis  represent  the  successive 
contributions  to  the  field  reflection  coefficients  given 
by  succeeding  iterations  of  the  integral  equation  (83) 
(These  are  physically  equivalent  to  multiple  reflections.). 
If  the  series  representation  of  the  reflection  coefficient 
is  convergent,  we  may  attach  to  each  term  a  measure  of 
turbulent  strength  e,  so  that  we  may  write  the  field 
reflection  coefficient  as 


D 

“v 


c°Z 


(36) 
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The  most  convenient  measurable  parameter  of 
interest  is  the  average  power  reflection  coefficient 
defined  by 

R  =  <RVR*>.  (87) 


Using  (86)  and  retaining  terms  to  order  e2  in  smallness,  we 
have 


r  *  <  |«0  |2>+<i?g/7l>+<i?;i?0>+<  |/?1  l2>+<.*?;*2>+</?0i?*> . 


(38) 


Each  higher  order  of  the  reflection  coefficient  R  is 
proportional  to  an  integral  over  a  higher  power  of  the 
random  quantity  n,  i.e., 


R o  “  constant 


*1 


/■lOOnCa’Ms’ 


(89) 


R2 


f  2(2')n(a,)d3 


/ 3  (“ ')  n  (s”)cV 


The  medium  properties  were  split  into  mean  and  fluctuating 
components  by  (58)  so  that  <n>  =  0.  Thus  <i?3>  =  RQ  and 
</?!>  =  0.  The  reflection  coefficient  to  second  order  is 


R 


2T.c  t 


(90) 


then 


* 
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By  convention,  the  reflection  coefficient  is  broken  up 
into  coherent  and  incoherent  components,  where  the  coherent 
reflection  coefficient  is  defined  as 

-  l<vl2.  t91) 


while  the  incoherent  coefficient  is  given  by 


WT  .  <  \R  .  <R  >  | 2> 
I  1  y  v  1 


(92) 


To  second  order  in  e, 

wc  •  \h0\z  *  :r.e{/^<-2>} 

Wj  =  <|flj2>. 


(93) 


We  shall  also  make  use  of  the  special  reflection 
coefficient  <i?'>  used  by  Hochstim  and  Martens  [1969;  1970] 
defined  as 


<R'>  *  R 


(94) 
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To  second  order  then, 


<fl’>  «  <|i?l|2>  +  2Re{/?*</?2> } .  (95) 


Two  iterations  of  the  integral  equation  (83) 

(which  corresponds  to  .*?0  +  Rl  +  R}  in  the  Neumann  series) 
then  yield 

♦  («)  ■  »*C-irlS/!([z01W-C-i)t5z0JW/V 

oo 

+  2 klNZ  |  ^(2,)n(2')[Z01(2')-(-l)i‘9fl2Z02(2‘)/B4]C(a'  ,s)dz' 
-  00 

(96) 

00 

+  (feo^t)2  j  ^(a,)n(s,)5(s'  ,z)dz' 

-  00 

00 

x  |  ^(2')n(a,,)[Z01(2")-(-l)i5F2Zo2(sM)/B4]ff(a",s')d3"} 

-  00 

where  ff(2',s)  is  given  by  (84). 

The  incoherent  power  reflection  coefficient  is 
then  given  by 


Wx  «  (mfe0/2)2  |//|  2;2{ 


l  g(z')Z*\U')dz'  <n(a,)n(3")>^(a,’)Zc?(?.")i2" 


+ 


♦  to 


4|02/B4|  J  GCe'JZoiCa'JZoiCa' 


<n(3')nC^,)>^(2”)z0l(3")z02(2")d2" 


|£2/Bj2  | 


9(z')Z02(3')dz' 


*  \  <11  (2  1  )n (2")>0 (2")Zo2 (2")df2" 


2Re  [  -  2  (  - 1)  (P2/B4)  j  ffCa'JZ^Ca')^,' 


(97) 


*  j  <n(2')n(2")>^(2M)Z01(2")z 


C-l)l25(B2/54)2  j  ^  (2 '  )zjx  (2  » )cfs  • 


x  J  <n  (2  '  )n  (2M)>^  (2")Zo2  (2")<f2" 


2(-l)  \B2/B^\  (32/54)  ^  &”(2,)ZcJ(5,)Z0J(3')cf2' 


j  <n(2 '  )n  (2")><7  (2M)z0f  (2")j2"] } . 


t 
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Thu  Lviri''  «.  uni  >•  <n>  have  vanished.  Similarly,  the 

average  inmlir-.-nt  transmission  coefficient  is  given  by 


’  v-V  /.  i5  2{ 


I 


<r)Cz'  )r)(zn)>g’(s")Z0l  (zn)Z02(z,,)dzn 


(98) 


<n(a')n  (z")>ff(z")Z02(z")dz" 


2?’  V->  )■  J  ^  (3  ’  )  2  o  t  (z')Z*02(z')dz' 


<n  (£,)r,(s,,)>?(2”)Z0‘(2")d2"]>. 


•  00 
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A  typical  tern  of  (98)  or  (99)  has  the  form 


(99) 


j  <n(x')n(x)>s(x)Zoi(x)Zo2(x)dx. 


Expressions  involving  similar  terms  arc  obtained  for  the 
coherent  fields  through  the  use  of  (93)  and  (96). 

The  correlation  <n(aOn(tf*  )>  will  be  assumed  to  be 
of  the  exponential  form* 


<  n  (~ ' )  n  (x )  > 


(ICO) 


Because  of  the  varying  format  of  the  functions 
Z o i  and  Z02 ,  the  integration  of  a  term  like  (99)  must  be 
broken  into  six  steps  as  follows 


*Such  a  correlation  function  is  mn.lytically  convenient  and 
yields  a  turbulence  spectrum  which  is  consistent  with 
experimental  results  [Guthart  o~  at.,  1966], 
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0 

|  g(x')Z<n-(.x')ZQ2-(x')e  bx  dx ' 


X 


a?' 

f  bx 

e  g  (x)Z0 1.  (x)Zq2- (.x)dx 


\  g C«’)s 


0l-C*#)*0a-C*#)«  dx ' 


|  e  iX?(^)Z01.  Car) 2 o 2 -  (x)<ix 


^(^,)z0i-C«,)z02-^x,5ciX  dx' 


-ix 


x  I  e  ^U)Z01  +  C«)Z02+(ar)dar 


^(=')Z01+(x,)Zo2+C:c')e  iX  dx ' 
0 


|  eix^(x)Z0l_ (x)ZQ2_ (x)ix 


(101) 


♦ 


5ft 


j  f* ' (*' )Zg2t (x* )e"^X 'dx* 


[  bx 

*  J  e  0(*)Zol+(x)ZO2+(x)dx 
0 


9  (*  )  Z  0  1  ♦  ^x  '  )  Z  0  2  ^x  '  )  ®  ^x  ' 


f  -ix 

*  J  e  g(x)ZQl+(x)ZQ2+(x)dx. 
x  * 


Wc  carry  out  the  integration  of  such  a  typical  term  in 
Appendix  B.  As  a  result  of  the  algebraic  complexity  and 
length  of  that  integration,  wc  defer  the  final  expressions 
to  Appendix  C,  where  the  results  are  given  in  terms  of 
certain  functions  defined  in  Appendix  B.  These  functions 
are  in  general,  sets  of  infinite  scries  which  must  be 
computer  evaluated. 

The  relative  simplicity  found  in  the  reflection 
and  transmission  coefficients  of  the  undisturbed  plasma 
U  *  0)  does  not  carry  over  to  the  turbulent  case,  even  for 
the  restrictive  situation  of  small  perturbations  from  the 


smooth  distribution. 


To  compute  the  turbulence  free 


reflection  coefficient  it  is  necessary  to  know  the  field 
only  at  where  the  asymptotic  forms  are  simple  exponen¬ 
tials.  In  the  disturbed  case  however,  the  zeroth  order 
field  must  actually  be  known  at  every  point  in  space  in 
order  to  sum  up  the  scattering  introduced  by  the  random 
inhomogeneities  throughout  this  space.  The  exception  to 
this  complication  is  the  case  for  S->0  in  the  transition 
layer.  This  distribution  then  becomes  that  of  free  space 
on  the  left  of  zero  and  a  half  space  of  plasma  of  constant, 
mean  electron  density  on  the  right  of  zero.  The  results 
of  this  case  are  of  sufficient  interest  to  promote  a 
separate  discussion  in  Section  6.1.  The  actual  results 
of  the  computer  evaluation  of  the  final  expressions  are 
discussed  in  Section  6.2. 
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6.  DISCUSSION  OF  RESULTS 

6.1  THE  TURBULENT  HALF-SPACE 

As  discussed  in  Section  5.2,  if  we  examine  the 

M 

ease  of  the  transition  layer  when  we  allow  S+ 0,  we  arrive 
at  the  problen  of  a  half-space  distribution.  The 
half-space  distribution  lias  already  been  studied  separately 
by  several  workers  [Rosenbaum,  1969;  Kupiea  et  at .  ,  1969], 
who  have  made  use  of  smoothing  equation  techniques.  It  is 
not,  however,  our  goal  to  compare  approximate  methods  of 
solution,  but  to  compare  our  solution  method  with  the 
exact  experimental  treatment  of  the  slab  geometry  [Hochstim 
and  Martens,  1969;  1970]  (where  applicable),  and  to  show 
that,  in  fact,  the  solution  restrictions  which  we  derive 
in  this  section,  convey  proper  information  about  the 
validity  of  our  method  of  solution. 

6.1.1  Comparison  with  the  Hochstim  and  Martens  Computer 
Experiment 


Since  the  only  experimental  solutions  available 
are  for  a  slab  geometry,  comparison  with  a  half -space  may 
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be  made  only  for  a  restricted  number  of  experiments. 

These  are,  those  eases  in  which  the  attenuation  of  a  wave 

across  the  width  of  the  slab  is  so  great,  that  there  is 

effectively  no  contribution  to  the  reflection  coefficient 

made  by  the  second  jump  in  medium  properties  (i.c.,  the 

rear  of  the  slab).  We  have  arbitrarily  chosen  the  criteria 

that  an  experimental  slab  can  be  compared  to  t lie  half- space 

distribution  (all  other  medium  properties  being  the  same) , 

if  a  wave  suffers  10  dB  of  attenuation  in  traveling  across 

the  width  of  the  slab.  This  restricts  the  comparison  to 

the  experimental  cases  which  have  1.5>  n  >  0.8,  n  >  0.01. 

P  c 

As  5-*-0,  the  constants  occurring  in  the  Epstein 
profile  solution  B\-B 4  may  be  evaluated  by  using  the 
following  expansion  for  the  gamma  function  when  2  is  small 
[Abromowitz  and  Stegun,  1964] 

r(a)  =  1/s.  (102) 

These  constants  then  become 


*  B4  -  (1+X)/2X 

B2  =  V 3  •>  U-l)/2*, 


(105) 
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I, 

where  X  =  +  The  honogcncous  solution  for  the 

Lpstein  t  ransition  profile  can  then  be  easily  evaluated  to 
be 


rt 

U 


(-D 

(-D 


iS/*eikor. 
-IS/ 2  e~ikoz 


(X-l)e~vklz]/2X 
(X-l)e~iklS]/2X , 


(104) 


The  exact  integral  equation  is  then 


a<0  • 

♦  (*)  -  c*=ottlptl‘fne)/ j  «)<*«' 

0 

♦  ♦  [il-X)/(l*X)]e~ik°Z 

(105) 

3>0  • 

♦  (a)  *  cfe^(n2(i-?*n^)/(i+n2) J  J  n (*')  +  (*')<% (a • 

o 

*  [2/Cl**)!*^0**, 
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with  the  Green's  function  given  by 


Z  '<3 
a  '<0 

z '  >z 
a'<  0 

z '  >a 
a '  >0 


a'<a 
a  '<0 

a  '<a 
a.'  >0 

s  *  >a 
a’>0 


G_(z'  ,z)  - 

G  (a  '  ,  a)  - 
•» 

G_ (a'  ,a)  - 

G+(a'  »a)  ■ 
G+(a'  ,a)  « 
ff+C3'  ,a)  « 


(Uk„yl[eik  .2 

C2<fe0)-1  [eifco2' 


(2£&o) 

(2tfc0) 

(2«0) 


ik^Xz  +  ^-a^aj^^a' 

tfcoXz'  + 


The  expression  is  just  that  obtained  by  considering  only 
the  half-space  geometry  [c.f.,  Rosenbaum ,  1969].  By 
applying  the  asymptotic  forms  as  | s  | -►<»  to  the  expressions 
for  the  scattering  coefficients,  we  may  straight-forwardly 
obtain  the  results  for  the  half-space  which  arc  displayed 
in  Appendix  C. 

6.1.2  Validity  of  the  Solution  Method 

Following  the  method  of  analysis  given  by  Prise': 
[196S]  for  the  validity  of  the  Born  series,  we  may  derive 
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a  similar  validity  for  the  half- space  geometry  with  exact 
Green's  function.  Denoting  the  successive  terms  of  the 
Neumann  series  by  it  follows  from  the  integral  form  of 
the  Minkowski  inequality  that 

V 

,  (107) 

n2(l -inj  f 

-  IH2J  \G(z'  ,z)\*dz' . 

V 

Expression  (107)  may  also  be  obtained  from  a 
convergence  criteria  for  the  Neumann  series  derived  by 
Jarem  [1969].  The  expression  (107)  implies  a  finite  volume 
of  plasma,  whereas  a  half- space  of  plasma  extends  to 
infinity.  If  we,  however,  allow  collisional  damping,  any- 
propagating  wave  is  continuously  attenuated  and  an  effective 
volume  is  created.  The  final  form  of  the  convergence 
criteria  will  show  that  when  there  arc  no  collisions,  our 
method  of  solution  does  not  yield  even  an  approximate 
solution. 

From  (106) , 


n2(l-in  ) 
*  2  <  k2  JLS-  .  Q* 

n* l  ’  o5  (i+n2) 


G{z\z) 


e-ik^z  • 9ik0Xz 


(108) 
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and 


The  integral  in  (107)  is  then  given  by 


/  -  ll+jrl)"1 # 


so  that 


<  c 


np 


- 1 


8 


■|l+*l  U„l 


With  a  complete  function  space,  the  convergence  of  £  $ 

•  »“1 
follows  from  the  convergence  of  £  jn2.  Therefore,  if 

n*l 


(109) 


(110) 


dll) 


n 


c#*/Ui«opV|i«xll  <  l, 


(112) 
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the  rapidly  convergent  Neumann  scries  may  be  used  as  an 
approximate  solution  to  the  reduced  wave  equation. 

Expressing  the  bound  as  a  limit  upon  the  allowable  turbu¬ 
lence  level,  we  have  the  restriction 

C  <  S’(l+fi2)*[(l*a')2  ♦  (e’)2]W.  (113) 

c  p 

• 

The  limiting  cases  of  inequality  (113)  arc 
displayed  in  Fig.  (2).  It  can  be  seen  there,  that  the 
restriction  (113)  is  consistent  with  the  physical 
interpretation  of  the  scattering  integral  equations.  Since 
our  iteration  method  considers  only  a  finite  number  of 
scattering  events  (one  event  for  the  incoherent  coefficient; 
two  events  for  the  coherent  coefficients),  our  solution 
will  not  be  valid  when  the  medium  is  such  as  to  be  very 
multiple  scattering.  Such  a  medium  would  be  one  which  is 
not  very  absorbing  and  of  large  physical  extent,  so  that  a 
wave  would  not  diminish  greatly  in  magnitude  following  a 
scattering  event,  and,  in  fact,  could  be  scattered  again 
wit  out  appreciable  attenuation.  Obviously,  a  theory  which 
allows  for  at  most  two  scattering  events  cannot  be  expected 
to  faithfully  predict  the  scattering  coefficients  for  such 
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a  media.  In  a  plasma,  the  attenuation  coefficient  goes  to 

zero  when  fi  <1  and  0  ->0.  As  can  be  seen  in  Fig.  (2),  for 

p  o 

this  situation,  the  turbulence  level  must  satisfy 

c  <  (nc/2)*.  (114) 

Or  for  a  non-attenuating  medium, 

C  <  0  ,  (115) 

i.e.,  the  Neumann  series  cannot  give  a  solution  for  any 
level  of  turbulence,  no  matter  how  small. 

In  Figs.  (3) -(7)  we  have  plotted  the  results  of 
the  reflection  coefficients  for  incoherent  reflection,  and 
the  special  reflection  coefficient  <R'>  for  the  experimental 
conditions  of  the  Hochstim-Martcns  experiment. *  Superimposed 
are  their  experimental  results  for  parameter  values  where 
the  slab  and  half  space  may  be  compared  according  to  the 
attenuation  diteria  previously  discussed  (Section  6.1.1). 

It  can  be  seen  that  for  the  largest  value  of  collision 


*The  only  parameter  not  given  explicitly  in  these  figures 
is  the  correlation  length.  In  the  computer  experiment  this 
was  taken  to  be  variable  according  to  the  form  kaar./2?.  . 
This  same  form  has  been  used  in  our  calculations.  ^ 
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frequency  (nc=1.0),  the  tv.o  solutions  arc  in  exact  agreement 
even  for  large  values  of  turbulence  level,  c.  As  the  plasma 
becomes  less  attenuating  for  fic<l  (by  decreasing  the 
collision  frequency),  our  result  begins  to  differ  from  *hc 
computer  experiment  at  the  higher  levels  of  turbulence*, 
and  in  fact  for  ftc*0.01  physically  unrealizable  incoherent 
reflection  coefficient  arc  obtained  (i.c.,  greater  than  one). 
Such  a  divergence  can  be  allowable  only  if  the  plasma 
parameters  are  such,  as  to  render  the  inequality  (113)  not 
valid. 


As  an  example,  i/c  have  examined  the  convergence 
situation  for  0^*1. 0,  the  value  of  plasma  frequency  for 
which  there  is  the  greatest  disparity  between  solution  and 
experiment.  In  Fig.  (S) ,  we  have  plotted  the  ratio  of  the 
perturbation  solution  incoherent  reflection  coefficients  to 
those  of  the  computer  experiment,  as  a  function  of  turbu¬ 
lence  level,  c  and  collision  frequency,  In  addition, 

we  have  plotted  the  bound  given  by  (113).  It  can  be  seen 
that,  if  we  do,  in  fact,  restrict  the  region  of  validity  to 
that  given  by  (113),  then  the  reflection  coefficient  will 


*This  variance  is  because  of  our  solution  method,  not  because 
the  half  space  and  slab  cannot  be  compared.  The  10  dB 
attenuation  criteria  is  still  met. 
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differ  by  no  more  than  1.5  dB  from  the  exact  solution. 

Speculation  has  arisen  over  the  reason  why  the 
reflection  coefficient  <tf’>  takes  on  negative  values. 

Boohstim  and  Martens  [1969]  report  that  this  is  due  to 
negative  fluctuations  which  near  decrease  the  mean 

reflection  coefficient  more  than  positive  fluctuations  are 
able  to  increase  it.  Additional  speculation  has  attributed 
these  negative  values  to  the  integral  wavelength  slab 
thickness  chosen  by  Hochstim  and  Martens.  Nc  do  however 
sco  in  Figs.  (6)  and  (7),  that  the  same  phenomenon  is 
observed  for  the  half-space,  negating  this  last  presumption. 

It  is  our  opinion  that  this  negative  value  is  due 
to  the  average  additional  path  length  traveled  by  a  wave 
which  has  undergone  two  random  events,  compared  to  the  wave 
which  has  been  scattered  by  the  mean  distribution.  By 
examining  (90),  we  sec  that  since  a  double- scattering 

term  2  is  the  lowest  order  term  which  can  yield 

a  negative  value.  Since  this  term  is  of  the  same  order  as 
the  incoherent  term,  <|/?jJ2>,  it  is  not  surprising  that 
their  sum,  <A’>  takes  on  negative  values  for  some  region  of 
parameter  space. 

The  following  necessary  condition  for  the  validity 
of  a  perturbation  technique  for  the  slab  geometry  has  also 
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been  derived  by  Jnreiu 


and  'lavicnz  , 


1970] 


c2  <<  [ (l-n2)2  +  n2]/n4 .  (116) 

P  0  P 

Such  a  Kccocc.zry  condition  is  of  limited  use,  however,  since 
expression  (116)  can  only  tell  us  that  the  technique  night 
be  g<  od .  As  an  example  of  the  confusion  such  an  expression 
raises,  let  us  examine  (116)  for  '.'C,<<1  and  rip < ^  1 . 

Inequality  (116)  can  be  seen  not  to  restrict  the  volume  size 
in  any  manner  as  the  plasma  becomes  nonattenuating ,  a  ease 
for  which  the  perturbation  solution  obviously  breaks  down 
with  a  large  enough  volume. 

6.1.3  Calculations 

In  Figs.  (9) -(11),  we  have  mapped  the  incoherent 

reflection  coefficient  as  a  function  of  plasma  frequency, 

with  collision  frequency  and  correlation  length  as 

parameters.  The  plots  arc  all  similar  in  form,  although 

there  is  a  shift  in  the  value  of  plasma  frequency  for 

maximum  reflection  (This  value  decreases  to  n  =1  as  n  -*-0.). 

P  a 

The  shift  may  be  related  to  the  onset  of  surface 
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scattering*  with  the  increasing  opacity  of  the  plasma, 
as  the  plasma  frequency  increases.  A  psuedo- surface  is 
formed  which  transmits  a  decreasing  amount  of  illumination 
as  the  plasma  frequency  is  increased,  because  of  absorption 
and  the  shielding  of  the  inner  plasma  from  illumination 
(and  subsequent  scattering) . 

For  values  of  plasma  frequency  greater  than  that 
of  maximum  incoherent  reflection,  the  incoherent  reflection 
coefficient  is  relatively  insensitive  to  changes  in 
correlation  length.  Where  this  docs  occur,  the  shin  depth, 
of  the  plasma  is  shorter  than  the  correlation  length,  and 
the  fluctuations  in  electron  density  arc  relatively 
uncorrclatcd  (because  of  the  large  absorption  over  the 
distance  of  a  correlation  length).  Wo  thus  expect  a 
saturation  with  increasing  correlation  lengths.  As  can  be 
seen  from  Figs.  (9)- (11),  at  a  given  value  of  plasma 
frequency,  departure  from  the  saturation  values  of 
incoherent  reflection  coefficient  first  occurs  for  the 


*Surfacc  scattering  is  scattering  primarily  from  a  thin 
layer  of  plasma,  with  very  little  transmission  through 
the  layer. 

tSkin  depth  is  the  distance  to  which  a  wave  will  travel 
before  its  amplitude  decreases  to  1/e  of  its  original 
magnitude  [ Corson  and  Lovvain,  1962], 
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lowest  value  of  correlation  length,  which  is  consistent  with 
t he  above  explanation. 

6.2  The  Turbulent  Transition  and  Symmetric  Layers 

In  Figs.  (12)  and  (13),  we  give  the  incoherent 
transition  layer  reflection  coefficient  for  two  values  of 
collision  frequency  (n0  -  1.0,  0.1)  with  the  layer  thickness 
S  as  a  parameter.  The  correlation  length  is  taken  to  be 
ka  *=  1.0.  From  Fig.  (13)  for  ■  0.1,  it  can  be  seen  that 
the  effect  of  increasing  the  transition  layer  thickness  is 
to  decrease  the  reflection  coefficient  in  the  underdense 
regime,  while  shifting  the  plasma  frequency  of  maximum 
reflection  to  higher  values.  Fig.  (12)  also  displays  the 
decrease  in  underdense  reflection  with  increasing  layer 
thickness,  and  from  a  comparison  with  Fig.  (10),  it  can  be 
assumed  that  at  higher  values  of  plasma  frequency,  the 
reflection  coefficient  also  peaks  and  drops  off  for  a 
collision  frequency  of  *1  *  1.0.  The  decrease  in  incoherent 

reflection  with  increasing  layer  thickness  can  be  attributed 
to  a  decrease  in  the  overall  electron  density  distribution 
as  the  layer  thickness  is  increased  [c.f.,  Fig.  (1)].  The 
shift  of  peak  reflection  may  be  due  to  the  incident  wave 
propagating  farther  into  the  plasma  (because  of  the  reduced 
level  of  electron  density  at  a  given  point),  and  effectively 
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seeing  a  larger  volume  of  plasma  before  encountering  surface 
scattering. 

In  examining  the  symmetric  layer  in  Figs.  (14) 
and  (15),  we  see  that  the  incoherent  reflection  coefficient 
increases  with  increasing  layer  thickness,  both  for  the 
underdense  and  overdense  plasma.  This  variation  is 
consistent  with  the  observation  that  the  reflection 
coefficient  should  approach  zero  as  S-*-0,  for  then  the  plasma 
distribution  occurs  at  only  one  spatial  point. 

Although  we  have  presented  results  only  for  one 
value  of  correlation  length,  wo  feel  that  Figs.  (9) -(11)  for 
the  halfspace,  may  be  used  in  conjunction  with  Figs.  (12)- 
(15)  to  obtain  a  feel  for  the  results  in  the  more  general 
case. 

6.2.1  Validity  of  the  Solution  Method 

In  Section  6.1.2,  v.rc  presented  a  validity 
condition  for  the  plasma  half-space.  Since  that  case 
represents  the  largest  inhomogencity  in  electron  density 
(the  finite  step),  it  seems  reasonable  that  the  restriction 
(113)  should  insure  a  valid  transition  layer  solution  when 
S>0,  although  it  nay  then  represent  a  conservative  estimate. 
For  the  symmetric  layer,  ^ 1 1 3 )  represents  an  even  more 
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conservative 
dealing  with 


bound,  for  in  the  symmetric  ease,  we  arc 
a  volume  which  is  almost  finite. 
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7.  CONCLUSIONS  AND  RECOMMENDATIONS 

In  this  research  paper,  we  have  attempted  to 
advance  the  field  of  turbulent  scattering  into  the  area  of 
generalized  plasma  distributions.  The  price  paid  was  the 
necessity  of  dealing  with  a  Green’s  function  which  was  not 
of  the  convolution  form,  and  because  of  the  complexity  of 
this  function,  we  restricted  our  analysis  to  the  most 
obvious  perturbation  solution  method.  Iv'c  have  shown  that 
even  such  a  simple  method  applied  to  the  problem  of  a 
randomly  inhomogeneous  medium  yields  results  which  must  be 
computer  evaluated  for  examination.  Over  the  parameter 
ranges  which  we  considered,  we  obtained  no  results  which 
could  not  be  plausibly  explained  on  a  physical  basis,  and 
where  comparison  with  experiment  was  possible,  it  was  much 
better  than  one  might  casually  anticipate  in  applying  a 
perturbation  technique. 

Because  of  the  complexity  of  the  lowest  order 
perturbation  solution,  we  do  not  recommend  further  investi¬ 
gation  using  such  a  procedure.  Although  additional  computer 
calculations  might  be  made,  we  believe  that  the  section  of 
parameter  space  for  which  we  have  presented  calculations 
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includes  those  regions  where  unusual  effects  Might  be 
observed,  and  for  which  our  solution  is  valid. 

We  do  however,  believe  that  the  Epstein  profile 
offers  sufficient  advantages  in  studying  inhomogeneous  and 
random  media  to  demand  further  research  effort.  The 
relatively  small  parameter  range  allowed  by  our  bounding 
inequality  docs  not  permit  a  valid  solution  for  the  useful 
ease  of  a  high  altitude  re-entry  wake  (in  which  (1^=0).  The 
smoothing  equation  discussed  in  Section  2.3  has  been 
successfully  applied  to  a  number  of  collisicnlcss  plasmas, 
while  a  modification  of  the  smoothing  technique  has  shown 
exceptional  promise  in  comparison  with  experimental  data 
I Jarem,  1970].  Although  these  techniques  arc  rather 
involved  for  even  simple  Green's  functions  (c.g.,  the  plasma 
slab),  we  feel  that  a  firm  understanding  of  randomly 
inhomogeneous  scattering  processes  dictates  that  one  of 
these  methods  be  applied  to  the  Epstein  profiles.  As  a 
possible  start,  investigation  could  be  made  into  determining 
a  transformation  for  the  Epstein  profile  equations  which 
would  have  benefits  equivalent  to  those  of  the  convolution 
forms  in  Fourier  transformation. 

Numerical  techniques  which  have  been  developed  to 
solve  the  lc.  cst  order  smoothing  equation  \Ua3cal ,  1971], 
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could  be  applied  for  special  cases  to  the  Epstein  profile. 
Such  a  technique  would  directly  make  use  of  the  Green's 
function  obtained  in  Section  S.2. 

An  additional  topic  which  could  be  usefully  studied 
with  the  assistance  of  the  Epstein  profiles  is  surface 
scattering  and  its  associated  effects.  Such  an  investigation 
would  be  valuable  even  without  including  the  effects  of 
turbulence.  The  choice  of  a  proper  effective  surface  for 
scattering  formulations  based  on  surface  scattering  [e.g., 
Jarem,  1964;  Waenecki,  1968]  has  been  carried  out  rather 
arbitrarily  (One  method  employs  a  skin  depth  criteria 
[Permanent  et  al.  ,  1967].),  and  we  thus  recommend  that  a 
more  comprehensive  investigation  of  nonturbulent  surface 
scattering  be  carried  out  using  the  Epstein  profiles.  A 
subsequent  study  night  investigate  turbulent  effects. 

Finally,  we  conclude  that  the  primary  value  of 
this  research  paper  lies  in  its  action  as  a  probation  to 
determine  the  difficulty  of  the  considered  problem.  We  arc 
hopeful  that  these  results  may  give  one  an  insight  into 
the  next  appropriate  approximation  to  apply  (e.g.,  by 
considering  the  asymptotic  functional  dependencies),  which 
might  eventually  produce  n  more  general  and  efficient 


solution. 
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HOT  REPKOOUCIBlt 

APPHNDIX  A 

REPUESHNTAT ION  OF  Till:  HYPIiKGKOMHT I 
Let  us  examine  the  integral  [;•* 

1953] 


^  rC 

where  the  path  of  Integration  is  chosen 
the  integration  path  extends  over  the  in 
the  poles  of  the  functions  r(c*i)  and  r( 
of  the  path  of  integration  and  the  pole: 
r(-t)  lie  to  the  right  of  the  path.  To 
remainder  of  the  integration  path  Jet  uf 
integrand  for  1 1 1 -►•*> . 

F~on  Abroraovitz  and  Steguu  [!'. 
form  of  the  gamma  function  is  given  >y 

r(t+l) 
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Using  the  relation  [Puij'r.t ,  19611 
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then 
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and 


(-a)  «  exp{A[cos  Sin  r-($-w)  sin  e] 

♦  irt[sin  eln  r  ♦  (♦*n)cos  0]), 

(sin  irt) " 


(A8) 
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Thus,  denoting  the  integrand  of  (Al)  by  I 
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(All) 
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Therefore,  for  In  2’  <  C  (|a|  <1),  and  - *1  / 2  <  a  <  n/2, 
closure  of  the  integral  in  (Al)  In  the  right  half  plane 
dees  not  contribute  to  the  value  of  the  integral.  Thus 

i- 

(2.<)'1  <}  lg  -  (2«i)’'  J  Ig,  (A12) 

-t* 

and  evaluating  the  integral  of  (Al)  by  sunming  the  residues 
at  the  poles  enclosed  by  the  integration  path  of  (Al) ,  we 
obtain 


1  r(c»»i)r(i+«)  1  J  • 


(A13) 


By  writing  out  several  tcr;ns  of  the  series  in  (A13),  we  can 
recognize  that  the  tcra  within  the  brackets  is  just  the 
hypergcor.ictric  function, 
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or  we  then  have  the  representation 

I  -  [r(a)r(b)/r(e)]f’(a,t  ;e;a) .  (A15) 

Similarly,  closure  takes  place  in  the  left  hand 
plane  when  | s  | >  1 .  In  this  ease  evaluating  the  integral  by 
residues  gives  the  following  csult 

(A16) 

T(0‘b) 

The  two  series,  (A15)  and  (A16)  arc  then  representations  in 
different  regions  of  the  integral  function  in  (Al). 
Equivalently,  v:c  may  consider  the  two  series  of  (A13)  and 
(A16)  to  be  the  analytic  continuations  of  one  another,  and 
we  may  write  the  formula  for  this  analytic  continuation  os 

F(a,b;c’,z )  •  f~(f)r lo~-~a)  aF{a,o'a-\\b-a-l\\/ z) 

(A17) 
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V 

The  homogeneous  solutions  Z01,  Z02  arc  given  in  (S3). 

Because  of  tlic  products  of  these  solutions  occur  in  the 
integrals  of  (B3) ,  the  number  of  individual  terms  is  in¬ 
creased  considerably  above  the  six  of  (53).  Hence,  only  a 
generic  term  will  be  evaluated  for  each  of  the  six  terms 
in  (B3) ,  and  will  serve  as  an  example  of  the  integration 
technique.  These  six  terms  arc  given  by 
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APPENDIX  B 

EVALUATION  OF  THE  SCATTERING  INTEGRALS 

In  tnis  appendix,  we  shall  evaluate  a  typical  terra 
which  occurs  in  the  expressions  for  the  reflection  and 
transmission  coefficients,  namely 
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«'i'l*'*,l,(*)Z01C*)Z||2{!c)(!x. 


(Bl) 


This  term  must  be  broken  down  into  the  six  terns  given  in 
(102).  In  that  equation,  we  make  the  substitution 


C 


..mx 


(B2) 


so  that  (Bl )  becomes 

*  I  (1*0 4  z o i  CO z02 

(B3) 

♦ 


o 

.  (B4) 

1 

*  F(a  ,o-y+1;®*p+i  ;1/0<*s 

m 

|  (*,)(Y’2o'1‘2,/m)*(1-5,/V  Ca,o-Y+l,a-B+l;l/C') 

0 

*  F*(a,a-Y*l;o-B*l;l/C,)d5' 

•  1 

|  ib'mi i-o^^(«»p;Y,*«)f,Co-Y+1.3-Y*1;2-Y;Od« 

0 

«i 

J  i- !-&/»)*  (1-5'  /*F*(a,a*Y+l  Jo-B+ljl/C') 

0 

F*(aia-Y+l;«‘B4l;l/5')d^' 

j  5Y-2o-l^/m(1.0Xf(a>a.Y+1;a.&4l;1/0 

*  1 

*  F(a,a-Y  +  l;a-B+l;l/0<*5 

0 

*  F*(a,a-Y  +  I;n-B«>1;1A,)<*S' 

•  m 

|  ^Y  1  ^^m(l-5)4F(a,o-Y+llo*S'fl'*l/^) 

c 

F(a,a-Y  +  l*o_2*l;  1/0^5* 


x 
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The  first  integration  (over  4)  can  be  categorized 
into  two  types: 


Type  I  0  0 


Type  II 


C’  -l  -l 


CBS) 


-1  -1  4 ' 

-•  ?'  — 

B.l  TEE  TYPE  I  IHT EG UAL 


From  Erdilyi  [1953]  , 

P(a,b;o;a )  •  (l-*)“fcF(b,e-a;o;  [-«/(!-*)]) .  (B6) 


The  integral  representation  for  the  hypergcometric  function 
discussed  in  Appendix  A  will  be  used  to  represent  the  two 
hypergcometric  functions  occurring  in  (B5),  i.e., 


F(a,b;o;a) 


r(a)r(b) 


<• 

1  r  (rj+t) 


r(-t) ( -z)tdtf 


(B7) 


x 
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where  the  contour  of  integration  is  taken  so  that  one  leg 
of  the  integration  path  extends  over  the  interval  (-t®,  £®) 
and  the  poles  of  the  functions  r(a+i)  and  lie  to  the 

left  of  the  path  of  integration  and  the  poles  of  the 
function  r(-fc)  lie  to  the  right  of  the  path.  Employing  a 
liberal  exchange  of  integration  order,  we  may  write  (B5)  as 


X(2n i) 


-2 


rCE+tQrCg-a+t'l 

r(e+t') 


r (-t')dt’ 


-t® 


t® 


|  r(-o« 

-i® 

J  c  (l-o  dK. 


(B8) 


The  integration  over  *  is  then  given  by  the  integral 


I  • 


(89) 


From  Gradohtycn  and  Ryshik  [1965] , 


0 

|  xv-1  (l-.x)’viix  *  -«wu’1f(v,p  ;l+n  ;u)  ,  (BiO) 

u 
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and  by  using  (B6) ,  this  becomes 

0 

^l-x)  Vdx  «  -y  1  [u/  (l-u)]p 
V 


(Bll) 


*  F(y , 1+y- v, 1+y , [ -w/ (1-tf) ] ) . 

0 

The  integral  over  the  definite  limits  /  is  just  a  special 
ease  of  (lilO) ,  while  the  integral  over  the  limits  fCj  can 
be  broken  into  two  parts  as  follows 


(B12) 


All  three  integrals  of  Type  I  can  thus  be  treated  in  the 
same  manner.  Our  Type  I  integral  thus  reduces  to  the 
following  double  integral 


I\  -  -K(2.£)-2[C'/(l-C,)]e*1 


X 


r(b+tDr(o-a+tD 

r («♦*') 


i- 


|  mi 


r(-t)[?'/(l-€')l 


(B13) 


(B+l+ht'J^FCBtinn'  ,A  +  n+2-b’B;B+2+t+t' ;  [-5  7 (!-«’) ])«**• 


X 
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Evaluation  of  these  integrations  can  be  carried  out  in 
the  manner  of  Appendix  A .  Following  the  first  integration, 
we  have 


1 ® 
e 


KU.ifU'/ii-e’))8*1 


J  r  (  - « ’ )  K '  /  0-  C )  ]  * ' 


-t® 


X 


m 

i 

l*  o 


r(e+nr(Y-a+n 

(y+D 


[*0/(i*0]1 


(B14) 


*  (B+l  +  tW)'1F(/?+l  +  t,*l,>!  +  0+2-i-S;0+2+fc'  +  l;[-c/(l-O])<ft  • 


By  integrating  again  after  freely  interchanging  operation 
orders,  we  have  the  result 


zl 


'  jj"*1 


X 


?  r(fl*nr(v»»n 

ji0  r(y*i)r(m) 


[•e/U-Ol* 


r  r (2?+n)r(o-a*n) 
*  *  r(e+n)r(l+n) 

n*° 


[-c/d-01n  . 


(BIS) 


B.2  THE  TYPL'  II  1 11  TEG  HAL 


Using  the  formula  of  (B6)  ,  we  may  write 


F(a  ,  &  ;c ;  I/O  ■  (-5/ (l-?)]F[i  ,a-a  ;  (1-5) 


(B16) 
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Again  using  the  integral  representation  of  the  hypcrgcomctr ic 
function,  and  following  a  change  in  integration  order,  the 
Type  II  integral  becomes 


(B17) 


•t" 


“  }  cB*fc+6(l-0 A‘b‘*‘tmt' dt.. 


To  integrate  over  £,  it  is  necessary  to  transform  the 
integral  of  (BIO). 

Let 


*  ■  -1/t, 


(B18) 


and  we  obtain 


tv  w  vdi  ■  ("l)1  v  1F(v,»i;1*p;-m)  .  (B19) 
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Now  let 


S  ■  -1/h» 


(B20) 


to  obtain 


j.  jVw* 1  -  (-1)''V%mfCv,i.;1»i*;1/s)- 


(B2I) 


Again  using  (M).  "'c  obtain  t,,c  rCS“U 


I  ■  c-n^V'ci-s)'* 

1  -1, 


(B22) 


The  remainder  of  the  integration  proceeds  as  in  the  Typo  I 
evaluation  with  the  result 


A  .  .  i  ^ 

*  •  J«0  *■  ^ 


r  (b)w(o*a)r:  (i-C) 

n-0  C<J5n(1)n 


-» 


(B23) 


X 


The  integral  !  1  can  be  broken  up  into  the  two 


parts 


(B24) 


The  integration  which  rcuains  (over  r,')  can  now  be  evaluated 
exactly  as  has  been  done  for  the  4  integration.  The  results 
arc  as  follows 


I 


(-1)*(B+BA) 


(2)- 2) 


m 

i 

lm  o 


(T)t(Dj 


nm  0 


(*)  U) 

n  n 


X 


X 


(A  +  B*2-b-S)k  J  (6*)i(Y*-o*)f  w 

(B+l+l+n+k)(l)k  (y^CD* 

w  *  (d**B*+2-£*-$*)m 

(o*)*(l),  *  (B*B**2*l*n+k*i+j*m)  ( 1 ) 

J  J  IH"  0  m 


(B25) 


2 


-m 


(B+B*)2>t+/t*+B+B*+2 


-  Cg)l(Y-«)l  t  -  CMw(c-a)w 

im  <*W  “V1** 


-  l*(B+2>+B)]fc 

-  J(  (J*„**-*-8-l)7I)~ 


£■0 


Cv*3<Cl)i 


(B26) 


m 

i 

im  o 


(»•)/«* 

— c5*)TuT7J‘2 

j 


I-(B**i**B*)L 

j  r  _ _ _ 2 _ 2  m 

*  (i+j  +  :  +  n+k+n-/1*-/l-B*-B-2)  (1) 

m*o  m 


x 
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In  order  to  nahe  the  expression  for  the  scattering  co¬ 
efficients  as  compact  as  possible,  the  results  of  these 
integrations  have  been  defined  in  terns  of  certain  functions 
as  follows 


<M o.B.Y.a.i.e.Fj  -  (-l)**2‘<fl+2>  l  2'1 

1*0  <^(1)* 


-  (*>)(<?-<*)  „  “  (B+A+2‘b-&) 

p  _  n  _  jt  ~‘tt  r  _ _ K  K 

*  o  (*)„(!)„  o  (*♦!♦*♦«♦*)  (Dk 


(B  27 ) 


G9(a,e,y,a,b,c,P)  *  (-1)  2 


B0A  +  B+2  r  C6) ^  (Y-a)j  - 


lie  (yMJ>Z 


(B28) 


-  Cb\l(c-c)r  -  _  k 

nl 0  k»0  U**«-A-B- !>(!)* 


^  —  —  »  «.  •  r pi p]  •  f ) 

#1  fa  .  6  ,Y  .a  »h  .e  »a  .  B  .  Y  .«  tb  ■  (*1)  2 


ti'-miMi  X 

p  (F+4+2-&-B)*  -* 

*  *£0  (*♦!♦*♦"♦*)(!)* 

-  (S*)j(7*-a*)t  ^  vJ 

i~o  Cr*)t.(l  )i  j«  o 

•  (I**F**2-F*-?*)„, 

T  - - - - - -  2~m 

rjBO  (F*F*+2+Z+>:+fc+i+-**)(l)w 


(B29) 


X 


■  r — i| 
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®  2  »  ®  * 


*  I 
<-0 


y  ,a,b,o,a,J,y,a,b,o,B,B] 


^iy(B*B)2A*A*B*B*2 


(r)z(Dl 


•  (i)  (e-a) 

I 


n-fl  f®)  <*> 

”  o  n  n 


-  2*n 


[-(B*fc+e)]fc 


"  J(  (1)^ 


($*) . Cy^-o*)  . 


j  j 

f**)/11* 


(B30) 


•i 

i 

mm  o 


[•(B*+b*+B*)]m 

{i+j+k+Um+n-A*-A‘B*‘B-2){l)m 


X 


97 


Qla  ,b  ,c  ,d  ,e  ,f  ,a  ,b  ,c  ,d  ,e  ,f  ,U }  ,U  2  ,U  ^,U 
“  Hl{a  ,b  ,c  ,d  ,e  ,f ,  (a  ,b  ,o  ,d  ,e  ,/)  ,d-a*y-l*l/m ,  (d-a*y-  1-b/m)*] 

*  G l[a  ,b  ,c  ,d  ,c  ,f  ,d-a*y‘l]G [a  ,b  ,c  ,d ,e  ,f  ,d~a*y  1  +b/m) 

*  H  ^la ,b ,c  ,d  ,e  ,f , (a ,b ,c ,d ,e ,f)* ,  (d-a+y-l^b/m)* ,d-a+y-l+b/m] 

♦  2Rc  (-1)  f‘a<72  (1 ,4  ,S,1 ,4 ,5,B3) 

♦  t/lt/4(-l)'8'^2(l,4, 5,2,6,7,54) 

♦  i/253(-l)‘e‘d(72(2,6,7,l,4,5,54) 

*  ♦  f/2f/4(-l)a'<i'2B(72(2,6,7,2,6,7,S5)] 

*  (a  ,b  ,o  ,d  ,o  ,f.,d-a+y-l+b/m)  } 

♦  C-l)'<f*e,C2(l  ,4  ,S,1 ,4  ,S,B6J 

♦  ®C2  (1,4,S,2,6,7  ,«7J  (B31) 

♦  1'2»,(-1)'‘,'S<!j(2.6.7,1,4,5,«7) 

♦  BjI',(-l)‘,,'a‘JSC2(2.«.7.2,6.7,S8)) 

«  ((/,(/,  (-l)'0-^  (1,4, 5, 1,4, 5. S3) 

♦  t/lt/fc(-l)*#'<iff2(1.4,S,2.6,7p,s4J 

♦  iy'3(-l)'|!'*«2(2,6.7,1.4,S.»4) 

♦  !/2B4(-l)'°’rf‘2Sff2(2, 6,7.2,6,7,85))* 

•  I|t'1l2|!/3|2(-l)'iRc(i*o) 

»  ff2 [1, 4, S, 1,4, 5, (1,4, 5, 1,4, 5)* ,B6 ,B3*] 

*  |B1l2I'3I/4*(-l)'<i'<i!'*0-e* 

*  H2  [1,4, 5, 1,4, 5,  (1,4, 5, 2, 6, 7)*, 56, 54*] 

♦  I  ^3  I  2  C*  0* 

*  n2 [1,4, 5, 1,4, 5, (2, 6, 7, 2, 6, 7)*, 56, 54*) 

♦  L'xl'sU2*L\*  (-l)'2  2*-2Re(a+tl) 

#2  [1,4, 5, 1,4, 5, (2, 6, 7 ,2,6,7)*  ,56 ,55*  J 


X 
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♦  \ux\2uku2*(- 

♦  UM2IM2(- 

« 

♦  oJffl»|o,l2(- 

♦  VjV'l.**' 

♦  I»2IjI»3I2(- 

♦  U2UliUl*U2*(‘ 

♦  02VKI*(-: 

♦  I^2  |2£/4^3*(- 

♦  |tf2l2KI2(- 


•  a* - d-d* 

x  fl2 [1,4, 5, 2, 6, 7, (1,4, 5, 1,4, 5)*, 87, 83*] 
-2Re(g+£f) 


x  #2 [1  ,4, 5, 2, 6, 7,  (1  ,4, 5, 2, 6, 7)*. ,87, 84*] 
-2Re($+d) 

x  H2 [1  ,4, 5, 2, 6, 7, (2, 6, 7, 1,4, 5)*, 87, 84*] 
a* -£?-<?* -0- 2  C* 

x  82 [1  ,4, 5, 2, 6, 7, (2, 6, 7, 2, 6, 7)*, 87, 85*] 
~  a* - 8~  d-d* 

x  82 [2, 6, 7, 1,4, 5, (1,4, 5, 1,4, 5)*, 87, 83*] 
-2Re(<i+e) 

«  ff2 [2, 6, 7, 1,4, 5, (1,4, 5,2, 6, 7)*, 87, 84*] 
-2Re ( d+e ) 


x  #2  [2, 6, 7, 1,4, 5, (2,6,7,]  ,4 , 5)*  ,87 ,34*] 
•dmd*+a*- 8-2B* 

x  82[ 2, 6, 7, 1,4, 5, (2, 6, 7, 2, 6, 7)*, 87, 85*] 
a’a*  ~dmd* ’ 28 

x  *2 [2, 6, 7, 2, 6,7, (1  ,4, 5, 1  ,4, 5)*, 88, 83*] 
a-d'd* • 28- 8* 

x  82[ 2, 6, 7, 2, 6, 7, (1,4, 5, 2, 6, 7)*, 88, 84*] 
a- d-d* -  28- 8* 

x  82[ 2, 6, 7, 2, 6, 7, (2, 6, 7, 1,4, 5)*, 88, 84*] 
2Re(a-<f-  B) 


82[2,6,7,2,6,7,(2,6,7,2,6,7)*,B8,35*]> 


X 
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♦ 


Bb  +"*  Bl 

Bl  BA  • 

B  8  B  5  4  i 


ami  with  the  following  abbreviated  notation 


Bl  *  y*l *b/n 
B 2  e  yl-v/r: 

F3  c  yl-2<jt-b/n 
B4  *  Y*l*a*S*b/w 
B5  B  y*1*23 -l/n 
Bb  *  Y*l*2a +l/n 
Bl  *  Y*l*o-fl *b/n 
B8  *  y*1*2 $+b/n 
B9  *  bl n 
BIO  »  * b/r, 

Bll  *  1*y *b/m 
B12  •  1*y -b/m 


1  ■  a 

2-6 

3  -  Y 

4  *  1  +  Ci-y 

5  *  l+ct*B 

6  ■  1  +  B* y 

7  *  1*  B  -  a 

8  -  2- Y . 


The  synbol  {•*-*}  represents  an  interchange  in  the  immediately 
preceding  braces. 
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APPENDIX  C 

SC  ATT  F  RING  COEFFICIENTS  OF  T!»E  TURBULENT  EPSTEIN  LAYER 

c.i  Tin:  TUFEumniT  half  space 

From  the  theory  in  Section  6.1.3,  the  incoherent 
and  coherent  reflection  coefficients  are  yiven  ny 


WT  «  c24fi4(l+n2)*  1  [l+2a'*(a'*+e'2)]‘2 

i  p  o 

*  [(fc0a)‘4  +  S(a'2-a,2)/(fe0a)2  ♦  16(o'2+e'2)] 

(Cl) 

x  {[(l/26’fc0a)-l]/(fe0a) 

♦  4 (a ' 2+B 1 2)  ♦  2(a,2*3,2)/6,^0a} 

We  -  |  (l-X)/(l«\Jr)  |  2 

♦  2Re|c2 [ (x*-l)/ (y* • 1)3 [n2  (l-intf)/ (l+n2) ] 2 
«  [2/  (1+X)  ]  *  (4tf)  ”  1  [  (2.Y)2+  (fcQa)  ”2]  ~ 1 

X  { [(*-l)/(X+l)]  (t/2fe0aA-lJ  (C 2 ) 

•i2k0aX-l  ♦  IC2Y)”1  (4A’2fe0a  ♦  2/k0a  *  i 2^) } ) 
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where 


a’  -  Red-;;) 


8'  -  Imu-;;} 


(C3) 


f*  *7*  'T  f  *  7?  f '  * r  r  *  *  *  *  C  V ' ' ' ' '  M  f  .*•  /,•**»  r»p  .•  "f  rc»  r^»i  r  4  V  •  •  r*  r* 

<•  •  £  •  ••  .  ij i'. ,  ,  »  ^ ^  h%*  i)  +  i\r.  *  u  J.  i  1 0 LtA  Y V«*i£ 


Fre;  the  theory  in  Section  5.2,  anil  u s in^i  the 
notntio*'  introduced  in  \npcn<Ji;-;  F,  the  incoherent  reflection 
coefficient  is  i;i  ven  hv 


-  c  h  o  /  2rr )  2  I  // 1  2 1 2  { 

«[l,2,3,l,2,3,(l,2,3,l,2,3)*,Blffl2,flltB2] 

♦  4|V*J2 

*  C[li“chh^)h(l»2|J|4|6,8)*,flj,B2,Bj, B 4 ] 


♦  |P2/\I4 


(C4) 


(J[  4 ,6,8 ,4,6,8  ,  (4 ,6,8,4 ,6,8)*  ,P3,fl4  ,£3  ,Bj 


X 
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2Re(-2(-l)iS(i!2/Sli) 

*  C[1»2,5,4,6,8,(1,2,3,1,2,3)*,Bi,B2»53»54] 

♦  C-l)i25(B2/B4)  2 

*  G[4,6,8,4,6,8,(1,2,3,1,2,3)*,B3,B4,£>3,B4] 

-  2(-l)i5|B2/B4|2(B2/B4) 

x  C[4,6,8,4,6,8,(l  ,2,3,4 ,6,8)*  ,B3  ,B4  ,B3  ,/?„]  )j 

Likewise,  tl>c  incoherent  transmission  coefficient  is 

Vx  »  (m7;2/2) 2  | ’//fc1B4  |  2c2| 

C[1,2,3,4,6,8,(1,2,3,4,6.S)*,B1,B2,B3,B41 

♦  I  2?2/b4  | 2  CCS) 

x  g(4,6,3,4,6,8,(4,6,8,4,6,8)*,B3,B4,B3,B4] 

♦  2Re{(-l)*5(B2/B4) 

x  C[4,6,8,4,6,8,(l,2,3,4,6,8)*,B3,B4,B3,B4]}j 


I 


4 


* 


(ft  p/ ftpo) 
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109 


w,/* 
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wr/f 


115 


(a*  -  .i) 


I04 

Ticir.vr 


w,  /  5 


w ,  A 
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